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V.  BIOTA 


This  chapter  is  the  second  of  three  chapters  which  discuss 
various  aspects  of  the  inpact  assessment  methodologies  developed 
in  Phase  I  of  the  Biota  Assessment.  While  the  previous  chapter 
focused  on  defining  baseline,  this  chapter  discusses  the  selec¬ 
tion  of  study  species  and  the  distributional  mapping  of  those 
study  species.  Chapter  VI,  which  follows,  works  through  the 
development  of  conceptual  and  mathematical  ecological  models. 

One  of  the  first  steps  in  selection  of  study  species  was  to 
delineate  the  various  organisms  by  major  functional  groups. 

These  have  already  been  discussed  in  Chapter  III  and  elsewhere, 
but  are  defined  here  in  terms  of  their  relationship  to  selection 
of  study  species. 

A.  MAJOR  GROUPS 

To  simplify  our  ecosystem  analysis  and  study  species  selection, 
we  have  grouped  the  living  components  of  the  Chesapeake  Bay  eco¬ 
system  into  seven  cagegories  which  reflect  both  function  and 
habitat  as  follows. 

Phytoplankton:  These  are  microscopic,  unusally  single-celled 
plants  which  represent  several  divisions  of  algae.  Functionally, 
the  group  comprises  both  net-  and  nanno-plankton ,  the  latter 
being  species  less  than  10  ym  diameter  (Van  Valkenburg  and  Flemer 
1974) .  Phytoplankton  are  pelagic,  and  are  moved  about  by  actions 
of  currents  and  tide.  Some  workers  further  distinguish  ultra¬ 
plankton,  species  less  than  2  to  3  ym.  This  category  would  in¬ 
clude  most  planktonic  bacteria,  which  are  heterotrophs  and,  as 
a  group,  not  well-known  in  the  Bay. 

Submerged  Aquatic  Vegetation :  These  are  plants  -  usually  rooted  - 
which  live  submerged  below  the  water's  surface.  Submerged 
aquatics  in  Chesapeake  Bay  are  chiefly  angiosperms  (seed  plants) , 
although  some  species  (e.g.  Nitella)  are  macroalgae. 

Emergent  Aquatic  Vegetation:  These  are  plants  which  grow  par¬ 
tially  submerged,  regularly  or  occasionally  flooded,  or  in 


wet  soils.  They  make  up  the  bulk  of  vegetation  in  marshes  and 
other  wetlands. 

Zooplankton:  These  are  usually  small,  sometime  microscopic, 
animals  from  several  phyla.  The  group  is  composed  of  holo- 
p lank ton,  species  which  are  planktonic  throughout  life,  and 
meroplankton ,  species  which  spend  only  part  of  their  life  cycle 
in  the  plankton. 

Benthos :  The  benthos  is  comprised  of  organisms,  mainly  inver¬ 
tebrates,  which  live  associated  with  the  substrate.  These  may 
be  epi fauna  -  species  which  live  attached  on  or  above  the  bottom  - 
or  infauna,  species  which  burrow  into  the  substrate.  Some  spe¬ 
cies,  such  as  crabs,  are  benthic  oriented,  but  are  motile  or 
vagile,  capable  of  considerable  swimming.  Other  species  are 
benthic  only  at  some  stage  of  their  life  cycle,  such  as  the  sea 
nettle  (Chrysaora  quinguecirrha) .  The  benthos  is  often  divided 
into  macro  -  ( greater  than  0 . 5  mm) ,  meio  -  (0.5  -  0.1  mm) ,  and 
micro  -  benthos  (less  than  0.1  mm)  (Coull  1973).  Of  these,  only 
the  macro-benthos  is  well-known  in  Chesapeake  Bay  (Lippson  et 
al.  1979) . 

Fish:  Fish  make  up  the  bulk  of  the  Chesapeake  Bay  nekton-species, 
which  exhibit  well-developed  powers  of  movement.  Dermersal 
fishes  are  those,  such  as  spot  (Lieostomus  xanthurus)  which  are 
associated  with  the  bottom  and  feed  chiefly  on  benthic  organisms. 
Pelagic  species  feed  in  the  water  column,  chiefly  on  fish  and 
macroinvertebrates . 

Wildlife  (Amphibians,  Reptiles,  Birds,  and  Mammals) :  Numerous 
waterfowland  shorebirds,  and  some  water-oriented  amphibians, 
reptiles,  and  mammals  use  the  estuary  and  adjacent  wetlands  for 
food,  shelter,  and  breeding  areas.  The  dependence  on  the  estu¬ 
ary  varies  from  species  to  species.  Many  of  the  birds  are  migra¬ 
tory  and  use  the  estuarine  resources  seasonally. 

B.  DISTRIBUTION  OF  ESTUARINE  ORGANISMS 

It  is  generally  acknowledged  that  the  distributions  of  estuarine 
organisms  do  not  have  the  sharp  boundaries  implied  by  the  Venice 


System  (Wolff  1973,  Boesch  1977) .  Rather,  each  species  popula¬ 
tion  itself  exhibits  a  continuum  or  cline  along  the  estuarine 
complex-gradient;  the  point  of  maximum  abundance  reflects  both 
innate  physiological  tolerances  and  effects  of  biotic  and  abio¬ 
tic  features  of  the  environment  (Whittaker  1970) .  Apparent 
discontinuities  may  result  from  substrate  changes,  effects  of 
competition,  or  the  overlapping  of  similar  distribu¬ 
tions  along  the  estuarine  gradient  (Boesch  1971,  Wolff  1973) . 

Nevertheless,  there  is  evidence  that  some  correlation  exists 
between  the  Venice  System  boundaries  and  estuarine  biotic 
zonation.  Dahl  (1956)  reviewed  numerous  studies  from  brackish 
environments  and  identified  three  zones  of  more  rapid  biotic 
change,  at  0.1  -  0.5  7^  ,  5.0  -  8.0°/^  ,  and  15.0  -  20. 0°/^^ 
Khlebovich  (1969)  characterizes  the  salinity  range  between  5.0 
-  8 . 0 °/  as  the  "critical  salinity";  below  this  range,  hyper¬ 
osmotic  regulation  is  required  to  prevent  internal  salinity  from 
dropping  below  4.0  -  5.0  %o  (at  which  point  serious  tissue 
damage  occurs).  Kinne  (1963,  1964)  divides  animals  into  four 
groups  based  on  their  osmoregulatory  abilities,  although  there 
is  variation  within  each  group.  These  groups  are: 

•  stenohaline  osmocon formers,  with  little  or 
no  capacity  for  regulation, 

•  euryhaline  osmoregulators ,  which  can  regulate 
in  water  of  reduced  salinity,  but  not  fresh 
water, 

•  holeuryhaline  osmoregulators,  which  can 
regulate  from  fresh  to  full  oceanic  salini¬ 
ties  ,  and 

•  oligohaline  osmoregulators,  which  can  regulate 
only  in  fresh  water  and  very  low  salinities, 
and  maintain  blood  hyperosmotic  to  the  external 
medium. 

The  estuarine  biota  is  thus  composed  of  organisms  which  show 
varying  amount  of  adaption  to  conditions  within  that  environ¬ 
ment  . 

Day  (1951,  1967)  first  recognized  five  major  components  in 
South  African  estuaries:  fresh  water,  true  estuarine,  eury¬ 
haline  marine,  stenohaline  marine,  and  migratory.  Carriker 
(1967)  relates  these  to  the  Venice  System  boundaries,  and  char¬ 
acterized  each  group  by  general  salinity  tolerances,  origin, 
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and  reproductive  requirements.  Boesch  (1971,  1977) further 
refines  this  classification  and  recognizes  the  following 
major  groups: 

Stenohaline  Marine:  These  organisms  are  characteristic 
of  euhaline  environments,  but  occasionally  penetrate  the 
estuary  to  25.0  Within  Chesapeake  Bay  they  are  re- 

stricted  to  the  polyhaline  zone.  Example:  the  cladooeran 
Evadne  tergistina  (Bryan  1977) . 

Euryhaline  Marine:  The  organisms  extend  from  the  euhaline 
zone  into  the  estuary,  sometimes  to  relatively  low  salin¬ 
ities  ( 15. Oppt) .  This  group  comprises  the  larges  part  of 
the  estuarine  biota  in  the  polyhaline  zone.  Some  species 
depend  on  recruitment  from  the  marine  environment,  but 
many  have  viable  reproducing  populations  within  the  estu¬ 
ary.  Example:  the  venerid  clam  Mercenaria  Mercenaria . 

Euryhaline  Opportuniate :  These  are  species  which  are  found 
from  the  euhaline  zone  (often)  to  the  oligohaline  zone. 
However,  they  are  generally  most  numerous  in  the  low  poly¬ 
haline  and  mesohaline  reaches  of  the  estuary.  Their 
reproductive  strategy  allows  them  to  colonize  rapidly 
disturbed  or  stressed  habitats,  as  well  as  salinity 
regimes  where  less  eurytopic  species  are  a  competive  dis¬ 
advantage  (Carriker  1967,  Boesch  1971,  1977).  Example: 
the  polychaete  Heteromastus  filiformis . 

True  Eetuarine  or  Estuarine  Endemics :  These  species  are 
restricted  to  the  estuary,  either  by  physiological  limit¬ 
ation  at  some  part  of  their  life  cycle  or  by  competition 
with  marine  species  in  offshore  environments  (Kinne  1966, 
Day  1967,  Jefferies  1967,  Diaz  1977).  Some  may  be  re¬ 
stricted  to  the  estuary,  however ,  by  substrate  or  circula¬ 
tion  requirements  (Day  1967,  Boesch  1971).  In  general, 
they  are  dominant  members  of  the  biota  below  15. Oppt. 
Example:  the  mactrid  clam  Rangia  cuneata. 

Tidal  Fresh  Water  and  Oligohaline :  These  species  are 

derived  from  fresh  water  forms,  but  are  able  to  tolerate 
varying  amounts  of  salt.  Most  are  restricted  to  salinities 
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below  0.1  %,,  but  some  persist  to  5.0%,  or  higher 
(Carriker  1967,  Diaz  1977).  In  general,  the  biota  of 
the  tidal  fresh  water  and  oligohaline  zones  consists  of 
a  mixture  of  a  few  eurytopic  fresh  water  species  and 
estuarine  endemics  (Diaz  1977) .  Example:  the  tubificid 
worm  Limnodrilus  hof fmeisteri. 

The  first  four  of  these  categories  of  organisms  are  of  marine 
origin,  and  have  adapted  to  the  estuary  to  varying  degrees. 

Boesch  (1977)  has  developed  a  general  model  of  estuarine  zona- 
tion,  based  on  distribution  of  benthic  invertebrates  along  the 
Chesapeake  bay  -  York  River  gradient  (see  Figure  III-14)  .  He 
found  faunal  changes  to  be  gradual,  with  zones  of  somewhat 
accelerated  change  in  the  3.0  -  8.0%>  and  the  15.0  -  20.0%> 
range.  There  is  a  gradual  decrease  in  species  richness  from 
the  polyhaline  to  the  oligohaline  zone.  Diaz  (1977)  found  the 
lowest  diversity  of  benthic  macro- invertebrates  in  the  oligo¬ 
haline  and  tidal  freshwater  zones  of  the  James  River,  with  an 
increase  in  diversity  in  non-tidal  freshwater  areas.  The  oligo- 
haline/tidal  freshwater  zones  were  high  energy  environments, 
which  were  turbid,  with  a  low  diversity  of  benthic  habitats. 

This  general  model  of  distribution  does  not  hold  as  well  for 
some  other  groups.  Planktonic  organisms  can  be  carried  into 
areas  outside  their  normal  salinity  range  by  actions  of  currents 
and  tides.  Typically,  "fresh-water"  phyto-  and  zooplankton 
are  found  throughout  the  oligohaline  zone,  and  often  carried 
into  the  low  mesohaline  by  normal  river flow  (Goodwyn  1970) . 

Species  typically  found  in  euhaline  or  polyhaline  environments 
may  be  carried  into  the  estuary  in  the  inflowing  deep  layers, 
later  being  admixed  into  upper  layers  of  lower  salinity  zones 
(Burrell  1972) .  Fish  are  able  to  move  with  relative  freedom 
within  their  range  of  physiological  tolerances.  Euryhaline 
marine  species  such  as  menhaden  will  be  found  much  further 
upstream  than  is  shown  in  Figure  III-14 , whereas  estuarine  endemics 


such  as  killifishare  found  further  downstream.  Larval  fishes, 
however,  ususally  have  more  restricted  ranges  due  to  their 
limited  osmoregulatory  abilities  (Polgar  et^  al.1976)  . 

It  is  a  generality  that  species  richness  and  diversity  is  low¬ 
est  in  the  oligohaline  zone,  although  this  is  not  equally  true 
for  all  groups.  Over  540  species  of  phytoplankton  were  re¬ 
corded  from  a  non-tidal  freshwater  site  on  the  Potomac  (ANSP 
1972)  ,  while  only  160  taxa  were  recorded  from  the  tidal  fresh- 
water/oligohaline  zone  by  Dahlberg  (1973)  .  Even  fewer  species 
(around  80)  were  found  at  Maryland  Point  on  the  Potomac,  a 
variable  oligohaline/low  mesohaline  area  in  the  estuarine  tran¬ 
sition  zone  (Mountford  1971) .  While  Morse  (1947)  identified 
over  200  from  the  mesohaline  region  at  the  mouth  of  the  Patu¬ 
xent,  total  species  number  obtained  from  several  surveys  in  the 
Bay's  polyhaline  zone  add  up  to  over  205  taxa  (Patten  et  al. 
1963,  Marshall  1967,  Mackiernan  1968). 

The  number  of  zooplankton  species  is  generally  high  in  fresh¬ 
water  areas,  in  part  due  to  the  diversity  of  rotifers  and  clad- 
ocerans  in  this  environment  (Dahlberg  et  a_l.  1973)  .  At  Douglas 
Point,  116  species  were  recorded,  of  which  33  were  rotifers  and 
25  cladocerans.  Similarly,  much  of  the  increased  species  rich¬ 
ness  of  macroinvertebrates  in  freshwater  environments  is  due  to 
the  high  proportion  of  insect  larvae;  at  the  non-tidal  Dickerson 
site,  387  out  of  452  species  were  insects,  while  38  species  out 
of  76  were  insects  at  Douglas  Point  (Dahlberg  et  al^.  1973)  .  In 
the  oligohaline  zone  the  majority  of  these  salt  intolerant  forms 
dissappear . 

Submerged  and  emergent  vegetation  have  their  highest  diversity 
and  species  richness  in  the  lower  salinity  areas,  although  for 
marsh  species  the  period  and  extent  of  tidal  inundation  is 
more  important  than  salinity  in  determining  distribution  (Boon 
et  al.  1977,  Orth  et  al .  1979).  There  are  relatively 
few  submerged  higher  plant  species  found  in  polyhaline  and 
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euryhaline  environments;  in  Chesapeake  Bay  only  Zostera 
marina  and  Rupp i a  maritima  occupy  this  zone  (Orth  et  al.  1979) . 

It  is  apparent  that  the  observed  biotic  zonation  of  the  estuary 
reflects  changes  in  importance  of  the  major  components  along  the 
ecocline.  A  generality  which  can  be  made  is  that  the  organisms 
of  marine  origin  are  limited  up-estuary  by  salinity  constraints, 
and  down-estuary  by  biological  interactions,  although  in  some 
cases  salinity  again  becomes  a  limiting  factor  (Larsen  1974, 

Boesch  1977,  Diaz  1977,  Heinle  et  al.  1978). 

Estuarine  endemics  and  euryhaline  species  often  show  greater 
range  of  salinity  tolerance  in  the  laboratory  than  is  realized 
in  the  environment  (Cain  1974,  Castagna  and  Chanley  1973).  They 
may  be  restricted  to  lower  salinities  by  predation  (Kinne  1966, 
Larsen  1974)  or  by  narrow  salinity  tolerances  of  a  particular 
life  stage  (Cain  1972,  Tagatz  1968).  Such  estuarine  species 
need  to  have  mechanisms  for  retention  of  their  larvae  and  juve¬ 
niles  within  the  estuary.  Some  brood  their  young,  or  have  non- 
planktonic  larval  stages  (eg.  the  oyster  drill  Urosalpinx  cinerea) . 
This,  however,  reduces  their  ability  to  colonize  new  habitat  or 
former  habitat  from  which  they  have  been  eliminated  (D.  Haven, 
personal  communication) .  Many  species  with  planktonic  larvae 
have  evolved  behavioral  mechanisms  which  take  advantage  of  estu¬ 
arine  circulation  to  enter  or  remain  within  the  estuary,  eg. 
oyster  larvae,  blue  crab  megalopes,  larval  fish  (Harrison  et  al. 
1967,  Wood  and  Hargis  1971,  Sandifer  1973,  1975). 


C.  INITIAL  STUDY  SPECIES  SELECTION 

The  Chesapeake  Bay  Low  Flow  Study  requires  that  "approximately 
fifty"  study  species  will  be  used  to  assess  the  biotic  effects 
of  low  fresh  water  inflows.  Selection  of  these  study  species, 
and  justification  of  their  choice,  is  an  extensive  task.  It  is 
estimated  that  over  2650  species  of  plants  and  animals  live  within 
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the  Chesapeake  Bay  (McErlean  et  al.  1972) .  These  include  com¬ 
mercially  and  recreationally  important  forms,  as  well  as  numer¬ 
ous  species  which  represent  major  links  in  the  estuarine  trophic 
structure. 


In  order  to  assess  the  effects  of  perturbations  on  the  tidal 
ecosystem,  study  species  should  represent  taxa  and  functional 
groups  from  the  major  Bay  habitats,  salinity  zones,  and  biotic 
subcomponents.  Swartz  (1972)  emphasized  that  selection  of  test 
organisms  should  be  based  on  their  relative  vulnerability  to 
change  and  stress,  ecological  significance,  distribution  within 
the  estuary,  phylogenetic  representation,  and  economic  signifi¬ 
cance.  To  this  should  be  added  the  practical  consideration  of 
data  availability  for  each  species. 

With  Swartz's  criteria  in  mind,  a  preliminary  list  of  167 
potential  study  species  was  generated  and  circulated  at  the 
November  1979  seminar  (see  Chapter  II) .  This  was  based  on  a 
variety  of  sources  , including  the  following  published  list  of 
major  Chesapeake  Bay  species: 

•  Chesapeake  Bay  Existing  Conditions  Report (USACE  1973) :  A 
list  of  110  species  recommended  for  bioassay  or  con¬ 
dition  indices,  in  order  to  assess  effects  of  environ¬ 
mental  stress  (Swartz  1972) . 

•  Chesapeake  Bay  Future  Conditions  Report (USACE  1977)  :  A 
list  of  126  important  species  and  genera,  based  on  a 
survey  of  Bay  researchers.  Species  were  included  on 
this  list  on  the  basis  of  15  criteria,  including  im¬ 
portance  to  trophic  structure,  and  distribution.  An 
attempt  was  made  to  include  species  representing  as 
many  Chesapeake  Bay  habitats  as  possible. 

•  Maryland  Department  of  Natural  Resources  list  of  44 
representative  species  from  tidal  and  non-tidal  waters. 
These  species  are  to  be  used  in  studies  assessing  impact 
of  discharges  into  natural  waters. 

In  addition,  numerous  reports,  papers,  and  data  sets  were  con¬ 
sulted  to  identify  major  species,  and  their  general  distribu¬ 
tion  in  regard  to  salinity.  The  following  were  major  sources 
used  to  generate  the  initial  study  species  list. 
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Phytoplankton:  Patten  et  al.  1965;  Marshall  1966,  lit-  ; 
Mackiernan  1968?  Mulford  1972;  Dahlberg  et  al.  1973;  Van 
Valkenberg  and  Flemer  1974;  Seliger  et  al.  1975;  Lear  and 
Smith  1976;  Mountford  1977;  Van  Valkenburg  et  al.  1978. 

Submerged  Aquatic  Vegetation:  Orth  1975,  Stevenson  and 
Confer  1978;  Anderson  1979  unpublished  data;  Orth  et  al 
1979;  Munro  1979;  Migratory  Bird  Habitat  Reseach  Labora¬ 
tory,  unpublished  data. 

Emergent  Aquatic  Vegetation:  Maryland  Wetlands  Survey,  DNR, 
1967  -  1968;  Keefe  1973;  Metzgar  1973;  Virginia  State 
Wetlands  Survey  Series,  VIMS,  1973  -  1978;  Boon  et  al. 

1977. 

Zooplankton:  Heinle  1966,  1969;  Herman  et  al.  1968; 

Bosch  and  Taylor  1968,  1973;  Goodwyn  1970;  Burrell  &  Van 
Engle  1976;  Dahlberg  et  al.  1973;  Heinle  et  al.  1975; 

Sage  et  al.  1976;  Bryan  1977;  Grant  1977;  Sage  and  Olsen 
1977;  Jacobs  1978;  Grant  and  Olney  1979;  Lippson  et  al. 

1979. 

Benthic  Organisms :  Corey  1967;  Pfitzenmeyer  1961,  1970, 
1973,  1975,  1976;  Boesch  1971,  1972,  1973,  1977;  Wass  et 
al.  1972;  Hamilton  and  LaPlante  1972;  Davies  1972;  Orth 
1973;  Larsen  1974;  Diaz  1977;  Mountford  et  al  1977;  Virn- 
stein  1977,  1979;  Haven  et  al.  1977,  1979;  Lippson  et  al. 
1979?  Reinharz,  Bricker  &  O'Connell  1979;  Cory  and  Dresler 

1980,  unpublished  data. 

Fish:  Hildebrand  and  Schroeder  1928;  Smith  et  al. 1966; 
Ritchie  1970;  Douglas  and  Stroud  1971;  National  Marine 
Fisheries  Service  Fishery  Statistics  of  the  U.S.  1976  “ 
1978?  Scott  and  Boone  1973;  Lippson  and  Moran  1974?  W.R. 
Carter,  unpublished  data?  NMFS  Current  Fishery  Statistics 
1975  -  1978. 


Birds  and  Mammals:  Dozier  1947;  Stewart  1962;  Willner 
et  al.  1975;  Perry  and  Uhler  1976;  Maryland  Dept,  of  Natural 
Resources,  Midwinter  Waterfowl  Surveys  1975  -  1980; 

Virginia  Fish  and  Game, Midwinter  Waterfowl  Surveys 
1975  -  1980;  Rawls,  unpublished  M.S. 

Eight  basic  selection  criteria  were  used  to  choose  species  for 
the  original  list.  These  were: 

Sensitivity  to  Salinity:  Salinity  tolerances  for  candidate  species  (see 
App.  B)  were  evaluated  from  several  sources:  laboratory  studies; 
field  studies,  and  extrapolation  from  field  collection  data. 

Although  the  majority  of  estuarine  organisms  tend  to  be  rather 
euryhaline,  many  exhibit  greater  stenotopy  at  certain  stages 
of  their  life  cycle;  eg.  Rangia  cuneata,  the  larvae  of  which 
require  salinities  between  2  -  10  to  survive  (Cain  1974, 

Hopkins  et  al.  1973) .  Laboratory  studies  commonly  demonstrate 
a  wider  range  of  salinity  tolerance  than  the  species  exhibits 
under  field  conditions  (Castagna  and  Chanley  1973) .  This  may 
reflect  interaction  of  salinity  with  some  other  factor  such  as 
temperature  or  substrate,  range  restriction  due  to  predation  or 
competition,  or  a  stenotopic  life  stage  (Kinne  1966,  Van  Engel 
1958,  Boesch  1977). 

Sensitivity  to  Other  Factors:  Chief  among  these  would  be  fac¬ 
tors  which  themselves  might  be  affected  by  salinity  or  low 
fresh  water  inflows: 

Circulation :  A  partially  mixed,  moderately  stratified 
estuary  such  as  Chesapeake  Bay  is  characterized  by  a 
net  seaward  flow  of  lower  salinity  upper  layers  and  a 
net  upstream  flow  of  higher  salinity  deep  water  (Prit- 
ard  1956,  1967) .  In  general,  the  outflow  at  the  surface 
is  the  driving  force  for  the  rate  of  inflow  of  higher 
salinity  bottom  water  (Pritchard  1967,  Tyler  and  Seliger 
1978) .  Many  organisms  use  the  upstream  movement  of 
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water  at  depth  to  transport  themselves  into  and  maintain 
themselves  within  the  estuary  (Haven  1957,  Harrison  et  al. 
1967,  Wood  and  Hargis  1971,  Sandifer  1973,  1975) .  Reduced 
freshwater  inflow  could  alter  the  rate  of  transport,  and 
allow  breakdown  of  density  stratification,  particularly 
in  the  tributary  rivers  (Schubel  1972) .  This  could  have 
effects  beyond  simple  salinity  tolerances,  if  important 
commercial  species  such  as  blue  crabs,  croaker  and  spot 
are  prevented  from  reaching  their  upstream  nursery  areas, 
or  oysters  their  upstream  beds. 

Temperature:  The  synergistic  effects  of  temperature  and 
salinity  have  been  described  by  Kinne  (1963,  1964)  and 
others.  Temperature  stress  can  narrow  the  salinity  toler¬ 
ance  zone  for  many  organisms,  and  vice  versa.  For  example, 
in  lower  salinities,  the  copepod  Arcartia  tonsa  has  a 
competitive  advantage  over  the  congeneric  clausi  at 
temperatures  from  11  to  18°  C,  as  it  is  less  affected  by 
the  salinity  stress  (Jeffries  1962) .  Chesapeafc 3  Bay 
represents  the  maximum  northward  range  extension  of 
several  southern  species  such  as  Rangia  cuneata,  ana  the 
southernmost  extension  of  others  such  as  Mya  arenaria. 
Adverse  salinities  during  cold  or  warm  periods,  respec¬ 
tively,  could  have  a  more  severe  effect  than  that  produced 
by  salinity  alone. 

Food:  Some  species,  themselves  euryhaline,  are  dependent 
on  a  more  stenotopic  food  source.  For  example,  the  red¬ 
head  (Ay thy a  americana)  feeds  extensively  on  Potamogeton 
spp.  (pondweeds) ,  plants  restricted  to  oligohaline  and 
low  mesohaline  areas  (Stewart  1962,  Stevenson  and  Confer 
1978) . 

Substrate :  Although  most  benthic  organisms  show  a  cer¬ 
tain  eurytopy  as  to  substrate,  sediment  preferences  do 
exist  (Kinner  et  al.  1974,  Maurer  et  al.  1978).  For 
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example,  sandy  substrates  are  most  numerous  in  the  lower 
Bay,  particularly  near  the  Bay  mouth,  restriction  of  cer¬ 
tain  species  to  this  section  of  the  estuary  is  less  a 
reflection  of  their  stenohalinity  as  it  is  of  their  psam- 
mophilic  nature  (Boesch  1971,  1977) .  Changes  in  fresh 
water  inflow  might  not  only  alter  the  areas  of  certain 
substrate  within  a  particular  salinity  zone,  but  could 
change  sedimentation  rates  and  sediment  types  in  parts 
of  the  Bay  and  its  tributaries  (Hart  and  Fuller  1972, 

Schubel  1972,  Sharaf  el  Din  1977,  Snedaker  et  al.  1977). 

Affected  by  Biological  Interactions :  As  discussed  earlier, 
these interactions  include  predation,  parasitism,  competition, 
and  disease.  Many  estuarine  endemics  and  euryhaline  opportun¬ 
ists  find  the  estuary  a  refuge  from  predation  and  competition 
(Kinne  1966,  Hodgkin  and  Rippingdale  1971,  Boesch  1977).  There 
are  numerous  examples  of  euryhaline  species  having  restricted 
ranges  due  to  incresed  predation  in  higher  salinities.  For 
example,  the  oyster  Crassostrea  virginica  is  predated  in  salin¬ 
ities  above  15  %o  by  the  oyster  drill  Urosalpinx  cinerea,  and 
suffers  heavy  mortalities  in  salinities  above  about  12  %o  due 
to  the  protozoan  parasites  Minchinia  nelsoni  (MSX)  and  Perkin- 
sus  marinus  ("dermo")  (Carriker  1955,  Gunter  1955,  Andrews 
1967,  Sprague  et  al.  1969,  Haven  et  al.  1978).  Although  pre¬ 
dation  seems  to  be  the  most  important  factor,  at  least  for  ben¬ 
thic  forms  (Virnstein  1977,  1979) ,  Evidence  exists  that  compe¬ 
titive  exclusion  may  restrict  ranges  of  some  species,  eg. 

Macoma  balthica  versus  tenta  (Boesch  1971) . 

Represent  Key  Trophic  Links:  Certain  species,  because  of  their 
numbers,  productivity,  or  distribution,  represent  major  links 
in  the  Chesapeake  Bay  food  web.  Results  of  caging  experiments, 
stomach  analyses,  and  laboratory  feeding  studies  have  been  used 
to  identify  major  food  items,  food  selectivity,  ingestion  rates 
and  vulnerability  to  predation  for  candidate  study  species 
(Heinle  1966,  1974;  Burrell  1972;  Perry  and  Uhler  1976;  Homer 


and  Boyton  1978;  Holland  et  al.  1979;  Rawls,  unpublished  M.S.; 
and  others) .  Some  abundant  species  are  numerous  because  they 
have  evolved  means  to  avoid  predation,  and  are  thus  not  key 
trophic  links  (Virnstein  1977,  1979) .  However  they  may  be  impor¬ 
tant  for  other  reasons,  such  as  substrate  modification  or  nutrient 
cycling. 

Perform  Key  Ecosystem  Processes:  These  functions  might  include 
nutrient  recycling,  substrate  modification  or  habitat  produc¬ 
tion.  Benthic  organisms  (particularly  by  the  meio-  and  micro¬ 
components)  ,  as  well  as  zooplankton  and  fish,  excrete  nitrogen¬ 
ous  and  phosphorous  containing  compounds;  these  can  be  utilized 
by  phytoplankton  and  rooted  aquatics  for  primary  production 
(Coull  1973,  Hale  1976,  Durbin  1976,  Taft  and  Taylor  1976, 

Kremer  1977,  McCarthy  et  al.  1977) .  Modification  of  substrate 
can  be  positive  or  negative.  Certain species ,  particularly 
polychaetes,  produce  tubes  which  bind  loose  sediments  and 
stabilize  the  bottom,  allowing  colonization  by  other  organisms 
(Kinner  and  Maurer  1978,  Virnstein  1979).  However,  bioturba- 
tion  by  benthic  infauna,  as  well  as  accumulation  of  fecal  mater¬ 
ial  can  create  a  loose  flocculant  substrate  inhibitory  to  many 
species (Rhoads  and  Young  1970,  Levinton  1977).  Prey-seeking 
behavior  by  fish  and  crabs  can  also  disrupt  the  substrate, 
reducing  numbers  and  diversity  of  species  found  (Orth  1975, 
Virnstein  1977,  1979) . 

Certain  species  so  physically  dominate  their  environment  that 
they  themselves  constitute  the  habitat.  In  Chesapeake  Bay, 
major  examples  are  the  oyster  reef  and  its  associates,  and 
submerged  aquatic  vegetation  beds.  Density  and  diversity  of 
species  in  these  habitats  are  greater  than  in  surrounding  sand 
or  mud  bottoms,  and  productivity  can  be  significantly  higher 
(Marsh  1973,  Orth  1973,  Bahr  1974,  Larsen  1974,  and  Penhale 
1977) .  Many  researchers  consider  these  species  associations  to 
represent  biocoenoses,  with  complex  interactions  between  their 
biotic  components  (Wells  1961,  Marsh  1973).  However,  Larsen 


(1974)  cautions  that,  at  least  in  the  case  of  the  oyster  com¬ 
munity,  few  of  the  associated  species  are  obligates;  rather, 
the  physical  structure  of  the  reef  provides  hard  substrate 
for  a  number  of  epifaunal  species,  as  well  as  shelter  for  a 
variety  of  infauna. 

Commercially  or  Reareationallu  Important  Species:  Organisms 
which  are  harvested  by  man,  or  which  provide  non- consumptive 
recreation  are  the  measure  by  which  the  public  tends  to  gauge 
the  "health  and  productivity"  of  the  estuary.  Most  of  these 
species  (eg.  sport  fish,  crabs,  and  waterfowl)  are  large  and 
conspicuous;  many  feed  fairly  high  on  the  trophic  chain. 

Their  continued  abundance  depends  on  the  integrity  of  the 
trophic  web  supporting  their  populations. 

Threatened  and  Endangered  Species:  A  number  of  threatened  and 
endangered  species  inhabit  the  Chesapeake  Bay  area  (U.S.  Fish 
and  Wildlife  Service  1979) .  Some  of  these  are  birds  and  mammals 
which  are  more  or  less  water-oriented,  and  may  depend  on  the 
estuary  seasonally  or  for  some  aspect  of  their  needs  (food, 
shelter,  etc.).  Examples  are:  Bald  Eagle  (Haliaeetus  leuco- 
cephalus) ,  Delmarva  Fox  Squirrel  (Sciurus  niger  cinereus) . 

Others  are  fish  or  reptiles  which  have  been  known  to  enter  the 
Chesapeake  Bay,  or  which  at  one  time  were  resident.  Examples 
are:  Shortnose  Sturgeon  ( Acipenser  brevirostrum) ,  Maryland 
Darter  (Etheostoma  sellare)  (U.S.  Fish  and  Wildlife  Service 
1979)  . 

Availability  of  Data:  Some  organisms  are  important  based  on 
the  above  criteria,  but  lack  adequate  biomass,  distribution, 
tolerance,  and  trophic  information  to  be  useful  study  species. 
Only  those  that  have  been  well  studied  in  terms  of  distribu¬ 
tion  were  included  in  the  preliminary  list. 


D.  INTERMEDIATE  AND  FINAL  SPECIES  SCREENING  CRITERIA 


1.  Intermediate  Screening 

Reduction  of  the  preliminary  list  of  167  species  required  two 
subsequent  screening  steps.  Selection  of  the  final  study 
species  proceeded  as  follows. 

The  original  list,  with  some  additions  suggested  by  reviewers, 
was  reevaluated  using  eight  criteria  listed  below.  To  facilitate 
this  intermediate  screening,  a  series  of  charts  was  developed 
for  consolidation  of  data  from  numerous  sources.  Information 
gathered  on  these  species  included: 

1.  Salinity  range  and  tolerance,  both  in  the  field  and 
from  laboratory  studies,  for  each  potentially  sensi¬ 
tive  life  stage.  If  the  study  was  from  an  area  other 
than  Chesapeake  Bay,  this  was  noted. 

2.  Temperature  tolerances,  both  from  field  and  laboratory 
information.  Of  particular  importance  were  lethal 
temperatures,  and  temperature  ranges  affecting  periods 
of  reproduction  and  growth. 

3.  Biomass  and  abundance  information,  from  Chesapeake  Bay 
and  other  areas.  Seasonality,  as  it  affected  biomass, 
etc.  was  noted. 

4.  Physiological  rates,  including  respiration,  growth, 
and  production  (of  plants) .  Variation  in  these  rates 
as  correlated  with  salinity  or  temperature  were  noted 
when  available.  Many  of  these  rates  were  taken  from 
studies  conducted  on  candidate  species  in  areas  other 
than  Chesapeake  Bay. 

5.  Preferred  substrates  for  species  (when  applicable) . 

6.  Trophic  relationships,  including  preferred  food  or 
prey,  major  predators,  feeding  rates  and  predation 
rates,  both  from  Chesapeake  Bay  and  from  other  areas. 
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Changes  in  rates  due  to  salinity  of  temperature  were 
noted,  when  available. 

7.  Competitors,  disease  predators,  and  other  limiting 
biotic  factors.  Information  from  areas  other  than 
Chesapeake  Bay  was  taken  when  the  same  species  were 
involved  (i.e.  candidate  species  and  its  competitor 
or  disease) .  Historical  changes  in  distribution  of 
important  diseases  or  predators  due  to  salinity  changes 
were  noted. 

8.  Other  limiting  factors  of  a  physical  nature,  such  as 
light,  depth,  turbidity,  etc.  were  noted  when  the 
information  was  available. 

The  task  of  filling  out  these  charts  for  each  of  the  167  candi¬ 
date  species  provided  a  chance  to  evaluate  the  adequacy  and 
scope  of  information  for  each.  Gaps  in  the  literature  were  so 
extensive  that  some  organisms  were  immediately  eliminated. 
Chiefly  on  the  basis  of  availability  of  data,  as  well  as  appar¬ 
ent  ecosystem  importance,  and  sensitivity  to  a  variety  of  fac¬ 
tors,  the  initial  list  was  reduced  to  a  second  list  of  81  species 
and  associations.  Because  of  the  nature  of  the  available  data, 
associations  rather  then  species  were  used  for  phytoplankton 
and  emergent  aquatic  vegetation.  Year  to  year  variability  in 
dominance  of  individual  phytoplankton  species,  but  relative 
stability  in  the  overall  seasonal  associations  led  to  this  deci¬ 
sion  in  the  case  of  phytoplankton.  Difficulty  in  resolving 
differences  in  Maryland  and  Virginia  wetlands  surveys  necessi¬ 
tated  the  use  of  recurrent  plant  associations  rather  than 
individual  species. 

In  developing  this  intermediate  list,  an  attempt  was  made  to 
include  representative  organisms  from  the  various  Venice 
System  zones,  and  from  the  major  ecological  groups  (i.e., 
estuarine  endemics,  euryhaline  marine,  etc.). 
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2.  Final  Study  Species  Screening 

Final  reduction  of  the  intermediate  species  list  involved  con¬ 
struction  of  a  species  screening  matrix.  The  81  species  and 
associations  were  evaluated  against  eighteen  weighted  factors 
(Table  V-l)  species  being  ranked  on  a  scale  of  0  -  4  in  each 
of  these  criteria.  The  ranking  value  and  the  weighted  value 
were  multiplied  to  give  a  score  for  each  factor,  and  these 
Slimmed  for  a  final  score  for  each  species.  Cutoff  values  were 
assigned  to  these  scores,  and  species  with  scores  above  these 
levels  became  final  study  species  (Table  V-2) .  Cutoff  values 
varied  between  functional  groups  (i.e.  zooplankton , benthos) 
because  all  eighteen  screening  factors  did  not  apply  to  every 
group. 

Considerable  discussion  entered  into  the  assigning  of  weight¬ 
ing  values  to  the  final  screening  criteria,  and  these  values 
generated  a  predictable  amount  of  comment  from  reviewers. 
Selective  judgments  had  to  be  made  in  assigning  weight  to 
the  screening  factors,  and  in  ranking  each  species  against 
them.  However,  it  is  hoped  that  at  least  some  of  the  bias 
inevitable  in  developing  any  list  of  "major  species"  has  been 
avoided. 

Screening  criteria  used  were  essentially  those  discussed  in 
Section  V-C.  However,  these  were  expanded  into  18  major  com¬ 
ponents.  That  is,  the  category  "Performs  Key  Ecosystem  Func¬ 
tions"  was  broken  down  into  "Important  to  Nutrient  Cycling', 
"Affects  Water  Quality",  and  "Modifies  Habitat  for  Other  Spe¬ 
cies"  .  The  greatest  weight  was  given  to  factors  which  could 
be  affected  by  low  freshwater  inflow  (eg.  salinity  sensitivity, 
dependence  on  estuarine  circulation)  or  which  measured  a  species' 
importance  to  the  ecosystem  (trophic  dominance,  biomass,  major 
predator,  etc.).  Also,  by  necessity,  the  availability  of  data 
needed  for  mapping  assessment  of  known  and  potential  habitat, 
and  trophic  information  for  ecosystem  analysis  were  heavily 
weighted.  Values  for  other  criteria  were  assigned  by  comparing 
their  relative  importance  with  the  heavily  weighted  factors 
discussed  above. 


TABLE  V-l 


FINAL  SPECIES  SCREENING  CRITERIA 

Eighteen  screening  factors  used  in  species  selection 
natrix,  with  weighting  values  for  each. 


Sensitive  to  Salinity  4 
Sensitive  to  Circulation  Changes  3 
Sensitive  to  Substrate  2 
Important  to  Nutrient  Cycling  2 
Affects  Water  Quality  1 
Major  Bianass  Contribution  3 
Wide  Distribution  in  Bay  1 
Rare  or  Endangered  Species  1 
Trophic  Importance  4 
Specialized  Food  Requirements  3 
Major  Predator  3 
Major  Competitor  1 
Economic  or  Social  Importance  1 
Opportunistic  Colonizer  1 
Modifies  Habitat  for  Other  Species  2 
Distribution  Data  Available  4 
Trqpluc  Data  Available  3 
Sensitive  Life  Stages  2 
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TABLE  V-2 

Final  Study  Species  List 


Winter/Spring 

Associations 


Summer/Fall 

Associations 


PHYTOPLANKTON  ASSOCIATIONS 


Cyclotella  meneghiniana/Melosira  granulata 
tidal  freshwater  association 

Katodinium  rotundatum/ Skeletonema  costatum 
oligohaline,  low  mesohaline  association 

Asterionella  japonica/ Skeletonema  aostatum 
dominated  mesohaline  association 

Nitsahia  pungens  atlantica/ Skeletonema  aostatum/ Chaeto- 
oevos  spp . 

dominated  polyhaline  association 

Anacystis/ Microcystis 
tidal  freshwater  association 

Gyrmodinium  spp ./Prorocentrum  minimum 

dominated  oligohaline,  low  mesohaline  associations 

Gyrmodi  ni  um/Chae  toce  ros/Ske  le  tonema 

dominated  high  mesohaline  polyhaline  associations 


SUBMERGED  AQUATIC  VEGETATION 


Ceratophyllum  dermersum 
Potamogeton 
Ruppia  maritima 
Zanichellia  palustris 
Zostera  marina 


hornwort 
pondweeds 
widgeon  grass 
horned  pondweed 
eelgrass 


SI 


EMERGENT  AQUATIC  VEGETATION  ASSOCIATIONS 


Tidal  Freshwater  Associations 
Spartina  spp. 

dominant,  brackish  tidal  marsh 

Juncus  roemerianus 
dominant,  brackish  tidal  marsh 


ZOOPLANKTON 


Ctenophora 

Cnidaria 

Rotifera 


Mnemiopsis  leidyi 
Chrysaora  quinquecirrha 
Brachionus  calcyi floras 


ctenophore 
sea  nettle 
rotifer 
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TABLE  V-2 

(Cont . ) 

ZOOPLANKTON  Cont. 

Crustacea 

Aaartia  clausi 

Aoai'tia  tonsa 

Eurytemova  affinia 

Soottolana  canadensis 

cope pod 

H 

H 

H 

Bo8mina  longirostris 

Evadne  tergestina 

Podon  polyphemoides 

cladoceran 

II 

II 

BENTHOS 

Annelida 

Limnodri  lus  hoffmeisteri 

oligochaete  worm 

Heteromastus  filiformis 
Pectinaria  gouldii 
Saolecolepides  virdis 
Streblospio  benedioti 

polychaete  worm 

II 

II 

•• 

Mollusca 

Urosalpinx  cinerea 
Crassostvea  vivginiaa 

Macoma  balthica 

Mercenaria  mercenaria 

Mulinia  lateralis 

Mya  arenaria 

Rangia  auneata 

oyster  drill 

oyster 

Baltic  macoma 

hard  clam 

coot  clam 

soft  clam 

brackish  clam 

Crustacea 

Ampelisca  abdita 

Balanus  improvisus 
Callinectes  sapidus 

Cyathura  polita 

Gammrue  daiberi 

heptocheiru8  plwmlosus 

amphipod 

barnacle 

blue  crab 

isopod 

amphipod 

It 

Palaemonetea  pugio 


grass  shrimp 


TABLE  V-2 
(Cont . ) 


4 

► 


1 


FISH 

Alosa  sapidissima 

American  shad 

Alosa  pseudoharengus 

alewife 

Brevoovtia  tyrannus 

menhaden 

Anohoa  mit chilli 

bay  anchovy 

Leiostomus  xantkurus 

spot 

Menidia  menidia 

Atlantic  silver side 

Micropogon  undulatus 

Atlantic  croaker 

Movone  sax^tilis 

striped  bass 

Morone  ameviaana 

white  perch 

Pevoa  flavescens 

yellow  perch 

WILDLIFE 

(BIRDS) 

Anas  platyrhynchos 

mallard 

Anas  rubripes 

black  duck 

Aythya  valisineria 

canvasback 
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The  rationale  for  each  of  the  weighting  values  follows. 


•  Sensitive  to  salinity.  Since  the  major  antici- 
pacted  effect  of  low  flow  conditions  is  an  alter¬ 
ation  of  salinity  regimes/  this  factor  was  weighted 
"4". 


•  Sensitive  to  circulation  changes.  Changes  in 
circulation  due  to  low-flow  or  altered  salinity 
patterns  can  also  be  anticipated,  and  could  affect 
distribution  of  some  species.  For  this  reason 
this  factor  was  given  a  weight  of  "3". 

•  Sensitive  to  substrate.  Substrate  changes  are  not 
an  anticipated  major  effect  of  low-flow,  although 
the  area  of  specified  substrate  within  a  certain 
salinity  range  will  probably  change.  This  factor 
was  therefore  weighted  "2". 

•  Important  to  nutrient  cycling.  Although  nutrient 
cycling  is  an  important  ecosystem  function,  the 
role  many  species  play  in  it  is  not  well-known. 

To  reduce  bias  in  favor  of  a  few  well-studied 
forms,  this  factor  was  only  weighted  "2". 

•  Affect  water  quality.  A  few  species  can  cause 
deleterious  changes  in  water  quality  (eg.  algae 
blooms) ,  and  these  might  be  enhanced  by  reduc¬ 
tion  in  flushing  rates  due  to  flow  flow.  Since 
these  effects  will  probably  be  local,  this  cri¬ 
terion  was  only  weighted  "1". 

•  Major  biomass  contributor.  Biomass  is  not  only  a 
measure  of  a  species  importance  or  dominance  in 
the  ecosystem,  a  certain  minimum  level  of  abun¬ 
dance  is  necessary  for  a  species  to  be  useful  as 
an  indicator  organism.  Thus,  this  factor  was 
weighted  "3". 

•  Wide  distribution  in  Bay.  This  actually  means 
that  a  species  with  very  restricted,  localized 
distribution  may  not  be  a  useful  indicator.  To 
minimize  bias  for  very  widespread  eurytcJpic 
species,  however,  this  criterion  was  only  weighted 

It  1  •• 


•  Rare  and  endangered  species.  Because  of  the  re¬ 
stricted  ranges  and  usually  minor  ecosystem  impact 
of  these  species,  this  factor  Wc. .  rated  "1".  How¬ 
ever,  it  was  suggested  that  because  of  these  organ¬ 
isms  legal  importance,  the  entire  group  should  be 
handled  as  an  entity  in  the  assessment  (see  Chapter 
V,  Section  E) . 


•  Trophic  importance.  Certain  organisms  are  extremely 
important  to  the  production  and  flow  of  energy  through 
the  estuarine  ecosystem.  Disruption  of  these  species 
could  have  severe  impact  on  other  levels  of  the  trophic 
web.  For  this  reason,  "Trophic  Importance"  was  ranked 
"4". 


•  Specialized  food  requirements.  Species  with  restrict¬ 
ed  food  requirements  at  some  point  of  their  life  cycle 
could  be  more  severely  affected  by  environmental  per¬ 
turbations  than  less  specialized  forms.  Some  species 
may  themselves  be  eury topic  in  regard  to  salinity, 
etc.  but  rely  on  a  more  stenotopic  food  species. 

This  criterion  was  ranked  "3"  for  the  above  reason. 

•  Major  predator.  Predation  has  been  shown  to  be  an 
important  factor  limiting  distribution  of  many  or¬ 
ganisms.  Change  in  distribution  of  a  major  predator 
might  have  significant  effects  on  the  Bay  ecosystem, 
therefore  this  factor  was  ranked  "3". 

•  Major  competitor.  Competition  appears  not  to  be 
as  important  in  mediating  organism  distribution 
as  predation,  so  this  factor  was  ranked  "1". 

•  Economic  or  social  importance.  Although  these  are 
the  factors  through  which  the  public  perceives  the 
Bay's  health,  it  was  felt  that  the  fact  that  an 
organism  was  economically  important  was  not,  a 
priori,  a  measure  of  that  species'  sensitivity  to 
low  flow.  Many  of  these  species  do  have  life  stages 
sensitive  to  salinity  changes,  or  have  predators  or 
diseases  which  could  be  affected  by  low  flow,  but 
these  species  would  receive  high  scores  on  those 
particular  criteria.  For  these  reasons,  this  factor 
was  ranked  as  " 1" . 

•  Opportunistic  colonizer.  Species  which  are  adapted 

to  rapid  colonization  of  disturbed  habitats  may  res¬ 
pond  quickly  to  habitat  alterations  due  to  low-flow. 
However,  too  high  a  ranking  on  this  factor  might  bias 
the  selection  of  species  in  favor  of  estuarine  oppor- 
tunitists,  most  of  which  are  quite  eurytopic.  This 
factor  received  a  "1"  rating. 

•  Modifies  habitat  for  other  species.  Species  which 
provide  habitat  for  other  organisms  (or  conversely, 
which  unfavorably  alter  habitat) ,  can  have  signifi¬ 
cant  ecological  impact.  This  criterion  was 

rated  "2". 

•  Distribution  data  available.  Pragmatically,  it  is 
necessary  to  have  reasonably  accurate  and  complete 
distribution  information  on  a  species  to  either  map 
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it  or  to  assess  changes  in  distribution  due  to  low 
flow.  Thus  this  factor  was  heavily  weighted  as  "4" . 


•  Trophic  data  available.  Complete  and  accurate  infor- 
mation  on  a  species'  ecological  importance  is  needed 
to  assess  what  effects  changes  in  its  distribution 
might  have  on  Bay's  ecosystem.  Thus,  this  factor 
was  also  heavily  weighted,  as  "3". 

•  Sensitive  life  stages.  Many  species  have  a  period 
in  their  life  cycle  which  is  potentially  sensitive 
to  environmental  perturbation;  this  is  typically  a 
larval  or  juvenile  stage.  Although  species  which 
have  such  life  stages  will  also  score  high  on  other 
factors  (such  as  Sensitive  to  Salinity)  it  seemed 
better  to  also  augment  their  score  with  an  additional 
factor.  This  screening  criterion  was  ranked  "2". 

This  matrix-screening  process  produced  a  list  of  57  study 
species.  The  list  was  distributed  to  the  WESTECH  review  team, 
and  presented  at  the  March  20,  1980  conference  for  peer  review. 
Input  from  the  review  process  was  used  to  generate  the  final 
list.  In  particular,  the  fish  species  were  reevaluated  in 
response  to  comments  that  life  stages  should  each  have  been 
screened  independently.  Ranking  each  life  stage  separately 
changed  the  relative  order  to  some  of  the  candidate  species, 
resulting  in  additions  and  deletions  from  the  orginal  list. 

One  species  which  elicited  wide-spread  comment  was  the  American 
shad,  Alosa  sapidissima.  In  light  of  its  severely  depressed 
populations,  its  suitability  as  a  study  species  was  questioned. 
However,  the  apparent  current  stresses  on  this  fish  are  such 
that  additional  pressures  due  to  low  flow  might  prove  critical, 
if  such  effects  can  be  separated  out  and  evaluated.  Thus  the 
species  was  retained  as  a  study  species. 


Several  benthic  species  were  also  reevaluated  on  the  basis  of 
comments,  and  additions  and  deletions  were  made.  In  particular, 
two  species  important  in  the  oligohaline  zone,  the  area  where 
pronounced  effects  of  low  fresh  water  inflows  are  expected, 
were  added.  Gamma rus  dai beri  and  Cy thaura  polita.  Corbicula 
manilens is  was  omitted  due  to  its  limited  distribution,  on  the 
advice  of  reviewers. 
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The  bald  cypress  (Taxodium  distichun.)  was  not  retained  as  a 
study  species  due  to  its  relatively  restricted  distribution. 
Stands  of  cypress  which  may  be  impacted  by  low  flow  will  be 
evaluated  on  a  site-by-site  basis.  Similarly,  rare  or  endan¬ 
gered  species,  or  others  of  special  note,  will  be  assessed  on 
an  individual  basis. 

The  final  study  species  list  contains  fifty-seven  species  and 
associations,  including  seven  phytoplankton  and  three  emergent 
vegetation  associations,  five  submerged  aquatic  vegetation 
species,  ten  zooplankton,  nineteen  benthic  invertebrates,  ten 
fish,  and  three  waterfowl.  Each  of  these  has  been  mapped,  and 
their  distribution  in  regard  to  salinity,  season,  substrate, 
etc.  assessed.  Predators  (such  as  Beroe  ovata)  or  diseases 
(MSX,  "dermo")  which  are  not  themselves  study  species  will  be 
addressed  in  relation  to  the  study  species  they  impact. 

Although  fifty-seven  species  represent  approximately  3  percent 
of  the  total  Chesapeake  Bay  biota,  these  study  species  include 
many  of  the  major  organisms  in  the  estuary.  In  addition,  they 
are  representatives  of  various  salinity  zones  and  estuarine 
habitats,  and  can  serve  as  "models"  for  other  species  with  sim¬ 
ilar  requirements.  Thus  impact  of  low  freshwater  inflows  can 
be  assessed  in  a  specific  manner  for  the  study  species,  and  to 
a  certain  extent  extropolated  for  the  entire  Chesapeake  Bay 
ecosystem. 


E.  RARE,  UNCOMMON  OR  THREATENED  SPECIES 

Selection  of  study  species  for  the  Biota  Assessment  has  been 
carried  out  with  consideration  of  eighteen  factors  which  focus 
on  importance  of  the  organism  to  the  Bay  ecosystem  (see  Section 
V-D  above) .  Due  to  the  particular  requirements  of  the  Low  Flow 
Project,  the  weighting  criteria  emphasized  trophic  and  salinity 
categories  and  deemphasized  such  factors  as  economic  importance 
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and  rarity.  Many  of  the  economically  important  species  were 
included  in  the  final  study  species  list  because  these  species 
(mostly  fish  and  shellfish)  also  rated  high  in  other  categories. 
Most  of  the  Bay's  rare  or  uncommon  species,  however,  due  to 
low  scores  in  other  categor ies,  were  not  included  as  study 
species . 

There  are  several  rather  clear  reasons  that  uncommon  species 
did  not  rate  highly  on  a  system  geared  toward  salinity  and 
trophic  relationships: 

•  insufficient  data  (on  distribution,  feeding, 
salinity  tolerance,  etc.), 

•  not  coupled  tightly  to  estuarine  system, 

•  minor  quantitative  importance  in  food  web,  at  this  time, 

•  Sensitive  stages  not  well  known. 

Many  rare  or  uncommon  species  are  known  only  from  a  few  sightings. 
Aside  from  organisms  endemic  to  a  certain  portion  of  the  estuary 
(i.e.  bald  cypress,  Maryland  darter) ,  other  organism  distributions 
are  known  from  a  spotty,  incomplete  data  base.  Many  are  plants, 
birds  or  mammals  which  are  often  somewhat  independent  of  the 
estuary  and  estuarine  food  webs.  Plants  such  as  bald  cypress 
may  be  found  both  in  estuarine  waters  and  in  non-estuarine  fresh¬ 
water  swamplands  in  the  Bay  region.  The  southern  bald  eagle,  while 
preferring  estuarine  habitat  and  food,  does  also  feed  on  freshwater 
organisms  or  occasional  terrestrial  ones.  Because  of  their 
rarity,  these  species  cannot,  by  biomass  alone,  either  predate 
or  contribute  in  any  major  quantitative  sense  to  the  estuarine 
food  web.  The  presense  of  bald  eagles,  for  example,  may  be  an 
important  indicator  of  ecological  health  and  productivity. 

However,  these  raptors  cannot,  due  to  their  sparse  populations, 
cause  significant  differences  in  the  populations  of  fish  or 
amphibians  on  which  they  prey.  Additionally,  because  of  their 
scarcity,  and  the  justified  reluctance  of  scientists  to  collect 
specimens  of  rare  organisms,  the  life  stages  and  physiological 
tolerance  levels  of  many  of  these  organisms  have  not  been  well 
studied. 
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Lists  of  rare  or  endangered  plant  species  are  published  by  the  U.S 
Fish  and  Wildlife  Service.  However,  only  vascular  plants  are 
considered.  Maryland  lists  9  and  Virginia  43  (5  of  which  are 
overlapping)  rare  or  endangered  plants  (Federal  Register  40:127 
pp27858,  27883-84).  Of  these,  the  majority  are  non-estuarine. 

A  few  such  as  alders  (Alnus  maritima) ,  rushes  ( Juncus  caesarien- 
sis) ,  sedges  (Carex  bil tmoreana ,  Carex  chapmanii)  and  bulrush 
(Scirpus  f laccidifolius)  fall  in  the  category  of  emergent  aquatic 
vegetation;  however,  no  submergent  plants  or  non-vascular  species 
have  been  enumerated.  Submergent  plants,  while  declining  in  many 
areas  of  the  Bay,  are  not  yet  sufficiently  scarce  to  be  listed  as 
endangered. 

Another  significant  plant  species,  although  not  officially  listed 
as  threatened,  is  the  bald  cypress  (Taxodium  distich um) .  The 
northernmost  outliers  of  this  tree  species  occur  in  the  southern 
swamps  of  Maryland  and  Virginia,  including  the  Pocomoke  in 
Maryland  (Shea  1976)  and  the  Chicahominy  and  other  drainages  in 
Virginia  (Smithsonian  Institution  1978) .  It  is  possible  that 
low  flow  salinity  changes  may  affect  the  distribution  and  success 
of  the  small  outlier  populations  of  this  species  in  the  Bay  region 

Only  a  few  animals  are  listed  officially  as  rare  or  endangered 
by  federal  standards  in  Maryland  or  Virginia.  Of  these,  the 
species  associated  with  the  estuarine  ecosystems  of  the  Chesa¬ 
peake  Bay  include  the  Maryland  darter  (Etheostoma  sellare)  and 
the  southern  bald  eagle  (Haliaeetus  leucocephalus  leucocephalus ) . 
In  addition,  shortnosed  and  Atlantic  sturgeon  and  Atlantic 
green,  hawksbills,  and  ridley  turtles  have  been  historically 
important  but  recently  reduced  in  the  Bay.  The  other  species 
are  basically  terrestrial  in  nature,  although  some  (such  as  the 
Delmarva  fox  squirrel)  occupy  nearshore  environments.  In 
addition,  species  such  as  the  osprey  are  present  only  in  limited 
numbers  or  occupy  goegraphical  areas. 

Many  more  species  have  been  listed  as  rare  or  endangered  by  the 
states  themselves.  The  Corps  of  Engineers,  Chesapeake  Bay 
Future  Conditions  report  lists  many  of  tnese  species.  Others 
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are  contained  in  lists  of  endangered  species  published  by  the 
Maryland  and  Virginia  Departments  of  Natural  Resources.  These 
lists  expand  the  national  lists  through  the  addition  of  such 
species  as  shortnosed  sturgeon,  osprey,  eastern  brown  pelican, 
Atlantic  green  and  hawksbill  turtles  and  many  other  aquatic 
species  which  may  utilize  Bay  habitat. 

Uncommon  species,  such  as  those  listed  above,  were  not  included 
as  study  species  for  various  combinations  or  reasons  presented 
in  this  section.  However,  this  in  no  way  diminishes  their 
importance  or  the  vital  ecological  functions  these  organisms 
might  perform  if  their  numbers  were  to  be  increased.  The  Biota 
Assessment  has  used  present  conditions  (1960's  &  70's)  as  an 
environmental  baseline;  however,  rare,  uncommon  or  threatened 
species  usually  have  both  a  genetic  and  ecological  significance  . 
in  their  own  right.  While  not  discussed  further  here,  these  have 
been  treated  by  others  (see  for  example  Clark  1977) .  Effects  of 
low  flow  will  be  evaluated  on  a  species  by  species  basis. 


F.  DISCUSSION  OF  STUDY  SPECIES 

In  order  to  more  clearly  explain  the  characteristics  of  each 
of  the  57  study  species,  individual  species  discussions  have 
been  prepared.  Due  to  their  total  length,  these  discussions 
cannot  easily  be  included  in  the  text.  They  have  therefore 
been  appended  as  Appendix  A;  however,  the  reader  is  cautioned 
that  many  of  the  textual  discussions which  follow  may  depend, 
to  some  extent,  on  knowledge  of  these  species  discussions. 

Below  we  present  an  outline  of  the  main  points  included  in  the 
discussions  of  study  species. 

The  study  species  discussions  first  define  the  most  widely 
accepted  common  name(s)  and  identify  the  type  of  organisms 
(i.e.  calanoid  copepod,  etc.).  General  range  in  the  Bay  and 
any  seasonality  of  distribution  or  behavior  are  then  discussed. 
Sensitivity  to  salinity  or  other  potential  effects  of  low  flow 


conditions  form  the  focus  of  each  species  discussion. 

For  some  species,  relevant  aspects  of  "potential  habitat"  are 
discussed.  The  purpose  of  this  is  to  define  habitat  areas 
which  may  not  have  been  completely  documented  on  the  litera¬ 
ture,  but  which  possess  environmental  conditions  within  which 
the  species  can  survive.  This  is  followed  by  a  brief  discus¬ 
sion  of  the  species  trophic  importance.  The  discussion  ends 
with  a  recapitulation  of  the  particular  selection  factors  which 
were  of  importance  in  the  selection  of  the  study  species  in 
question.  See  Appendix  A  for  individual  species  discussion. 


G.  STUDY  SPECIES  DISTRIBUTION  AND  MAPPING 

Mapping  of  study  species  was  carried  out  on  1:250,000  scale, 
large-size  ("33"  x  54")  base  maps  of  the  Chesapeake  Bay .  These 
maps  are  in  the  Map  Atlas  which  is  on-file  with  the  Baltimore 
Corps  of  Engineers.  The  purpose  of  this  mapping  was  to  provide 
an  atlas  of  species  distribution  under  base  year  (1960  -  1961; 
see  Chapter  IV)  salinity  conditions.  The  upstream  limit  of  all 
mapping  was  defined  to  be  the  head-of-tide  in  each  tributary 
(see  Figure  IV-2  )  and  the  downstream  limit  was  the  Bay  mouth. 
Within  this  area,  basemaps  were  prepared,  using  differing  den¬ 
sity  shading  films  to  indicate  distribution  of  species  and/or 
populations  (depending  on  organism  type).  In  many  cases,  the 
data  permitted  the  elucidation  of  seasonal  shifts  in  distri¬ 
bution  or  abundance,  migration  and  other  pertinent  factors. 
Figures  V-  1  -  V- 2  are  reduced-scale  examples  of  some  of  the 

map  characteristics.  Mapping  technology  and  techniques  were 
to  a  considerable  extent  patterned  after  Lippson  (1973). 

In  some  cases,  it  was  possible  to  differentiate  "known"  versus 
"potential"  habitat  (see  section  II-E) .  Sparcity  of  detailed 
field  data  precluded  doing  this  with  more  species.  Fish,  for 
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example,  have  seldom  been  sampled  on  a  large-scale  geographic 
basis  with  consistant  methodologies.  In  many  cases,  extensions 
of  information  on  known  behavior,  location  and  salinity  tolerance 
were  made  into  areas  in  which  no  data  exists.  This  is  particu¬ 
larly  true  for  most  of  the  eastern  shore  tributaries,  which  have 
been  studied  very  little  in  comparison  with  western  shore  rivers. 

The  mapping  process  consisted  of  several  steps.  Species  dis¬ 
tribution  was  first  determined  from  existing  field  studies, 
aerial  surveys,  map  surveys,  and  other  literature  or  information 
bases.  Points  or  areas  where  organisms  were  present  or  absent 
were  recorded  precisely  on  base  maps  and  color-keyed  to  the 
source  study.  In  most  cases  compilation  of  the  Baywide  map 
for  each  study  species  was  the  result  of  the  juxtaposition  of 
many  studies,  each  on  an  individual  tributary  or  Bay  segment. 
Then,  since  each  individual  study  typically  came  from  differing 
years  or  time  periods,  it  was  sometime  necessary  to  adjust  the 
upstream  or  downstream  limits  of  species  distribution  to  the 
"base  year",  predicated  on  field  data  on  organisms  salin¬ 
ity  tolerances.  The  mapped  locations  and  abundances  thus  repre¬ 
sent  a  best  judgment  base-year  distribution.  Since  there  is  no 
one  year  in  which  all  species  have  sampled,  this  was  found  to 
be  necessary  standardization  procedure. 

The  salinity  information  which  formed  the  basis  for  this 
standardization  was  plotted  on  full  sized  basemaps  in  the  form 
of  lines  of  equal  salinity  values  { isohaline)  derived  from  the 
Chesapeake  Bay  Institute  slack-water  runs  up  the  Bay  main-stem 
as  reported  in  the  Chesapeake  Bay  Salinity  Atlas  (Stroup  and 
Lynn  1963) .  Tributary  salinity  values  were  obtained  from  state 
or  federal  data  bases,  from  biological  studies,  or  interpolated 
or  extrapolated  from  known  values.  Where  fish  or  benthic  sur¬ 
veys  also  reported  salinities  and  station  location,  these 
salinities  were  plotted  and  used  in  defining  isohalines.  Very 
limited  salinity  information  exists  on  eastern  shore  rivers 
during  the  base  year  (  1961  ) .  For  this  reason,  it  has 

sometimes  been  necessary  to  substitute  data  from  other  years. 
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Wherever  organism  distribution  is  mapped  from  non-base  year  data, 
this  has  been  noted  on  that  species  map. 

Depth,  substrate  and  presence  of  other  organisms  has  also  been 
used  to  define  species  distributions.  Species  have  been  classi¬ 
fied  into  known  suitability  for  certain  mappable  substrates 
(sand,  muddy  sand,  etc.)  and  for  certain  depth  categories  (0  -3, 
3-6,  6+  meters).  Other  species  are  known  to  coexist  with  such 
organisms  as  particular  groupings  of  submerged  aquatic  veteta- 
tion.  Base-maps  for  these  parameters  have  been  used  to  define 
or  adjust  potential  habitat  wherever  applicable. 

All  mapping  of  salinity  (as  well  as  depth  and  substrate  to  some 
extent)  represents  a  "snapshot"  taken  of  a  continuous  process. 

The  isohalines  shift  with  the  state  of  the  tide,  seasonally,  and 
from  year  to  year.  The  maps  show  species  boundaries  along  base- 
year  isohalines  (Venice  boundaries)  in  order  to  provide  means  of 
quantifying  the  distribution  of  study  species.  It  should  be 
understood  that  during  other  seasons  or  years  organisms  may  be 
found  outside  of  the  mapped  boundaries,  while  still  occupying 
the  potential  habitat  indicated  by  these  base-year  maps.  The 
maps  should  be  read  as  if  the  organisms  had  been  synoptically 
sampled  at  high-slack  water  during  a  particular  season  of  1960 
and  1961  (Water  Year  1960)  . 

Most  of  the  groups  mapped  were  found  to  have  certain  character¬ 
istics  peculiar  to  the  organism  group  or  the  ways  in  which  it 
has  been  historically  sampled.  Planktonic  organisms  are  often 
mapped  as  associations  since  distribution  patterns  of  many  of 
the  species  making  up  these  groups  are  little  known.  Since 
these  organisms  are  predominantly  affected  by  water  character¬ 
istics,  seasonal  maps  were  prepared  in  all  applicable  instances. 
Substrate,  depth  and  other  organisms  were  not  usually  important 
factors  for  determining  plankton  distribution. 

For  benthic  organisms  in  particular,  substrate,  depth  ,  salin- 


ity  and  other  organisms  are  often  all  important  factors  in 
determining  distribution.  Prior  history  of  a  site  or  area  can 
also  be  relevant.  It  should  be  remembered,  for  example,  that 
the  base-year  period  preceded  the  effects  of  Tropical  Storm 
Agnes,  and  contemporary  benthic  distributions  may  be  displaced 
downstream  from  mapped  distributions.  In  other  cases,  benthic 
surveys  are  even  farther  out-of-date.  The  Maryland  oyster 
grounds  survey  was  last  completed  during  1909  -  1913,  and  signif¬ 
icant  changes  have  occurred  in  the  location  and  extent  of  oyster 
bars  since  that  time.  Virginia  has  completed  an  updated  and 
much  more  accurate  resurvey  (Haven  et  al_.  1977,  1979,  1980)  and 
Maryland  is  attempting  to  resurvey  at  the  present  time. 

Submergent  and  emergent  aquatic  plants  have  both  had  large- 
scale  surveys  completed  in  recent  years.  A  Baywide  survey  of 
submerged  aquatic  vegetation,  using  aerial  photography  (remote 
sensing)  has  been  recently  completed.  These  data  serve  as  the 
basis  for  SAV  mapping.  Each  of  the  bay  states  surveyed  wet¬ 
lands  in  the  1960's  and  early  1970's.  Wetland  inventories  required 
two  years  (Maryland)  to  nearly  a  decade  (Virginia)  to  complete. 
Emergent  vegetation  was  mapped  from  these  county  by  county  wet¬ 
land  surveys. 

Nekton  (mainly  fish)  were  mapped  by  sequential  aggregation  of 
studies.  Maps  were  prepared  of  estuarine  segments  where  studies 
have  been  conducted.  Sampling  stations  were  plotted  and  coded 
by  presence  or  absence  of  the  species,  as  well  as  abundance 
categories  where  data  were  available.  When  density  or  abun¬ 
dance  information  was  collected,  the  sampling  points  were  color 
coded  into  density  ranges.  The  sectional  maps  were  then  pieced 
together  to  make  a  rough-copy  basemap.  The  rough  basemap  of 
known  distributions  was  next  examined  with  respect  to  Venice 
boundaries,  depth,  substrate  type,  etc.,  which  delineate  poten¬ 
tial  habitat  and  necessary  minor  adjustments  made.  The  final 
maps  were  then  created  from  the  rough  basemaps. 
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Aerial  surveys  by  state  management  agencies  (Maryland  Wildlife 
Administration  and  Virginia  Fish  and  Game  Commission)  were  used 
to  define  waterfowl  distribution.  Data  from  3  to  5  year  averages 
were  obtained  wherever  possible,  since  waterfowl  data  showed 
wide  fluctuations  from  year  to  year.  Birds  were  counted  by 
census  tracts  set  up  by  the  agencies.  These  tract  counts  were 
reaggregated  into  Bay  modeling  segments  and  mapped  by  density 
of  birds  per  100  square  kilometers. 

In  Phase  II,  the  Map  Atlas  will  be  supplemented  through  maps  of 
organisms  distributions  based  on  average  inflow  (modal  hydrograph) 
and  drought  scenario  salinity  data  from  Corps  of  Engineers  hydraulic 

model.  The  species  selection  and  mapping  described  in  this  chap¬ 
ter  will  not  only  provide  a  data  base  m  its  own  right,  but  will 
also  provide  a  reference  point  for  assessing  the  relationship 
of  hydraulic  model  data  with  actual  historical  base-year  dis¬ 
tributions  . 


-241- 


CHAPTER  V.  REFERENCES 


Academy  of  Natural  Sciences  of  Philadelphia.  1972.  Potomac 
River  Survey. (Dickerson)  1968  for  the  Potomac  Electric 
Power  Company.  90  pp. 

Andrews,  J.D.  1967.  "Interaction  of  Two  Diseases  of  Oysters 
in  Natural  Waters."  Proc .  Nat .  Shellfish  Assoc.  57:38-49. 

Bahr,  L.  1974.  Aspects  of  the  Structure  and  Function  of  the 
Intertidal  Oyster  Reef  Community  in  Georgia.  Ph ,D.  Disser¬ 
tation.  University  of  Georgia.  ~ 

Boesch,  D.F.  1971.  The  Distribution  and  Structure  of  Benthic 
Communities  in  a  Gradient  Estuary^  Ph .D.  Dissertation . 
College  of  William  and  Mary.  132  pp. 

_ .  1972.  "Species  Diversity  of  Marine  Macrobenthos  in  the 

Virginia  Area."  Ches .  Sci .  13 ( 3) : 206-211 . 

_ .  197  3.  "Classification  and  Community  Structure  of  Macro- 

Benthos  in  the  Hampton  Roads  area,  Virginia."  Mar.  Biol. 21 (3) 
226-244. 

_ .  1977  .  "A  New  Look  at  Zonation  of  Benthos  Along  the 

Estuarine  Gradient."  IN:  Ecology  of  Marine  Benthos.  B.C. 

Coull  (ed.).  University  of  South  Carolina  Press.  Columbia, 
South  Carolina. 

Boon,  J.D.  Ill,  M.E.  Boule,  G.M.  Silberhorn.  1977.  Delineation 
of  Tidal  Wetlands  Boundaries  in  Lower  Chesapeake  Bay  and 
Its  Tributaries.  Spec.  Report  No.  14 0  in  applied  marine 
science  and  ocean  engineering.  VIMS,  Gloucester  Point,  VA. 

127  pp. 

Bosch,  H.F.  and  W.R.  Taylor.  1968.  "Marine  Cladocerans  in  the 
Chesapeake  Bay  Estuary."  Crustaceana  15 (a) : 161-164 . 

_ .  1973.  "Distribution  of  the  Cladoceran  Podon  polyphemoides 

in  the  Chesapeake  Bay."  Mar.  Biol .  19:161-171. 


Bryan,  B.B.  1977.  The  Ecology  of  the  Marine  Cladocera  of 
Lower  Chesapeake  Bay.  Ph.D.  Dissertation .  University of 
VA.  10 1  pp. 


Burrell,  V.G.  1972.  Distribution  and  Abundance  of  Calanoid 
Copepods  in  the  York  River  Estuary,  Virginia,  1968  and 
1 9 70.  Ph .d:  Dissertation .  College  o£  William  and  Mary . 
234 pp. 


-242- 


Burrell,  V.G.  Jr.  and  W.A.  Van  Engel.  1976.  "Predation  by  and 
Distribution  Of  a  Ctenophore  Mnemiopsis  leidyi,  A.  Agassiz, 
in  the  York  River  Estuary."  Estuar.  &  Coastal  Marine  Sci. 
4(3): 235-242 .  - 

Cain,  T.D.  1974.  "Combined  Effects  of  Changes  in  Temperature 
and  Salinity  on  Early  stages  of  Rangia  cuneata."  Va.  J. 

Sci.  25(1) s 30-31 .  —  “ 


Carriker,  M.E.  1955.  "Critical  Review  of  Biology  and  Control 
of  Oyster  Drills,  Urosalpinx  and  Eupleura. "  USF&W  S.S.R:148. 

Carriker,  M.R.  1967.  "Ecology  of  Estuarine  Benthic  Invertebrates: 
a  Perspective."  IN:  Estuarines.  G.H .  Lauff  (ed.).  A.A.A.S. 

Publ .  83,  757  pp. 

Carter,  N.  1980.  Unpublished  data. 

Castagna,  M.  and  P.  Chanley  1973.  "Salinity  Tolerance  of 

Some  Marine  Bivalves  from  Inshore  and  Estuarine  Environments 
in  Virginia  Waters  or  the  Western  Mid-Atlantic  Coast." 

Malacoe  12(1): 47-96. 

Clark,  J.  1977.  Coastal  Ecosystem  Management.  John  Wiley  and 
Sons.  New  York. 

Cory,  R.L.  1967.  "Epifauna  of  the  Patuxent  River  Estuary, 

Maryland,  for  1963-1964.”  Ches.  Sci.  8(2):71-89. 

Cory,  R.L.  and  _  Dressier.  1980.  Unpublished  data. 


Coull,  B.c.  1973.  "Estaurine  Microfauna,  A  Review:  trophic 
relationships  and  microbial  interactions."  IN:  Estuarine 
Microbial  Ecology.  Stevenson,  L.H.  and  R.R.  Colwell  (eds.). 
Univ.  South  Carolina  Press.  Columbia,  SC.  536  pp. 

Dahl,  E.  1956.  "Ecological  Salinity  Boundaries  in  Poikilohaline 
Waters."  Oikos  7:1-23. 

Dahlberg,  M.D.,  P.A.  Dahlborg,  and  P.V.  Morgan.  1973.  An 
Ecological  Survey  of  the  Potomac  River  at  Douglas  Point, 
Maryland!  Pepco  Report,  by~NUS  Corp.,  Cyrus  Wm.  Rice 
Division.  Pittsburg,  PA.  365  pp. 

Davies,  T.T.  1972.  "Effect  of  Environmental  Gradients  in  the 
Rappahannock  River  Estuary  on  the  Molluscan  Fauna."  Geol . 

Soc .  Am.  Memoir  133:263-790. 

Day,  J.H.  1951.  "The  Ecology  of  South  Africa  Estuaries." 

Part  1:  General  Considerations.  Trans.  Roy.  Soc.  So. 

Africa  33:53-9’ 


-243- 


Day,  J.H.  1967.  "The  Biology  of  the  Knysna  Estuary,  South 

Africa."  IN:  Estuaries.  G.H.  Lauff  (ed.).  A.A.A.S.  Publ. 
83,  757  pp. 

Diaz,  R.J.  1977.  The  Effects  of  Pollution  on  Bentic  Communities 
of  the  Tidal  James  River,  Virginia.  Ph.D.  Dissertation 
University  of  Virginia.  158  pp. 

Douglas,  P.A.  and  R.H.  Stroud  (eds.).  1971.  A  Symposium  on 

the  Biological  Significance  of  Estuaries.  Sport  Fishing 
Institute.  Washington,  D.C.  ill  pp.  March  1971. 

Dozier,  H.L.  1947.  "Salinity  as  a  Factor  in  Atlantic  Coast 
Tidewater  Muskrat  Production."  Trans.  North  Am.  Wild. 

Conf .  12:398-420. 

Durbin,  A . G .  1976.  The  Role  of  Fish  Migration  in  Two  Coastal 

Ecosystems:  The  Atlantic  menhaden  (Brevoortia  tyrannus) 

TrT  Narragansett  Bay,  and  the  alewife  (Alosa  pseudoharengus ) 
in  Rhode  Island  ponds.  Ph.D.  Dissertation.  Univ.  of  Rhode 
I s 1 and ,  King s ton ,  RI . 

Federal  Register:  July  1,  1975.  Threatened  or  Endangered  Fauna 
or  Flora:  Review  of  status  of  vascular  plants  and  determina¬ 
tion  of  "Critical  Habitat'1.  U.S.  Fish  and  Wildlife  Service. 

40 : 127 .  Washington,  D.C. 

Goodwin,  F.  Jr.  1970.  "Zooplankton"  IN:  Gross  Physical  and 

Biological  Effects  of  Overboard  Spoil  Disposal  in  Chesapeake 
Bay.  (Interim  and  final  reports) .  Final  Report,  NRI  Spec. 

Rep .  #  3 . 

Grant,  G.C.  1977.  "Seasonal  Distribution  and  Abundance  of  the 
Chaetognatha  in  the  Lower  Chesapeake  Bay."  Estaur.  Coastal 
Mar.  Sci.  5:809-824. 


Grant,  G.C.  and  t.E.  Olney.  1979.  Lower  Bay  Zooplankton 
Monitoring  Program:  an  introduction  to  the  program  and 
results  of  the  initial  survey  of  March  1978.  VIMS  SSR 
#93.  92  pp. 

Gunter,  G.  1955.  "Mortality  of  Oysters  and  Abundance  of  Certain 
Associates  as  Related  to  Salinity."  Ecology  36 ( 4) : 601-605 . 

Hale,  S.S.  1976.  The  Role  of  Benthic  Communities  in  the 

Nitrogen  and  Phosphorus  Cycles  of  an  Estuary.  Univ.  Rhode 
Island  Mr.  Reprint  57~.  Reprinted  from:  Howell,  F.G. 

J.B .  Gentry  and  M.H.  Smith  (ed.).  Mineral  cycling  in 
southeastern  ecosystems.  ERDA  Symposium  series  1975 


Hamilton,  D.H.  and  R.S.  LaPlante.  1972.  Cove  Point  Benthic 
Study .  Annual  Report.  Chesapeake  Biol"!  Lab .  ,  Nat .  Res. 

Inst.  Ref.  No.  72-36.  66  pp. 

Harrison,  W.,  J.J.  Norcross ,  N.A.  Pore,  and  E.M.  Stanley.  1967. 
Circulation  of  Shelf  Waters  off  the  Chesapeake  Bight, 

Surface  and  Bottom  Drift  of  Continential  Shelf  Waters 
Between  Cape  Henlopen,  Delaware,  and  Cape  Hatteras ,  North 
Carolina"!  June  1 9 £> 3 -December  1964.  U .  S .  Dept .  of  Commerce , 
Prof .  Paper  3,  Environmental  Sci.  Ser.  Admin.,  Washington 
D.C.  82  pp. 

Hart,  C.W.  Jr.  and  S.L.  H.  Fuller.  1972.  Environmental 
Degradation  in  the  Patuxent  River  Estuary,  Maryland. 

Conti  No .  1  from  Dept  of  Limnology,  Academy  of  Natural 
Sciences  of  Philadelphia. 

Haven,  D.S.  1957.  "Distribution,  Growth,  and  Availability  of 
Juvenile  Croaker,  Micropogon  undulatus ,  in  Virginia." 

Ecology  38  (1):88-9T1 

_ .  1980.  Personal  communication.  Virginia  Institute  of 

Marine  Science.  Gloucester  Point,  VA. 

Haven,  D.S.,  W.J.  Hargis,  Jr.,  and  P.C.  Kendall.  1978.  The 
Oyster  Industry  of  Virginia:  its  status,  problems  and~ 
promise-!  S".  R.A.M.S.O.E.  No.  168  VIMS. 

Haven,  D.S.,  J.P.  Whitcomb  and  P.C.  Kendall.  1977.  The 

Present  and  Potential  Productivity  of  Baylor  Grounds  in 
Virginia.  Phase  I.  The  Reppahannock  River.  Annual 
contract  report  for  the  period  15  September  1976  to 
15  July  1977.  Contract  No.  3-265-R-l  to  the  National 
Marine  Fisheries  Resources  Commission.  VIMS.  Gloucester 
Point,  VA.  65  pp. 

Haven,  D.S.,  J.P.  Whitcomb,  and  P.C.  Kendall.  1979.  The 
Present  and  Potential  Productivity  of  Baylor  Grounds  in 
Virginia.  Phase  II.  The  Great  Wicomico  and  PiankatanE~~ 
Rivers.  Annual  Contract  Report  for  the  period  15  July 1977 
to  15  May  1978.  Contract  No.  3-265-R-2  to  the  National 
Marine  Fisheries  Service  throug  VIMS.  Gloucester  Point, 

VA.  pp.  108. 

Heinle,  D.R.  1966.  "Production  of  a  Calanoid  Copepod,  Acartia 
tomsa,  in  the  Patuxent  River  Estuary."  Ches .  Sci .  7(2)  59-74. 

.  1969.  Effects  of  Temperature  on  the  Population  Dynamics 

of  Estuarine  Copepods.  Ph. D .  Dissertation.  Univ.  of  Maryland 
132  pp. 

.  1974.  "An  Alternate  Grazing  Hypothesis  for  the  Patuxent 

Estuary."  Ches .  Sci .  15  ( 3 ): 146-150 . 


-245- 


Heinle,  D.R.,  H.S.  Millsaps,  and  C.V.  Millsaps.  1975.  "Effects 
of  Tropical  Storm  Agnes  on  Standing  Crops  and  Age  Structure 
of  Zooplankton  in  Middle  Chesapeake  Bay . "D117-D136 .  IN: 

The  Effects  of  Tropical  Storm  Agnes  on  the  Chesapeake” Bay 
Estuarine  System.  Davis,  Jr.  (ed.).  CRC  Publ.  No.  34. 

Heinle,  D.R.,  S.  Richman,  S.D.  Van  Valkenburg,  and  J.D.  Allan. 
1978.  Seasonal  Interactions  Among  Estuarine  Primary 
Producers  and  Herbivores.  Final  Tech.  Rep.  to  ERDA, 
tfMCEES  Ref.  78-29CBL. 

Herman,  S.S.,  J.A.  Mihursky,  and  A.J.  McErlean.  1968.  "Zooplank¬ 
ton  and  Environmental  Characteristics  of  the  Patuxent  River 
Estuary  1963-1965."  Ches .  Sci .  9(2):67-82. 

Hildebrand,  S.F.  and  W.C.  Schroeder.  1928.  Fishes  of  the 
Chesapeake  Bay.  Bull.  U.S.  Bureau  Fish.  V~.  53.  pt  1. 

(Reprint  1972  by  Smithsonian  Inst.  Washington,  D.C.) 

Hodgkin,  E.P.  and  R.J.  Rippingale.  1971.  "Interspecies 
Competition  in  Estuarine  Copepods."  Limnol.  Oceanogr . 
16:573-576. 

Holland,  A.F.,  N.K.  Mountford,  M.  Keigel,  D.  Cargo,  and  J.A. 
Mihursky,  MMC  and  CBL.1979.  Results  of  Benthic  Studies 
at  Calvert  Cliffs,  Final  Report"!  Md .  DNR  PPSP-MP-28, 

UMCEES  Re  FT  79 -1 1 9 -CBL . 

Homer,  M.  and  W.R.  Boynton.  1978.  Stomach  Analysis  of  Fish 
Collected  in  the  Calvert  Cliffs  Region,  Chesapeake  Bay  - 
197 7.  Final  Report  to:  Md.  DNR  Power  Plant  Siting  Program. 
Annapolis,  MD. 

Hopkins,  S.H.,  J.W.  Anderson,  and  K.  Horvath.  1973.  The  Brack¬ 
ish  Water  clam  Rangia  cuneata  as  Indicator  of  Ecological 
Effects  of  Salinity.'  U.S.  Army  Engineer  Waterways  Experi¬ 
ment  Station,  Vicksburg,  Miss.,  Contract  Rep.  250. 


Jacobs,  F.  1978.  Zooplankton  Distribution,  Biomass,  Biochemi¬ 
cal  Composition  and  Seasonal  Community  Structure  in  Lower 
Chesapeake  Bay".  Ph.D.  Dissertation.  University  of  Virginia 
T05~pp. 

Jeffries,  H.P.  1962.  "Copepod  Indicator  Species  in  Estuaries." 
Ecology  43(4) :730-733. 


246- 


Keefe,  C.W.  1973.  "Standing  Crop  of  Salt  Marshes  Surrounding 
Chincoteague  Bay,  Maryland-Virginia . "  Ches .  Sci.  14(2): 
117-123. 

Khlebovich,  V.V.  1969.  "Aspects  of  Animal  Evolution  Related 
to  Critical  Salinity  and  Internal  State."  Marine  Biol. 
2:338-345. 

Kinne,  0.  1963.  "Salinity,  Osmoregulation,  and  Distribution 

in  Macroscopic  Crustacea."  IN:  Marine  Distributions. 

M.J.  Dunbar  (ed.).  Roy.  Soc.  Canada  Spec.  Publ.  No.  5. 

_ .  1964.  "The  Effects  of  Temperature  and  Salinity  on 

Marine  and  Brackish  Water  Animals.  II.  Salinity  and 
Temperature  -  salinity  combinations."  Oceanogr.  Mar.  Biol. 
Ann.  Rev.  2:281-339. 

_ .  1966.  "Phvsiological  Aspects  of  Animal  Life  in 

Estuaries  with  Special  Refernce  to  Salinity."  Neth.  J. 

Sea.  Res.  3 (2) : 222-244 . 


Kinner ,  P.  and  D.  Maurer.  1978.  "Polychaetous  annelids  of 
the  Delaware  Bay  Region'.'  Fishery  Bull.  76  ( 1)  : 209-224  . 

Kinner,  P.  and  D.  Maurer,  and  W.  Leathern.  1974.  "Benthic 

Invertebrates  in  Delaware  Bay:  Animal -sediment  associations 
of  the  dominant  species."  Int.  Revue,  ges.  Hydrobiol. 

59(5) : 2-18 . 

Kreraer,  P.  1977.  "Respiration  and  Excretion  by  the  Ctenophore 
Mnemiopsis  leidyi . "  Marine  Biology  44:43-50. 

Larson,  P.  1974.  Quantitative  Studies  of  the  Macrofauna 

Associated  with  the  Mesohaline  Oyster  Reefs  of  the  James 
River,  Virginia.  Ph.D.  Dissertation  1974,  College  of  Will¬ 
iam  and  Mary.  Jl5  pp. 

Lear,  D.W.  and  S.K.  Smith.  1976.  "Phytoplankton  of  the 
Potomac  Estuary."  IN:  The  Potomac  Estuary:  Biological 
Resources,  Trends  and  Options.  Mason,  W.T.  and  K.C.  Flynn 
(eds) .  Proc.  Symp.  Interstate  Comm.  Potomac  River  Basin. 
Bethesda,  MD  Tech  Publ.  76-1,  140  pp. 

Levinton,  J.S.  1977.  "Ecology  of  Shallow  Water  Deposit-Feeding 
Communities,  Quisset  Harbor,  Massachusetts."  IN:  Coull  B.C. 
(ed.).  Ecology  of  Marine  Benthos.  Univ.  of  South  Carolina 
Press . 

Lippson,  A . J .  (ed.).  1973.  The  Chesapeake  Bay  in  Maryland; 

An  Atlas  of  Natural  Resources.  Johns  Hopkins  University, 
Baltimore,  MD.  55  pp. 


Lippson,  A . J .  and  R.L.  Moran.  1974.  Manual  for  Indentif ication 
of  Early  Development  Stages  of  Fishes  of  the  Potomac  River 
Estuary.  Power  Plant  Siting  Program,  Md.  Dept  Nat.  Res . 

2S2  pp. 

Lippson,  A . J .  ,  M.S.  Haire,  A.F.  Holland,  F.  Jacobs,  J.  Jensen, 
R.L.M.  Johnson,  T.  Polgar,  and  W.A.  Richkus.  1979. 
Environmental  Atlas  of  the  Potomac  Estuary.  Martin  Marietta, 
Baltimore,  MD.  279  pp . 

McCarthy,  J.J.,  W.R.  Taylor,  T.J.L.  Taft.  1977.  "Nitrogenous 
Nutrient  of  the  Plankton  in  the  Chesapeake  Bay,  1:  Nutrient 
Availability  and  Phytoplankton  Preferences."  Rimnol. 

Oceanogr.  22 ( 6) : 996-1011 . 

McErlean,  A.J.,  C.  Kerby,  R.C.  Swartz,  and  L.C.  Kohlenstien. 

1972.  "Conclusions  and  Recommendations"  Ches.  Sci.  13 
Suppl : 58-516 .  ” 

Mackiernan,  G.B.  1968.  Seasonal  Distribution  of  Dinof lagellates 
in  the  Lower  York  River,  Virginia"  M.A.  Thesis.  College  of 
Williams  and  Mary.  VIMS .  122  pp. 

Marsh,  G. A.  1973.  "The  Zostera  epifaunal  community  in  the  York 
River,  Virginia."  Ches .  Sci .  14:87-97. 

Marshall,  H.G.  1966.  "The  Distribution  of  Phytoplankton  along 
a  140-mile  Transect  in  the  Chesapeake  Bay."  Va.  Journal  of 
Science  17 (2)  : 105-119  . 

_ .  1967.  "Plankton  in  James  River  Estuary,  Virginia.  I: 

Phytoplankton  in  Willoughby  Bay  and  Hampton  Roads." 

Ches.  Sci.  8(2):  90-101. 


Maryland  Dept,  of  Natural  Resources.  1967-1968.  Unpublished 
Data  from  Wetland  Habitat  Inventory.  Water  Resources 
Administration,  Annapolis,  MD. 

Maryland  Dept,  of  Natural  Resoures .  1975-1979.  Maryland 

Midwinter  Waterfowl  Surveys,  Waterfowl  Population'  Estimates . 
MD ,  DNR ,  Wildlife  Adminstration ,  Annapolis,  MD. 

Maurer,  D.,  L.  Watling,  P.  Kinner,  and  C.  Wethe .  1978.  "Benthic 

Invertebrates  Assemblages  of  Delaware  Bay."  Marine  Biol. 
45:65-78. 

Metzgar,  R.G.  1973.  Wetlands  in  Maryland.  Pub.  #157.  Dept,  of 
State  Planning  DNR^  Dept  of  Economic  &  Comm.  Dev. ,  Baltimore, 
MD . 

Morse,  D.c.  1947.  "Some  oberservations  on  seasonal  variations 

in  plankton  population,  Patuxent  River,  Maryland,  1943-1945." 
Chesapeake  Bio.  Lab.  Pub]  .  No.  65,  May  1947  ,  21  pp . 


-248- 


Mountford,  K.  1977.  "Phytoplankton  and  Primary  Productivity" 
IN ;  Academy  of  Natural  Sciences  of  Philadelphia.  Morgantown 
Station  and  the  Potomac  Estuary:  a  316  environmental 
demonstration.  Vol.  II.  PEPCO  Report  by  ANSP,  Sept.  1977. 

Mountford,  N.K.,  A.F.  Holland,  and  J.A.  Mihursky.  1977. 

"Identification  and  Description  of  Macrobenthic  Communities 
in  the  Calvert  Cliffs  Region  of  the  Chesapeake  Bay." 

Ches.  Sci.  18  (4)  : 360-369 . 


Mulford,  R.A.  1972.  "An  Annual  Plankton  Cycle  on  the  Chesapeake 
Bay  in  the  Vicinity  of  Calvert  Cliffs,  Maryland."  June  1969- 
May  1970  .  Proc .  Acad .  Nat.  Sci .  Phila.  124  ( 3)  :  17-40 . 


Munro,  .  1979.  Unpublished  data.  Migratory  bird  habitat 

Research  Laboratory,  Patuxent,  Maryland. 


National  Marine  Fisheries  Service.  1976.  Fishery  Statistics 
of  the  United  States,  1973.  Statistical  Digest  No.  67 . 

U.S.  Government  Printing  Office,  Washington,  D.C. 

_ .  1977  .  Fishery  Statistics  of  the  United  States,  1974. 

Statistical  Digest  No.  68 .  U.S.  Government  Printing  Office 
Washington,  D.C. 

_ .  1978.  Fishery  Statistics  of  the  United  States,  1975. 

Statistical  Digest  No.  (f9^  U.S.  Government  Printing  Office, 
Washington,  D.C. 

NOAA .  1976.  Fisheries  of  the  United  States,  1975.  U.S.  Dept, 

of  Commerce,  NOAA/NMFS,  Washington  D.C.  Current  Fisheries 
Statistics  No.  6900.  100  pp. 

_ .  1977.  Fisheries  of  the  United  States,  1976.  U.S.  Dept. 

of  Commerce,  NOAA/NMFS,  Washington,  D.C.  Current  Fisheries 
Statistics  No.  7200.  96  pp. 

_ .  1978.  Fisheries  of  the  United  States,  1977.  U.S.  Dept. 

of  Commerce,  NOAA/NMFS,  Washington,  D.C.  Current  Fisheries 
Statistics  No.  7500.  112  pp. 

_ .  1979.  Fisheries  of  the  United  States,  1978.  U.S.  Dept. 

of  Commerce,  NOAA/NMFS,  Washington,  D.C.  Current  Fisheries 
Statistics  No.  7800.  120  pp. 

Orth,  R.J.  1973.  "Benthic  Infauna  of  Eelgrass,  Zostera  marina, 
beds."  Ches .  Sci .  14 (4) : 258-269 . 

1975.  "Destruction  of  Eelgrass,  Zostera  marina,  by  the 
Cownose  Ray,  Rhinoptera  Bonasus,  in  the  Chesapeake  Bay. 

Ches.  sci.  16TTr72'05-208  . 


Orth,  R.J.,  K . A.  Mcore ,  and  H.H.  Gordon.  1979.  Distribution 
of  Abundance  of  Submerged  Aquatic  Vegetation  in  the  Lower 
Chesapeake  Bay,  Virginia .  Final  Report  to  EPA,  NO.  6(50/ 
8-79-029/SAV  1,  199  pp . 

Patten,  B.C.,  R.A.  Mulford,  and  J.E.  Warinner.  1963.  "An  Annual 
Phytoplankton  Cycle  in  the  Lower  Chesapeake  Bay'.'  Ches.  Sci. 
4(1) : 1-20. 

Penhale ,  P.A.  1977.  "Macrophyte-epiphyte  Biomass  and  Productivity 
in  an  Eelgrass  ( Zostera  marina)  Community.”  J.  Exp.  Mar. 

Biol .  Ecology  26:211-224. 

Perry ,  M.C.  and  F.M.  Uhler.  1976.  Availability  and  Utilization 
of  Canvas-back  Food  Organisms  in  Chesapeake  Bay.  Pres .  at 
Atlantic  Estuarine  Res.  Soc .  Conf . ,  Rehobeth  Beach , 

Del.  6-8  May  1976. 

Pf itzenmeyer ,  H.T.  1961.  "Benthic  Shoal  Water  Invertebrates 
from  Tidewaters  of  Somerset  County,  Maryland"  Ches.  Sci. 

2(1-2) :89-94. 

_ .  1970 .  "Benthos IfJ:  Gi  jss  Physical  and  Biological  Effects 

of  Overboard  Spoil  Disposal  in  Chesapeake  Bay.  (Interim 
and  Final  Reports).  Final  Report,  NRI  Spec.  Rep.  #3. 

_ .  1973.  Hydrographic  and  Ecological  Effects  of  Enlargement 

of  the  Chesapeake  and  Delaware  Canal.  Appendix  III. "Benthos 
of  Marylancf~Waters  in  ancT  near  C&P  Canal."  NRI  Ref.  #7  3-113. 
CBL,  Univ.  of  Maryland,  CBI,  Johns  Hopkins  Univ. ,  College 
of  Mar.  Studies,  Univ.  of  Delaware.  38  pp. 

_ .  1975.  Benthos.  IN:  A  Biological  Study  of  Baltimore 

Harbor.  CEES,  Univ.  of  Maryland.  Spec.  Rep.  No.  6,  115  pp. 

_ .  1976.  "Some  Effects  of  Salinity  on  the  Benthic  Macro¬ 
invertebrates  of  the  lower  Potomac"  IN:  The  Potomac  Estuary 
Biological  Resources,  Trends  and  Options."  Mason,  W.T. 
and  K.C.  Flynn  (ed.).  Sympos.  Proc. ,  Interstate  Comm. 

Potomac  R.  Basin  Tech.  Publ .  76-2,  140  pp. 

Polgar,  T.T.,  G.M.  Krainak,  and  H.T.  Pf itzenmeyer .  1975.  A 
Methodology  for  Quantifying  the  Responses  of  Benthic 
Communities  to  Environmental  Perturbations.  Morgantown 
Monitoring  Program  Rep.  Ser.  REf.  No.  MT-75-2.  January 
1975.  35  pp.  &  appendices. 

Pritchard,  D.W.  1956.  "The  Dynamic  Structure  of  a  Coastal  Plain 
Estuary. "J.  Mar.  Research  15:33-42. 

_ .  1967.  "Observations  of  Circulation  in  Coastal  Plain 

Estuaries'.'  IN:  Estuaries.  G.H.  Lauff  (ed.).  American  Assoc. 
Adv.  Sci.  Publ.  83,  Washington,  D.C. 


-250- 


Rawls,  C.K.  In  press.  Food  Habits  of  Waterfowl  in  the  upper 
Chesapeake  Bay.  CBL.  unpublished  M.S.  — 

Reinharz,  E. ,  0.  Bricker,  and  A.  O'Connell.  1979.  Chesapeake 
Bay  Earth  Science  Study;  Animal-Sediment  Relationships. 

EPA  Chesapeake  Bay  Program  Data  Report.  8  pp . 

Rhoads,  D.C.  and  D.K.  Young.  1970.  "The  Influence  of  Deposit- 

Feeding  Organisms  on  Sediment  Stability  and  Community  Trophic 
Structure."  J.  Mar .  Res .  28:150-178. 

Ritchie,  D.E.  Jr.  1970.  "Fish."  IN:  Gross  Physical  and  Biolog¬ 
ical  Effects  of  Overboard  Spoil  Disposal  in  Chesapeake  Bay. 
(Interim  and  Final  Report)  Final  Report,  NRI  Spec.  Rep.  #3 

Rupp,  N.H.  1969.  Seasonal  and  Spatial  Distribution  of  Acartia 
tonsa  and  A.  clausi  in  Chesapeake  Bay.  MA  Thesis,  Johns 
Hopkins  Univeristy,  45  pp. 


Sage,  L.E.,  J.M.  Summerfield,  and  M.M.  Olson.  1976.  "Zooplank¬ 
ton  of  the  Potomac  Estuary'.'  IN:  The  Potomac  Estuary: 

Biological  resources,  trends  and  options.  W.T.  Mason  and 
K.C.  Flynn  (eds.).  Symp.  Proc.,  Interstate  Comm,  for  the 
Potomac  River  Basin,  ICPRB  Tech.  Publ.  76-2,  140  pp. 

Sage,  L.E.  and  M.M.  Olson.  1977.  "Zooplankton"  IN:  Morgantown 
Station  and  the  Potomac  Estuary:  a  316  environmental 
demonstration.  Vol.  II.  PEPCO  Report  prepared  by  ANSP. 

Sandifer,  P.A.  1973.  "Distribution  and  Abundance  of  Decapod 
Crustacean  Larvae  in  the  York  River  Estuary  and  Adjacent 
Lower  Chesapeake  Bay,  Virginia,  1968-1969."  Ches.  Sci. 

14(4) :235-257. 

_ .  1975.  "The  Role  of  Pelagic  Larvae  in  Recruitment  to 

Populations  of  Adult  Decapod  Crustaceans  in  the  York 
River  Estuary  and  Adjacent  Lower  Chesapeake  Bay,  Virginia." 
Estuarine  &  Coastal  Mar.  Sci.  3:269-279. 

Schubel,  J.R.  1972.  "  The  Physical  and  Chemical  Conditions  of 

the  Chesapeake  Bay."  J.  Wash .  Acad .  Sci .  62(2):56-87. 

Scott,  R.F.  and  J.G.  Boone.  1973.  Fish  Distribution  in  Various 
Areas  of  Maryland  Tidewater  as  Derived  from  Shore  Zone- 
seining  1956-1972.  Md.  Dept.  Nat.  Res.  Fish.  Admin.,  Annapolis 
MD. 


Seliger,  H.H.,  M.E.  Loftus,  and  D.V.  Subbarao.  1975.  "Dino- 
flagellate  Accumulations  in  Chesapeake  Bay."  IN:  V.R. 
LoCicero  (ed.)  Proc.  First  Intern.  Conf.  on  Toxic  Dino- 
flagellate  Blooms,  Nov.  1974,  Mass.  Sci.  &  Technol.  Found. 
Wakefield,  Mass. 


-251- 


Sharaf  El  Din,  S.H.  1977.  "Effect  of  Aswan  High  Dam  on  the  Nile 
Flood  and  on  the  Estuarine  and  Coastal  Circulation  Pattern 
along  the  Mediterranean  Egyptian  Coast."  Limnol.  Oceanogr. 
22:194-207. 

Shea,  G.B.  1976.  Evaluation  of  Pocomoke  River  Swamp,  Worcester 
CCunty,  Maryland,  for  Designation  as  a  National  Natural 
Landmark  National  Park  Service,  Philadelphia,  PA. 

Smith,  R.F.,  A.H.  Swartz,  and  W.H.  Massman.  1966.  A  Symposium 
on  Estaurine  Fisheries.  Am.  Fisheries  Soc.  Spec^  Publ . 

No.  3^  154  pp. 

Smithsonian  Institution.  1975.  Compendium  of  Natural  Features 
Information,  Vols.  I  and  II.  Center  for  Natural  Areas 
and  Maryland  Dept,  of  State  Planning.  Annapolis,  MD. 

Snedaker ,  S.,  D.  deSylva,  and  D.  Cottrell.  1977.  A  Review  of 
the  Role  of  Freshwater  in  Estuarine  Ecosystems.  tL  Miami 
School  of  Marine  Science,  Final  Report  to  S.W.  Fla. 

Water  Management  District,  May  1977,  126  pp. 

Sprague,  V.  E.A.  Dunnington,  and  E.  Drobeck.  1969.  "Decrease 
in  Incidence  of  Minchinea  Nelsoni  in  Oysters  Accompanying 
Reduction  of  Salinity  in  the  Laboratory."  Proc .  Nat .  Shell¬ 
fish.  Assoc.  59:23-26. 


Stevenson,  J.C.  and  N.M.  Confer.  1978.  Summary  of  Available 
Information  on  Chesapeake  Bay  Submerged  Aquatic  Vegetation 
Bioloqical  Services  Program,  FWS/OBS-78/66 .  MD.  Dept,  of 
Natural  Resources,  U.S.  EPA  Fish  and  Wildlife  Service, 
Annapolis,  MD,  Dept,  of  the  Interior,  Washington,  D.C.  335  pp. 

Stewart,  R.E.  1962.  Waterfowl  Populations  in  the  Upper 
Chesapeake  Region.  U . S.  Fish  and  Wildlife  Serv.  Spec. 

ScTT  Rep.  No.  65,  208  pp. 

Stroup,  E.D,  AND  R.J.  Lynn.  1963.  "Atlas  of  Salinity  and 

Temperature  Distributions  in  Chesapeake  Bay  1962-1961  and 
Seasonal  Averages."  1949-1961.  Graphical  Summary  Report 
2  Ref.  63-1,  CBI ,  Johns  Hopkins  Univ. ,  Baltimore,  MD.  410  pp. 

Swartz,  R.C.  1972.  "Biological  Criteria  of  Environmental 

Change  in  the  Chesapeake  Bay."  Ches .  Sci .  13  Suppl: 17-41. 

Taft,  J.L.  and  W.R.  Taylor.  1976.  "Phosphorus  Distribution 
in  the  Chesapeake  Bay."  Ches .  Sci .  17:67-73. 


Tagatz,  M.E.  1968.  "Biology  of  the  Blue  Crab,  Callinectes 
sapidus,  in  the  St.  John's  River,  Florida."  Fish.  Bull. 

6T:TT-"3  3 . 

Tyler,  M.A.  and  H.H.  Seliger.  1978.  "Annual  Subsurface  Trans¬ 
port  of  a  Red  Tide  Dinof lagellate  to  its  Bloom  Area:  Water 
Circulation  Patterns  and  Organism  Distributions  in  the 
Chesapeake  Bay."  L&O  23 (2) : 227-246 . 


-252- 


U.S.  Army/Corps  of  Engineers  -  Baltimore  District.  1973. 

Existing  Conditions  Report.  Appendix  C,  The  Bay,  Processes 
and  Resources.  Baltimore ,  MD. 

_ .  1977.  Chesapeake  Bay  Future  Conditions  Report.  Vol .  II, 

Biota.  Baltimore,  MD. 

U.S.  Fish  &  Wildlife  Service.  1979.  "List  of  Endangered  and 
Threatened  Wildlife  and  Plants"  Federal  REd.  44(12)  :3636- 
3654. 

Van  Engel,  W.A.  1958.  "The  Blue  Crab  and  Its  Fishery  in  Chesa¬ 
peake  Bay,  Part  I:  Reproduction,  Early  Development,  Growth, 
and  Migration."  Comm.  Fish.  Review  20  (6) :6— 17. 

Van  Valkenburg,  S.D.  and  D.A.  Flemer.  1974.  "The  Distribution 
and  Productivity  of  Nanoplankton  in  a  Temperate  Estuarine 
Area."  Estuarine  &  Coastal  Marine  Sci.  2:311-322. 

Van  Valkenburg,  S.D.,J.K.  Jones  and  D.R.  Heinle.  1978.  "A 
Comparison  by  size,  class  and  volume  of  detritus  versus 
phytoplankton  in  Chesapeake  Bay."  Estuarine  and  Coastal 
Mar.  Sci.  6:569-582. 


Virginia  Commission  of  Game  and  Inland  Fisheries  1975-1979. 

Virginia  midwinter  Waterfowl  Surveys,  Waterfowl  Population 
Estimates .  Va.  Comm,  of  Game  and  Inland  Fisheries,  Tappa- 
hannock, VA. 

Virginia  State  Wetlands  Survey  Series  1973-1979.  VIMS  SRAMSOE 
Reports:  45,  47,  53,  58,  59,  60,  62,  63,  64,  108,  118, 

137,  138,  139,  167,  169,  187,188  (various  authors). 

Virnstein,  R.W.  1977.  "The  Importance  of  Predation  by  Crabs  and 
Fishes  on  Benthic  Enfauna  in  Chesapeake  Bay."  Ecology 
58:1199-1217. 

.  1979.  "Predation  on  Estuarine  Infauna:  Response  Patterns 

~~  of  Component  Species."  Estuaries  2(2):69-86. 

Wass,  M.L.  (ed.)  et  erL.  1972.  A  Checklist  of  the  Biota  of 
Lower  Chesapeake  Bay.  Spec^  Sci.  Report  #65.  VIMS, 
Gloucester  Point,  VA.  290  pp. 

Wells,  H.W.  1961.  "The  Fauna  of  Oyster  Beds,  with  Special 
Reference  to  the  Salinity  Factor."  Ecol .  Mono  31:239-266. 

Whittaker,  R.H.  1970.  Communities  and  Ecosystems.  The  Mac¬ 
Millan  Co.,  New  York"  158  pp. 


-253- 


Willner,  G.R.,  J.A.  Chapman  and  J.R.  Goldsberry.  1975.  A  Stud; 
and  Review  of  Muskrat  Food  Habits  with  Special  Reference 
to  Maryland.  MD  Wildlife  'Admin . ,  Publ.  Wild.  Ecol.  H  2T  p 

Wolff,  W.J.  1973.  "The  Estuary  as  a  Habitat.  An  analysis  of 
data  on  the  soft  bottom  macrofauna  of  the  estuarine  area 
of  the  Rivers  Rhine,  Meuse,  and  Scheldt."  Zool.  Verhandel. 
126:1-243. 

Wood,  L.H .  and  W.J.  Hargis  Jr.  1971.  "Transport  of  Bivalve 
Larvae  in  the  Tidal  Estuary."  IN:D.J.  Crisp  (ed.).  Fourth 
European  Marine  Biology  Symposium,  Cambridge  University 
Press,  Cambridge. 


VI .  THE  CHESAPEAKE  BAY  ECOSYSTEM  MODEL 


This  chapter  describes  the  third  major  conceptual  component  of 
methodology  developed  during  Phase  I,  that  of  ecological  modeling. 
Such  modeling  can,  in  theory,  range  from  the  setting  of  a  con¬ 
ceptual  framework,  to  the  creation  of  a  complex  dynamic  mathe¬ 
matical  computer  model,  as  the  data  permit.  The  extent  to  which 
these  approaches  are  appropriate  depend  on  the  availability 
and  accuracy  of  data  on  ecosystem  interactions.  However,  some 
form  of  modeling,  be  it  comceptual  or  mathematical,  is  necessary 
to  order,  group  and  understand  the  numerous  complex  interactions 
which  comprise  an  estuarine  system  as  large  and  diverse  as 
Chesapeake  Bay.  For  this  reason,  concurrent  with  habitat  class¬ 
ification,  definition  of  salinity  tolerances  and  other  methodol¬ 
ogical  tasks  (see  Chapters  IV  and  V) ,  WESTECH  developed  first  a 
conceptual  f-amework  (model)  of  the  Chesapeake  Bay  ecological 
system  focusing  on  trophic  relationships.  From  this  conceptual 
basis,  drawing  extensively  on  the  scientific  Bay  literature, 
a  mathematical  simulation  model  of  the  Bay  was  developed. 

A.  INTRODUCTION 

One  of  the  objectives  of  the  Chesapeake  Bay  Low  Flow  Study  is 
to  quantify,  to  the  extent  possible,  differences  in  the  produc¬ 
tivity  and  functioning  of  biological  systems  due  to  changes  in 
salinity,  due  to  decreased  freshwater  inflow.  To  this  end, 

WESTECH  has  developed  impact  evaluation  strategies  based  on 
potential  habitat  differences  and  conceptual  modeling  and  ecosystem 
simulation .  In  this  chapter,  the  development  and  structure  of  the 
conceptual  and  mathematical  model  are  discussed.  The  structure 
and  capabilities  of  the  computerized  math  model  are  also  discus¬ 
sed,  as  are  the  limitations  of  ecosystem  simulation.  Several  one- 
year  simulations  of  the  Patuxent  estuary  using  the  Chesapeake 
Bay  Ecosystem  Model  (CBEM)  are  also  included  in  this  chapter. 
Finally,  several  simulations  of  alternative  salinity  scenarios 
in  the  Patuxent  estuary  are  included  to  demonstrate  the  utility 
of  the  CBEM  as  a  tool  in  understanding  salinity-based  ecosystem 
modifications . 
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A  computer  simulation  of  an  ecosystem  attempts  to  reproduce 
some  function  of  that  system,  while  at  the  same  time  basing  the 
simulation  on  a  simplification  of  the  ecosystem  structure  (Hall 
and  Day  1977) .  Structure  includes  the  biotic  as  well  as  the 
abiotic  components,  which  in  real  ecosystems  become  so  numerous 
and  the  relationships  so  complex  that  simplification  is 
essential.  Indeed,  it  is  doubtful  whether  all  the  species  in 
any  natural  ecosystem  have  been  identified  and  counted; 
certainly  all  the  relationships  between  components  are  not  known 
for  any  ecosystem.  However,  it  is  doubtful  that  such  complete 
detail  of  information  is  necessary  for  computer  simulation  to  be 
a  useful  tool  (Patten  1971) .  There  are  many  species  in  biologi¬ 
cal  systems  whose  abundance  is  such  that  they  seem  to  contribute 
little  to  the  trophic  schemes  usually  used  in  modeling,  although 
some  of  those  species  will  be  important  in  regulating  ecosystem 
processes  (Kuenzler  1961,  Connell  1961).  This  is  not  meant  to 
minimize  the  importance  of  such  species  in  long-term  ecosystem 
development  (Darwin  1859)  but  to  point  out  that  daily  and 
seasonal  processes  are  often  dominated  by  a  few  abundant  species. 
In  any  event,  computer  simulation,  as  well  as  much  ecological 
research,  involves  aggregating  species  into  trophic  assemblages 
where  more  or  less  general  storage  and  transfer  rates  are  used 
to  describe  the  assemblages  and  their  inter-relationships.  This 
trophic-dynamic  view  of  ecology  (Lindeman  1942)  can  support 
considerable  theoretical  analysis  (Ulanowicz  and  Kemp  1979) . 

A  computer  simulation  of  an  ecosystem  is  a  tool.  Reproducing 
ecosystem  function  while  being  able  to  manipulate  the  components 
can  provide  insight  in  several  directions.  Ecological  theory 
can  be  explored  with  general  simulation  models  (Kemp  and  Mitsch 
1979,  Geritsen  and  Strickler  1977)  while  simulation  of  a  specific 
ecosystem  aids  in  understanding  processes  in  that  system  as  well 
as  yielding  insights  into  the  more  general  aspects  of  ecosystem 
function  (Kremer  1978,  Anderson  and  Ursin  1977).  Finally, 
computer  simulation  can  be  used  to  predict  changes  in  the 
structural  or  functional  aspects  of  an  ecosystem  based  on  the 
manipulation  of  important  parameters.  This  is  probably  the 
most  difficult  use  of  computer  simulations  since  such  a  model 
should  be  run  using  independent  data  that  was  not  used  to 
initialize  the  model  (Ulanowicz  et  al.  1978). 
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A  number  of  simulation  models  directed  toward  biological  objec¬ 
tives  have  been  or  are  being  developed  for  the  Chesapeake  Bay- 
Simulation  models  are  being  developed  to  investigate  various 
aspects  of  submerged  aquatic  vegetation  in  the  Maryland  Bay 
(Stevenson  et  al.  1979)  and  in  the  Virginia  section  (Wetzel  et 
al .  1979)  as  a  part  of  the  U.S.  Environmental  Protection  Agency's 
Chesapeake  Bay  Program.  Other  EPA  related  model  development  in 
the  Chesapeake  Bay  includes  the  modeling  of  Bay  circulation 
(Shubinski  1979)  and  the  development  of  models  to  identify 
factors  affecting  eutrophication  in  the  Bay  (Ambrose  1979) . 
Ulanowicz  (1976)  has  reviewed  much  of  the  hydrologic  modeling 
literature  important  to  the  Chesapeake  Bay.  O'Connor  et  al. 

(1975)  applied  a  model  of  phytoplankton  dynamics  to  the  Potomac 
estuary.  Their  model  placed  much  emphasis  on  problems  of 
eutrophication  in  the  Potomac.  Numerous  biological  models  have 
been  developed  for  other  aquatic  ecosystems,  such  as  Narragan- 
sett  Bay  (Kremer  and  Nixon  1978) ,  the  North  Sea  (Anderson  and 
Ursin  1977) ,  the  Scramento-San  Joaquin  Delta  (Di  Tora  et  al. 

1971)  and  others. 

The  Chesapeake  Bay  Ecosystem  Model,  in  conjunction  with  mapped 
distributions  and  abundance  information,  is  being  used  in  the 
Biota  Assessment  to  help  understand  the  patterns  and  inter¬ 
relationships  of  species  under  modified  salinity  regimes.  CBEM 
is  not  meant  to  automatically  predict  the  changes  that  occur  when 
salinity  regimes  are  modified.  Rather  it  is  used  heuristically 
in  that  a  number  of  expected,  specific  changes,  such  as  respira¬ 
tion  or  predation  rates,  can  be  programmed  and  CBEM  will  inte¬ 
grate  them  into  the  system.  The  resulting  patterns  are  then 
analyzed  and  interpreted  in  the  light  of  the  information  known 
about  the  real  Chesapeake  Bay. 

To  properly  integrate  this  information,  CBEM  must  simulate  a 
number  of  locations  throughout  the  Bay  where  initial  physical 
and  biological  conditions  differ  substantially.  A  primary 
objective,  however,  is  to  keep  the  number  of  necessary  segments 
to  the  minimum  required  to  estimate  the  effects  of  salinity 
changes  (Ulanowicz  and  Neilson  1974).  Toward  this  end,  the 
Bay  has  been  divided  into  thirteen  segments  (Figure  Vl-1  ) . 
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Except  for  the  Susquehanna  River  segment  (#1),  which  includes 
the  flats,  each  of  the  major  western  shore  rivers  comprise  a 
segment.  The  remaining  segments  are  comprised  of  the  mainstem 
of  the  Bay,  including  the  eastern  shore  rivers.  Each  mainstem 
segment  has  only  one  western  shore  tributary  feeding  into  it. 
Locations  of  pertinent  salinity  changes  within  these  segments 
will  be  selected  for  simulation  in  Phase  II  of  the  Biota 
Assessment.  In  the  sections  below  we  summarize  the  development 
of  CBEM,  beginning  with  conceptual  bay  models. 


B .  CONCEPTUAL  MODELS 

Intensive  compartment  models  of  the  Chesapeake  Bay  were  for¬ 
mulated  through  a  process  of  studying  one  trophic  aggregation 
or  compartment  at  a  time  (i.e.  net-phytoplankton,  macrozooplank¬ 
ton,  demersal  fish,  etc.).  This  process  began  by  defining 
sources  of  food  and  energy,  predators,  life  stages,  seasonal 
migratory  behavior,  requirements  for  nutrients  and  other  inter¬ 
actions  for  each  compartment.  From  these  very  detailed  com¬ 
partment  models  (not  shown)  a  simplified  conceptual  model  of 
the  Chesapeake  Bay  was  constructed  (Figure  VI-2) .  The  symbols 
used  in  the  conceptual  model  are  defined  in  Figure  VI-3. 

The  conceptual  model  shown  illustrates  the  flow  of  energy  from 
the  sun  through  the  plants  and  animals  of  the  ecosystem.  Also 
shown  is  the  movement  of  nutrients  and  non-living  particulate 
matter  through  the  ecosystem.  In  the  model,  radiant  energy  is 
used  by  four  primary  producers  compartments  (net  phytoplankton, 
nannoplankton,  submerged  aquatic  vegetation,  and  emergent  aquatic 
vegetation)  to  produce  plant  tissue.  Two  of  these  plant  compart¬ 
ments,  net  phytoplankton  and  nannoplankton,  produce  material 
which  primarily  enters  a  grazing  food  web.  Zooplankton,  both 
macro- and  micro- (copepods ,  rotifers),  feed  on  these  plants,  as 
do  icthyoplankton,  invertebrate  meroplankton  (oyster,  barnacle 
larvae,  etc.),  forage  fish  and  menhaden.  Benthic  suspension 
feeders  also  graze  the  phytoplankton. 
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The  other  two  primary  producer  compartments,  emergent  aquatic 
vegetation  and  submerged  aquatic  vegetation  (which  includes  a 
smaller  epiphytic  community)  contribute  the  major  portion  of 
their  production  to  the  detrital  food  chain  (although  a  sub¬ 
stantial  amount  is  eaten  by  waterfowl) .  The  detritus  produced 
is  utilized  by  benthic  detritivores  (crabs,  etc.)  and  benthic 
suspension  feeders  (oysters,  clams,  etc.).  Macrozooplankters , 
forage  fish,  and  menhaden  also  utilize  detritus  to  a  certain 
extent . 


The  plankton  species  mentioned  above  (macro-  and  micro-zooplank¬ 
ton,  icthyoplnakton ,  and  invertebrate  meroplankton)  are  fed 
upon  by  ctenophores  and  cniderians  ( comb- jellies  and  sea-nettles), 
fish  such  as  menhaden  and  forage  fish  ( s ilversides ,  etc.),  and 
benthic  suspension  feeders  (such  as  barnacles) .  At  this  point 
in  the  conceptual  model  energy  flows  to  the  predator  layers  of 
the  system.  Pelagic  fish  (bluefish,  etc.),  demersal  fish 
(flounder,  etc.),  and  waterfowl  ( canvasbacks)  are  top  predators, 
and  the  energy  equivalent  of  the  food  they  eat  is  exported  from 
the  system  or  expended  in  feedback  controls  in  the  system. 

The  original  intent  of  the  mathematical  model  development  was 
to  simulate  the  interactions  shown  in  the  entire  conceptual 
model  through  the  use  of  "key"  species  and  interactions.  As 
refinement  of  the  conceptual  model  progressed,  two  factors 
emerged  to  modify  this  approach.  First,  it  became  clear  that 
the  number  of  necessary  key  species  throughout  the  Bay  would 
be  too  large  (on  the  order  of  50  species  or  more)  to  carry  out 
economical  computer  simulations.  Secondly,  community  variation 
in  different  geographic  regions  of  the  Bay  made  selection  of 
Baywide  key  species  difficult.  To  solve  these  two  problems 
and  based  on  suggestions  by  anchor  team  members,  we  elected  to 
focus  on  specific  areas  of  the  Bay,  modeling  mathematically  key 
communities  which  are  subcomponents  of  the  overall  conceptual 
model.  One  such  area  was  selected  for  development  and  cali¬ 
bration  of  the  mathematical  model.  As  previously  stated,  the 
mathematical  computer  model  designed  by  WESTECH  was  developed 
to  help  predict  the  responses  of  an  ecosystem  to  salinity 
changes . 
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The  Patuxent  River  has  been  found  to  have  one  of  the  most 
complete  and  consistant  biological  data  bases  in  the  Chesapeake 
Bay  (Mihursky  and  Boynton  1978) .  Hence,  this  area  provides  a  data 
source  with  which  to  calibrate  and  verify  the  model  with  accur¬ 
ate  historical  data  for  a  number  of  species  and  compartments. 
Therefore,  the  model  was  developed  based  on  selected  organisms  known 

to  inhabit  the  lower  (mesohaline)  reach  (see  Figure  IV-16)  of  the 
Patuxent  estuary.  During  Phase  II  of  the  Study,  the  model  is  expected 
to  be  applied  to  other  zones  and  other  key  geographical  areas 
of  the  Chesapeake  Bay. 

The  final  conceptual  model  used  as  a  basis  for  the  mathemati¬ 
cal  model  is  shown  in  Figure  VI- 4  .  Considerable  simplification 
has  occurred  in  this  model  as  compared  to  the  conceptual  model  in 
Figure  IV-2,  but  the  biological  components  most  critical  in 
estuarine  analysis  have  been  included.  The  phytoplankton, 
zooplankton  (represented  by  Acartia  tonsa  and  A.  clausi)  and 
the  ctenophores  (represented  by  Mnemiopsis  leidyi)  are  fully 
interactive  components  of  the  model.  The  remaining  compartments 
are  present  as  forced  drains,  although  when  feedback  from  the 
interactive  components  to  a  forced  compartment  is  considered 
critical  it  will  be  made  interactive.  In  effect,  the  forced 
compartments  represent  a  built  in  flexibility  of  the  CBEM, 
since  simulations  of  varied  Bay  segments  can  be  accomplished 
relatively  rapidly  by  modifying  the  forced  compartment  and 
the  relevant  abiotic  drivers  (discussed  below) . 


4 
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MATHEMATICAL  AND  COMPUTER  MODELS 


1.  The  Mathematical  Representation  for  the  CBEM 

CBEM  has  been  programmed  in  Fortran  V  for  use  on  the  UNIVAC 
1108  system.  All  model  runs  have  been  conducted  using  the 
University  of  Maryland  1108  Reentrant  Algorithmic  Processor 
(designed  RALPH)  compiler  system.  Some  subroutines  are 
programmed  in  Univac  Assembler  coding.  Standard  Univac 
graphical  packages  are  coupled  to  CBEM  for  display  of  data. 

The  theoretical  basis  for  CBEM  assumes  that  the  main  aspects 
of  the  instantaneous  state  of  the  ecosystem  are  represented 

by  a  vector  X  containing  components  x2 ,  x  , . x^, 

each  of  which  represents  the  energy  or  biomass  of  a  species 
at  one  instant  of  time  (Patten  1971) .  Changes  in  the  biomass  of 
each  species  (growth  or  decline)  are  governed  by  factors  of  food 
(or  energy)  availability,  excretion  and  death  under  the  envir¬ 
onmental  conditions  present.  In  a  mathematical  model,  changes 
in  a  species  x  are  represented  by  the  symbol  x1#  which  generally 
depends  on  the  species  itself  (X3  ),  other  species  (x2,  x3 ,  etc.) 
and  in  some  cases  physical  driving  factors  directly  (i.e.  sunlight, 
nutrients) .  These  relationships  can  be  summarized  in  a  differ¬ 
ential  equation  of  the  form: 

x  =  f  +  ax,  +  bx„  +  .  .  .  .  nx 

where  the  coefficients  a,  b,  etc.  are  each  made  up  of  rate  rela¬ 
tionships  for  assim  lation,  respiration  and  other  factors.  In  this 
context,  the  state  of  the  entire  system  becomes  a  set  of  coupled 
differential  equations  which  represent  the  changes  over  time  of 
x1#  x  .  .  .  .  x  .  In  many  cases,  the  terms  in  these  equations 
are  found  to  depend  on  more  than  one  organism  (i.e.  terms  of  the 
form  ax1  x3etc.).  This  effect  makes  the  differential  equations 
non-linear  and  hence  somewhat  difficult  to  deal  with  computationally 


In  CBEM,  we  use  a  system  of  coupled  quasi-linear  differential 
equations  of  the  form: 


*1  = 

£ol  +  an  *1 

+ 

3 12 

*2  +  a!3  *3  + 

.  ...  a,  x 

In  n 

*2  = 

a2I  X1 

+ 

a22 

X2  + . 

.  .  .  .  x 

2  n  n 

*3  = 

anl  X1 

+ 

an2 

X2  + . 

.  .  .  .  a  x 

nn  n 

each 

coefficient 

"a 

"  is 

determined  by 

combining  rates 

of  assimilation,  respiration,  excretion,  death  and  other  factors 
(import,  export,  etc.)  derivable  from  the  literature.  In  cases 
where  there  is  dependence  on  more  than  one  species  per  term 
(non-linearity) ,  a  quasi-linear  approach  to  the  numerical  inte¬ 
gration  has  been  used.  The  non-linear  term  is  "disguised"  as  a 
constant  within  the  system  of  linear  differential  equiations 
but  is  periodically  recomputed  outside  of  the  equation  set. 

Hence  these  disguised  constants  are  reset  periodically  but  do  not 
introduce  direct  nonlinearity  into  the  equations.  The  constants 
are  manipulated  to  form  part  of  the  a_  or  part  of  the  driving 

functions  f .  .  . 

13 

The  differential  equation  set  is  solved  by  numerical  integration 
techniques  using  a  canned  Runge-Kutta  numerical  integration 
package.  The  advantages  offered  by  this  particular  method  in¬ 
clude  use  of  a  variable  time  step  and  a  flexibility  in  calculating 
the  future  condition  of  the  system  independent  of  its  past  evo¬ 
lution.  This  permits  rapid  solution  of  the  equations  even  during 
periods  of  rapid  growth  or  reduction  of  certain  species. 


2.  The  Model  Drivers  (Light,  heat  and  nutrients) 

Light,  heat  and  nutrients  are  the  basic  abiotic  parameters  in  any 
ecosystem.  Light  is  used  by  photosynthetic  organisms  to  fix 
carbon  into  organic  molecules.  The  plant  material  produced  is 
the  basis  for  the  entire  trophic  structure  of  the  ecosystem. 

Plant  material  is  eaten  and  the  energy  released  is  used  to  build 
animal  tissue  in  both  grazers  and  carnivores.  Fecal  material  and 
dead  animal  and  plant  cells  supply  energy  to  decomposers. 
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The  ambient  temperature  of  the  environment  directly  or  indirectly 
controls  the  rates  of  biological  (and  non-biological)  reactions. 
Photosynthesis,  respiration,  and  ingestion  can  be  controlled  by 
temperature,  as  can  various  reproductive  processes.  Adaptations 
to  long-term  (seasonal) ,  predictable  changes  in  temperature  by 
the  biota  contribute  to  large  scale  patterns  in  ecosystems 
(Slobodkin  and  Sanders  1969) . 

Nutrients  are  components  of  a  great  many  molecules  critical  to 
the  physiology  of  living  organisms.  Nitrogen,  found  in  proteins, 
lipids,  nucleic  acids,  chlorophyll  and  other  molecules;  and 
and  phosporus,  a  component  of  proteins,  lipids,  nucleic  acids 
and  high  energy  molecules  such  as  ATP,  are  often  limiting  to  the 
organisms  in  an  ecosystem. 


In  CBEM ,  nutrient  values  are  input  for  the  particular  river 
basin  from  data  either  as  initial  values  which  are  then  acted 
upon  by  the  biota,  or  as  periodically  corrected  values  based 
on  seasonal  nutrient  values.  Nutrients  are  presently  represented 
as  total  nitrogen  and  phosphorus  or  as  the  limiting  nutrient 
although  it  is  possible  to  use  subcategories  (i.e.  ammonia 
nitrogen)  where  data  exist.  Nutrient  values  change  through 
time  due  to  seasonal  changes  in  input  or  through  interactions 
with  other  biotic  or  abiotic  compartments  in  CBEM. 


Light:  For  the  Patuxent  River  simulation,  long-term  incident 

solar  radiation  data  taken  at  the  Patuxent  Naval  Air  Station 
was  used  (Cinquemani  et  al_.  1978)  .  A  smooth  curve  was  genera¬ 
ted  from  this  data  to  provide  daily  radiation  inputs  into  the 
phytoplankton  compartment.  With  such  a  smooth  average  curve, 
daily  fluctuations  in  incident  radiation  due  to  cloud  cover  are, 
of  course,  lost.  These  fluctuations  appear  to  be  most  impor¬ 
tant  when  phytoplankton  light  acclimation  is  considered  (Kremer 
and  Nixon  1978) . 
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The  effects  of  variable  light  intensity  on  phytoplankton  popu¬ 
lations  in  the  field  is  unclear  however  (Steel  1974),  and  not 
all  apparent  light  acclimation  is  due  to  changes  in  the  photo¬ 
synthetic  response  to  changed  light  intensities  (Yentsch  and 
Lee  1966) .  Because  light  acclimation  by  phytoplankton  is  not 
a  part  of  the  CBEM,  no  daily  variability  has  been  provided 
in  the  sunlight  input. 

The  rate  at  which  light  is  attenuated  as  it  moves  through  the 
water  column  is  affected  both  by  living  and  non-living  material. 

The  living  component  is  assumed  to  be  chlorophyll  and  chlorophyll- 
related  substances,  while  the  non-living  component  includes  both 
the  water  and  substances  in  it  (Riley  1975).  Riley's  (1956) 
equation  relates  phytoplankton  chlorophyll  to  the  total  extinc¬ 
tion  coefficient  of  water: 

k  =  k  +  0 . 054C • 667  +  0.0088C 

w 

where  k  =  total  extinction  coefficient  of  water  (m-1) 

=  extinction  coefficient  of  water  with  no  chlorophyll  (m  *) 
C  =  chlorophyll  a  concentration  (mgm-3) 

k^  values  are  important  in  the  model  when  calculating  the  total 
extinction  coefficients  (k)  produced  as  a  result  of  increases  in 
the  phytoplankton  compartment  (this  self-shading  factor  is 
discussed  below) .  To  determine  k  values  back-calculations  were 

W 

done  using  the  above  equation  and  pubiisned  values  for  chloro¬ 
phyll  a  and  extinction  coefficients  (Flemer  and  Olmon  1971, 

Whaley  et_  al.  1966)  . 

Nutrients:  Nutrient  levels  in  the  sub-estuary  depend  upon 
loading  rates  from  run-off,  inputs  from  the  sediments,  and  uptake 
and  excretion  by  Lhe  organisms  present.  Nutrient  loadings  of 


nitrogren  and  phosphorus  are  comprised  of  point  and  non-point 
sources  (Correll  1976) .  Point  source  loadings  remain  relatively 
constant  over  time  while  non-point  source  loadings  may  be  closely 
tied  to  run-off  (Clark  e_t  al .  1973)  . 

Dilution  and  freshwater  run-off  effects  on  nutrients  provide  a 
means  of  manipulating  the  nutrient  input  to  the  phytoplankton 
compartment  during  base  year  flow  and  low  flow  simulations. 

However,  published  nutrient  budget  data  for  the  Chesapeake  Bay 
and  individual  sub-estuaries  are  not  yet  adequate  to  accomplish 
this  in  a  satisfactory  way  (Heinle  et^  al.  in  press,  Mihursky 
and  Boynton  1978);  however,  best  possible  estimates  will  be 
used.  Since  the  lower  Patuxent  and  other  Bay  sites  are  two¬ 
layered  systems  during  part  of  the  year,  we  plan  to  modify  a  flusing 
equation  based  on  salinities  (Kremer  and  Nixon  1978)  to 
develop  the  capability  to  simulate  nitrogen  and  phosphorus 
mixing  in  and  between  the  model  segments . 

Nitrogen  and  phosphorus  regeneration  from  the  sediment  can  also 
significantly  affect  the  overlying  water,  especially  in  shallow 
estuarine  systems  (Pomeroy  1975) .  Regeneration  is  correlated 
with  temperature,  and  variations  in  nutrient  output  between 
disparate  community  types  may  be  minimal  (Hale  1975) .  Sediments 
may  supply  approximately  30%  of  water  column  nitrogen  demand 
(Boynton  et;  al^.  1977)  .  Benthic  nitrogen  regeneration  in  the  CBEM 

is  based  on  temperature: 

N  =  Nq  +  116.5  (8.66T  -  34.9) 

although  other  facets  of  the  nitrogen  budget,  such  as  nitrogen 
fixation  and  denitrification,  are  still  being  examined. 


Tempera  lure ;  Temperature  is  important  in  the  CBEM,  acting  as 
an  external  interactive  variable  in  almost  all  compartments. 
Temperature  acts  as  a  switch  for  a  number  of  organismic  processes, 
including  the  onset  of  feeding  and  reproduction , and  the  develop¬ 
ment  time  of  eggs  and  larvae  (see  below) .  The  temperature 
inputs  to  the  model  are  from  long-term  data  taken  in  the  Patuxent 
estuary  (Ritchie  and  Genys  1975). 


3.  The  Model  Compartments 

1  inkton :  The  phytoplankton  compartment  is  driven  by 
temperature  to  reproduce  at  a  maximum  rate.  This  maximum  repro¬ 
ductive  rate  is  then  decreased  by  light  and  nutrient  limitations 
(Kremer  and  Nixon  1978,  DeToro  et  al.  1971).  Eppley's  (1972) 
growth  equation: 


log^o  u  =  0.0275  Temperature  -  0.070 


where  u  =  divisions  per  day  (maximum)  at  the  specified  tempera¬ 
ture,  is  converted  to  the  base  of  natural  logarithms  in  order  to 
get  an  instantaneous  growth  coefficient  (Kremer  and  Nixon  1978) : 


rmax 


0 . 59e 


0.0633T 


where  rmax  =  the  maximum  instantaneous  growth  rate,  at  the 
specified  temperature. 


In  the  CBEM  the  equation  has  been  changed  somewhat  to  decrease 
the  growth  rate  below  5°C  in  order  to  account  for  the  warmer 
water  species  present  in  the  Chesapeake  Bay  (Patten  e_t  a_l.  1963)  . 

The  maximum  growth  rate  (rmax)  occurs  when  there  are  no  other 
factors  limiting  the  population.  Realized  growth  rate  occurs 
when  there  are  other  environmental  factors  acting  upon  the 
population  (Krebs  1972,  Odum  1971).  The  two  most  important 
factors  (besides  grazing)  are  light  and  nutrients. 
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The  maximum  growth  rate  (rmax)  is  reduced  by  the  following 
equation  to  account  for  the  effects  of  variable  light  intensity: 

I  e-KZ  _  I 

I  I 

Opt  __e  opt 

where  r  =  realized  growth  rate 

K  =  extinction  coefficient 
Z  =  depth 

F  =  photoperiod  as  a  fraction  of  the  24-hour  day 
I  =  incident  radiation 
I  t  =  optimal  light  intensity 

This  equation  integrates  the  effect  of  non-optimum  light  con¬ 
ditions  both  with  depth  and  with  the  time  of  day.  In  a  column 
of  water  there  will  be  a  point  where  the  amount  of  light  reaching 
the  phytoplankton  is  optimum;  however,  above  or  below  this  point 
photosynthesis  will  decrease  due  to  inhibition  by  too  much  light 
or  decrease  due  to  insufficient  light.  This  depth  of  optimum 
light  will  change  as  the  day  progresses.  I0pt'  the  inten¬ 

sify  at  which  the  photosynthetic  rate  is  maximum,  ranges  from  0.1 
to  0.2  calories/cm2  per  minute  from  spring  to  late  summer  in 
the  model,  to  account  for  the  increasing  dominance  of  dinoflag- 
ellates  over  diatoms  at  this  time  (Lehman  et  a^.  1975,  Eppley 
and  Strickland  1968). 

The  increased  light  attenuation  in  the  water  column  due  to  phyto¬ 
plankton  growth  provides  immediate  self-regulating  feedback  to 
the  phytoplankton  compartment  (Kiefer  and  Austin  1974)  .  In  the 
model,  the  equation  of  Riley  (1956)  previously  discussed  in 
relation  to  calculating  kw  (the  extinction  coefficient  of  water 
without  chlorophyll)  is  used  to  calculate  a  daily  k  (total 
extinction  coefficient  of  water)  based  on  the  phytoplankton 
compartment  size.  This  k  then  is  used  in  the  variable  light 
equation  (see  above) . 


'max 
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Nutrient  limitation  has  been  addressed  in  the  model  by  using  a 
Michael is-Menten  type  equation: 


r 


r 

max 


n 


where  r  =  realized  growth  rates 

n  =  the  concentration  of  the  nutrient  being  considered 

kg  =  the  concentration  of  the  nutrient  at  which  growth 
is  H  the  maximum. 


This  type  of  equation  is  widely  used  in  phytoplankton  nutrient 
relationships  (Eppley  and  Strickland  1968,  Steel  and  Frost  1977, 
DiToro  1980).  The  equation  describes  the  increase  in  growth  as 
the  concentration  of  nutrient  increases  from  much  smaller  to 
much  larger  than  k  . 

Nitrogen  is  the  nutrient  currently  included  in  the  Patuxent 
estuary  simulation,  since  the  P:N  input  values  seem  to  indicate 
the  limiting  roles  of  the  nutrient  (Heinle  et  erl.  in  press, 
Mihursky  and  Boynton  1978) .  CBEM  has  the  capability,  however, 
to  test  for  the  nutrient  most  limiting  to  phytoplankton  growth 
(eitner  N  or  P)  and  use  it  to  reduce  the  maximum  growth  rate. 


The  instantaneous  rate  of  phytoplankton  growth,  rmax  is  reduced 
by  multiplying  r  by  the  nutrient  and  light  equations: 
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which  can  be  integrated  over  the  depth  of  the  water  column. 


CV’pcq  :j>du :  The  copepods  in  the  lower  Patuxent  estuary  are  dom¬ 
inated  by  Acar t ia  tonsa  and  A^.  Clausi .  A.  tonsa  is  most  abundant 
in  the  summer  and  A.  clausi  in  the  late  winter,  but  A.  tonsa  is 
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found  throughout  the  year  (Heinle  1974).  Rates  of  ingestion, 

respiration,  and  to  a  certain  extent  reproduction  are  controlled 

by  temperature  in  the  model.  The  copepod  ingestion  rate  is 

calculated  as  a  function  of  the  rate  of  water  filtered  and  the 

concentration  of  food  in  that  water.  However,  since  copepod 

feeding  rates  are  dependent  upon  the  concentration  of  food  (Lamm 

and  Frost  1976)  a  Michaelis-Menten  type  equation  has  been  added 

(DiToro  1971)  .  The  filtration  rates  for  the  two  species  of  Acartia 

are  based  on  temperature.  The  following  equations  were  calculated 

from  published  data  (Anraku  1964)  and  are  similar  to  the  final 

equations  used  in  CBEM;  although  slight  adjustments  were  found 
to  be  necessary  during  calibration: 

•  Acartia  tonsa 

_  ..  nc_  o.15T 

Fm  =  0.053  e 

•  Acartia  clausi 

Fm  =  0.076  e  °*13T  Ti15°C 

Fm  =  0.534  T> 15  °C 


where  Fm  =  maximum  filtering  rate  (liters  per  calorie  of  copepod 
per  day) 

e  =  base  of  natural  logarithms 
T  =  temperature  °C 


The  above  equations  calculate  the  maximum  filtering  rate  at  that 
temperature.  This  rate  is  then  reduced  by  the  following 
equation: 


F 


=  Fm 


1- 


P 

k  +  P 
s 


where  F  =  realized  filtering  rate  (liters  per  calorie  of  copepod 
per  day 

Fm  =  maximum  filtering  rate  (liters  per  calorie  of  copepod 
per  day 

P  =  phytoplankton  concentration  (calories  per  liter) 

k  =  phytoplankton  concentration  at  which  the  filtration 
rate  is  h  the  maximum  (calories  per  liter) 
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The  amount  of  phytoplankton  filtered  per  calorie  of  copepod  per 
day  is  then  calculated  by  multiplying  the  realized  filtering  rate 
by  the  concentration  of  phytoplankton. 

Basic  respiration  rates  of  Acartia  tonsa  and  A.  clausi  were  also 
calculated  from  Anraku  (1964).  The  respiration  rates  are,  of 
course,  dependent  upon  Q1f)  (physiological)  values  which  can  change 
depending  upon  various  ecosystem  parameters.  Respiration  rate 
equations  may  be  adjusted  in  the  model. 


Copepod  reproduction  is  the  difference  between  assimilation  and 
respiration  (Petrusewicz  1967)  .  This  energy  is  stored  first  as 
eggs,  then  as  juveniles,  until  a  temperature  determined  development 
time  has  elapsed,  at  which  time  they  enter  the  adult  compartment. 
Reproduction  for  Acartia  clausi  is  programmed  to  occur  between 
4  -  20 °C,  wnile  A.  tonsa  has  a  lower  limit  to  reproduction  of 

o  — 

10_C  (Jeffries  1962).  These  values  reflect  in  general  the 
temperature  division  between  these  two  species. 


Hatching  times  for  the  eggs  of  both  species  are  based  on  the 
following  equation  (McLaren  1966,  Mixon  and  Kremer  1978). 


-0 . 1  IT 

e 


In  the  model,  both  Acartia  tonsa  and  A.  clausi  follow  this  equa¬ 
tion,  which  was  calculated  from  a  number  of  species.  At  the  time 
of  hatching  approximately  half  of  the  egg  weight  becomes  nauplii 
weight  (Landry  197  5)  .  The  juveniles  enter  a  development  array 


where  the  projected  development  time,  based  on  temperature,  is 
calculated  for  each  daily  cohort.  Equations  for  development 
were  calculated  from  Heinle  (1966),  Miller  et  a^.  (1977),  and 
Landry  (1975)  . 

The  development  of  juveniles  from  nauplii  to  adults,  encompasses 
a  total  weight  gain  of  approximately  6ug  (Heinle  1966,  Miller 
et  al .  1977)  .  Since  more  than  75%  of  the  weight  gain  by  develop¬ 
ing  juveniles  occurs  in  the  last  25%  of  the  development  time 
(Miller  et  al.  1977)  ,  an  equation  relating  development  stage  to 
weight  has  been  calculated.  Thus  the  juveniles  gain  a  certain 
amount  of  weight  each  day.  This  increase  in  weight,  plus  a 
calculated  respiration  rate  (based  on  temperature)  is  the  juven¬ 
ile  assimilation.  With  an  assimilation  efficiency  of  80% 

(Petipa  1978) ,  the  ingestion  of  the  juveniles  can  be  calculated 
and  subsequently  subtracted  from  the  phytoplankton  compartment. 

Eggs  and  juveniles  are  reduced  by  a  forced  mortality  rate  in  the 
model  so  that  only  a  small  percentage  (less  than  5%)  of  the 
original  juveniles  are  alive  to  become  adults.  Mortality  is 
calculated  daily.  In  addition,  the  juvenile  copepods  are  grazed 
by  the  adult  copepods  in  direct  proportion  to  the  daily  ratio  of 
juveniles  to  phytoplankton.  This  ration  is  usually  insignificant 
except  when  the  phytoplankton  compartment  becomes  small. 

Ctenophores :  The  ctenophore  Mnemiopsis  leidyi  is  one  of  the 

primary  zooplankton  carnivores  in  the  Chesapeake  Bay  (Miller 
1974)  ,  as  well  as  in  other  estuaries  (Kremer  1979) .  In  the 
Chesapeake  Bay,  Mnemiopsis  is  present  in  low  abundance  in  the 
winter  and  early  spring,  but  increases  greatly  in  the  summer 
(Miller  and  Williams  1972) . 

In  the  model,  ingestion  by  Mnemiopsis  is  based  upon  the  amount 
of  water  filtered.  This  is  a  temperature  based  function, 
intially  calculated  from  Kremer  (1979)  : 
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0.51T 


liters  filtered _ 

calorie  of  ctenophor 


per  day  = 


0.08 


e 


Unlike  the  copepods ,  the  ingestion  rate  of  Mnemiopsis  increases 
with  increasing  concentrations  of  prey  (Kremer  1979)  so  that 
the  ingestion  rate  is  calculated  by  multiplying  the  filtering 
rate  by  the  concentration  of  prey. 


Basic  respiration  rates  of  Mnemiopsis  were  calculated  from 

Kremer  (1978) : 


calories  respired 
calorie  of  ctenophore 


per  day  =  0.0068  e 


0.125  T 


Benthic  compartment:  The  model,  at  this  time,  essentially 
simulates  the  water  column  above  an  oyster  (Crassostrea  virginica) 
community.  Three  different  levels  of  oyster  density  are  used, 
depending  on  the  substrate  type  (Dexter  Haven,  personal  communi¬ 
cation)  : 

2 

Rock  bottom  -  500  bushels/acre  -  49  oysters/m 

2 

Sand/Shell  bottom  -  200  bushels/acre  -  20  oysters/m 

2 

Mud/Shell  bottom  -  78  bushels/acre  -  8  oysters/m 

In  the  model,  oysters  remove  phytoplankton  from  the  water  column 
at  a  rate  proportional  to  the  pumping  rate.  Although  the  size 
of  the  phytoplankton  cells  available  will  influence  the  amount 
removed  from  the  water  column  (Haven  and  Morales-Alamo  1967), 
a  generalized  phytoplankton  compartment  renders  such  a  differen¬ 
tiation  unnecessary. 

The  pumping  rate  used  in  the  model  is  the  average  rate  of  6.3 
1/hr  per  grarr;  dry  weight  of  oyster  reported  by  Langefoss  and 
Maurer  (1975).  Complete  clearing  of  the  pumped  water  by  the 
oysters  is  assumed,  but  since  the  oysters  function  in  the  model 
is  to  graze  the  phytoplankton  and  supply  bursts  of  meroplankton, 
no  partitioning  of  the  filtered  material  is  necessary.  Pumping 
by  the  oyster  requires  a  temperature  threshhold  of  8°C  (Galtsoff 
1964),  above  which  maximum  pumping  begins. 


Above  a  temperature  threshhold  of  20  C  oyster  reproduction  begins 
(Galtsoff  1964) .  This  seasonal  input  into  the  plankton  may  have 
significant  effects  on  planktonic  relationships.  Rates  of 
oyster  reproduction,  larval  biomass  and  metabolic  functions  are 
taken  from  the  literature  (Dame  1967,  Galtsoff  1964,  Rodhouse 
1979,  Grant  and  Olney  1979). 

Other  carnivorous  zooplankton:  Fish  larvae  can  be  significant 
predators  on  the  zooplankton  in  Chesapeake  Bay  due  to  their 
abundance  and  residence  time  in  the  Bay.  Larval  stages  of 
hogchoker  (Trinectes  maculatus) ,  spot  (Leiostomus  xanthurus) , 
croaker  (Micropogon  undulatus) ,  and  post  larvae  of  menhaden 
(Brevoortia  tyr annus)  are  common  in  the  Bay  at  certain  times 
of  the  year.  Fish  larvae  are  assigned  a  daily  ingestion  rate 
calculated  as  a  percentage  of  the  dry  weight.  A  temperature 
relation  was  calculated  to  increase  the  daily  ingestion  rate  to 
a  maximum  of  50%  of  the  dry  weight  at  maximum  temperatures 
encountered  (Laurence  1975) .  Larval  fish  become  most  important 
in  the  lower  salinity  reaches  of  the  Bay  tributaries . 

Although  Mnemiopsis  leidyi  is  probably  the  major  carnivorous 
zooplankter  in  the  Bay  (see  above) ,  two  other  zooplankters 
could  be  important  predators.  Beroe  ovata,  a  ctenophore , 
and  Chrysaora  guinguecirrha ,  a  coelenterate ,  feed  upon  Mnemiopsis 
as  well  as  other  organisms  (Cargo  and  Shultz  1967) .  There  is 
evidence  that  B.  ovata  which  feeds  heavily  on  Mnemiopsis  and 
other  ctenophores  (Swanberg  1974)  can  decimate  populations  of 
Mnemiopsis  (Kremer  and  Nixon  1976)  .  For  the  most  part,  infor¬ 
mation  is  lacking  on  these  two  species,  although  several  biomass 
and  feeding  reports  have  enabled  us  to  use  these  organisms 
as  a  possible  drain  on  Mnemiopsis  leidyi  and  the  copepods 
(Kremer  and  Nixon  1976,  Miller  1974,  Swanberg  1974). 


Menhaden  (Bnevoortia  tyrannus):  Adult  (but  sexually  immature) 
and  juvenile  (up  to  one  year  old)  menhaden  are  extremely  numer¬ 
ous  in  the  Chesapeake  Bay,  supporting  a  large  fishery  in  the 
Virginia  portion  (Bell  and  Fitzgibbon  1978,  Frye  1979).  Al¬ 
though  menhaden  have  long  been  assumed  to  feed  to  a  large 
degree  on  phytoplankton,  recent  evidence  points  toward  a  more 
zooplankton-oriented  diet,  especially  where  the  phytoplankton 
is  predominantly  smaller  than  15p  (Durbin  and  Durbin  1975)  . 

Feeding  d ominated  by  zooplankton  prey  may  also  be  true  in  the 
Chesapeake  Bay,  where  up  to  80%  of  the  phytoplankton  have  been 
reported  to  be  less  than  10y  in  diameter  (Van  Valkenburg  et_  al . 

1978)  .  Although  biomass  values  for  menhaden  in  the  Chesapeake 
Bay  are  few,  Durbin  (1976)  has  given  biomass  values  and  meta¬ 
bolic  relationship  for  menhaden  in  Narragansett  Bay.  Carter 
(personal  communication)  has  taken  biomass  values  in  the  Choptank 
River,  Maryland.  In  the  Chesapeake  Bay  model,  menhaden  can  act 
as  a  forced  drain  on  either  the  zooplankton  or  the  phytoplankton, 
or  both,  and  the  impact  of  such  a  drain  investigated.  The  form 
and  intensity  of  possible  menhaden  drains  were  tested  as  shown  in 
section  6,  but  were  generally  found  to  be  realistic  in  the  5-10%  range. 

4 .  Computer  Model  -  Program  Structure 

CBEM  is  programmed  in  Fortran  V  for  use  on  the  UNIVAC  1108 
system  (see  above) .  The  main  program  and  approximately  21 
sub-routines  comprise  CBEM.  The  main  program  coordinates  the 
linear  flow  of  the  simulation  by  calling  the  appropriate 
subroutines.  Subroutines  fall  into  four  general  categories: 

1.  Input  oriented  subroutines 

2.  Internal  subroutines 

3.  Output  oriented  subroutines 

4.  Salinity  oriented  subroutines 

Functions  not  calculated  in  the  subroutines  (i.e.  phytoplankton 
biomass)  are  calculated  in  the  main  program. 
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The  following  lists  names  and  describes  the  primary  subroutines 
in  CBEM: 

1)  Input  oriented  subroutines 

•  CSOLAR  -  calculates  the  daily  input  of  solar 

radiation  (calories  per  square  meter  per  day) . 

•  CTEMP  -  calculates  the  daily  water  temperature. 

•  DATA  -  contains  the  input  data  needed  to  run  the  pro¬ 

gram.  Initial  compartment  values,  Michaelis  - 
Menten  constants,  the  length  of  time  the  program 
is  to  simulate,  and  other  values  are  contained 
here.  Coded  instruction  for  various  output 
formats,  such  as  graphs,  are  also  included  in 
this  subroutine. 

•  FOTO  -  the  maximum  phytoplankton  growth  rate,  rmax»is 

calculated,  as  are  the  limiting  values 
based  on  nutrients  and  light.  r,  the  realized 
growth  rate  (before  predation)  is  calculated 
in  the  product  of  these  three  values. 

•  MONTH ,  SEASON  -  these  two  subroutines  locate  the 

simulation  output  in  useful  time  frames. 

2)  Internal  Subroutines 

•  CALCMX  -  calculates  the  basic  matrix  from  which  the 

differential  equations  are  calculated. 

•  DEGJUV  -  the  initial  amounts  of  time  left  as  copepod 

eggs  and  juveniles  are  calculated. 

•  DELAY  -  calculates  the  daily  hatching  and  growth  rates 

of  copepod  eggs  and  juveniles,  and  mortality 
rates,  including  predation  by  adults. 

•  DEPROD  -  calculates  the  total  weight  gain  of  copepod 

juveniles  and  reduces  the  phytoplankton 
by  an  amount  proportional  to  this  gain. 

•  FUNC  -  creates  the  differential  equations  from  the 

program  matrix. 

•  NUTRNT  -  the  nutrient  compartments  are  increased  and 

decreased  in  this  subroutine. 

3)  Output  oriented  subroutines 

•  GRAPH  -  graphs  the  requested  statistics. 

•  PRTDIF  -  the  differential  equations  are  printed. 

•  PRSTAT  -  the  requested  block  of  statistics  are  printed 

4)  Salinity  oriented  subroutines 

•  SALINE  -  contains  the  basic  changes  applied  when 

salinity  scenarios  are  shifted.  Species  presence, 
predation,  respiration  and  other  factors  change 
with  salinity. 
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5.  Model  Ouput 

The  basic  simulation  period  of  CBEM  currently  being  used  is  one 
year.  Compartment  values  can  be  printed  daily,  or  any  number  of  days 
may  be  omitted  from  the  printout.  Generally  one  day  and  five  day 
increment  printouts  proved  most  useful  in  model  development. 

Three  printout  "packages"  are  basic  to  the  model  and  provide 
most  of  the  information  required  to  understand  the  simulation. 

The  first  printout  contains  the  standing  crop  (cal/m^)  of  the  inter¬ 
active  compartments  and  their  ingestion  and  respira  i<  a  rates.  The 
month,  day,  and  water  temperature  are  also  lifted  in  ^his  print¬ 
out,  as  are  the  values  for  nutrients  and  the  ilized  phytoplank¬ 
ton  growth  rate.  An  example  of  this  printout  format  is  shown 
in  Table  VI-  1. 

The  second  printout  package  is  useful  in  analyzing  the  results 
of  a  simulation.  The  two  limiting  values  ( see  p.277)  calculated  from 
nutrients  and  light  are  listed  in  this  printout,  as  are  the 
values  for  the  extinction  coefficient  of  the  water  column,  the 
depth  of  the  euphotic  zone,  daily  insolation,  chlorophyll  a,  and 
the  ratio  of  daily  insolation  to  optimum  insolation.  The  month 
and  day  are  also  listed  in  this  printout,  as  is  the  maximum  rate 
of  phytoplankton  growth  (calculated  from  temperature)  (Kremer 
and  Nixon  1978)  . 

The  third  printout  format  lists  the  productivities  (cal/day)  of  the 

phytoplankton,  copepods,  and  ctenophores,  along  with  the  basic 
standing  crop  information  for  those  compartments.  An  example 
of  this  printout  format  is  shown  in  Table  VI- 2  . 

Host  of  the  above  data  can  be  output  in  graphical  format,  at 
any  interval  desired.  Graphs,  of  course,  show  patterns  that 
might  be  difficult  to  observe  in  standard  numerical  printouts. 
(Examples  of  graphical  model  output  are  shown  in  subsection  6.) 
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Initial  conditions  have  been  defined  by  utilizing  the  best 
possible  values  of  winter  data  (on  or  about  Dec.  15)  in  the 
literature  on  standing  crop  for  each  organism.  In  some  cases 
this  data  is  corrected  seasonally,  as  when  a  predator  is  known 
to  enter  the  estuary  in  the  spring  (and  may  not  be  present  in 
winter;  i.e.  Chrysaeora) .  The  data  come  from  a  wide  variety 
of  sources  which  have  been  discussed  in  Chapter  III,  many  of 
which  are  cited  in  this  chapter.  Complete  documentation 
on  any  particular  species  is  available  from  WESTECH;  however, 
the  prime  sources  for  the  species  modeling  include: 

•  phytoplankton  -  Heinle  1974,  Mackiernan  (unpubl.  data) 

•  Acartia  spp.  -  Heinle  1966,  Stross  and  Stottlemeyer  1965 

•  Ctenophore  -  Mihursky  and  Boynton  1978,  Kremer  and  Nixon 

1968 

•  predators  -  Cargo  and  Shultz  1967,  Durbin  and  Durbin 

1975,  Homer  and  Boynton  1978,  Miller  1974, 
Swanberg  1974 

Physical  data  was  derived  from  actual  solar  measurements  in 
the  Bay  area  (Cinquemani  et  al.  1978)  as  are  temperature 
measurements  (Ritchie  and  Genys  197  5)  .  Nutrient  data  was  derived 
from  USGS  files  and  Mihursky  and  Boynton  1978. 

In  many  cases,  species  feeding,  respiration  or  other  parameters 
were  not  known  exactly  for  the  Patuxent  estuary.  In  other  cases, 
data  existed,  but  across  some  range  of  values.  The  model  was 
calibrated  by  selecting  and  running  CBEM  with  various  biologically 
realistic  values  within  such  ranges  of  values,  or  through  making 
biologically  defensible  judgements  for  minor  alterations  in 
values.  Such  judgements  were  based  on  other  data  whenever 
possible.  For  instance,  phytoplankton  growth  temperature  depen¬ 
dence  was  determined  in  some  cases  by  using  analysis  of  Chesa¬ 
peake  Bay  chlorophyll  data  to  correct  productivity  calculations 
originally  made  on  Narraganset  Bay  (Kremer  and  Nixon  1978) . 

The  most  biologically  realistic  runs  were  then  selected  to 
provide  the  calibrated  values. 
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6.  Agreement  with  Observed  Data  (Model  Validation) 


Although  some  of  the  basic  data  used  to  program  the  responses 
of  organisms  in  CDEM  comes  from  laboratory  studies  or  studies 
on  ecosystems  other  than  the  Chesapeake  Bay,  the  physical 
parameters  which  drive  and  control  the  system  must  come  from 
the  system  being  modeled.  The  physical  factors  currently  being 
utilized  in  the  model  are  those  from  the  lower  Patuxent  estuary. 

Whenever  possible,  biolgical  data,  such  as  seasonal  standing 
crop,  were  taken  from  studies  done  in  this  river  reach.  When 
values  were  not  available,  such  as  for  the  seasonal  abundance 
of  the  ctenophore  Mnemiops is  leidyi ,  research  conducted  at  other 
locations  was  used  for  data. 

Similar  constraints  were  present  for  the  historical  time  per¬ 
iod  in  which  the  data  was  taken.  Ideally,  the  data  used  should 
be  as  current  as  possible  to  refract  the  present  state  of  the 
area  being  modeled.  Even  better,  the  modeling  effort  should 
be  integrated  into  a  research  program  in  which  the  two  approaches 
can  complement  each  other  and  indicate  the  best  directions 
for  both  the  research  and  modeling.  CBEM ,  however,  has  been 
developed  to  utilize  the  present  knowledge  of  the  Bay  and  is 
not  associated  with  a  research  program.  We  have  used  data  from 
a  number  of  historical  periods  (in  general  different  from  that 
data  used  for  calibration)  to  compare  with  the  output  from  CBEM 
for  purposes  of  model  validation.  Details  of  the  results  of 
several  simulations  are  discussed  below  and  compared  with  reported  data. 

Phytoplankton  values,  as  simulated  by  the  model,  in  general 
showed  a  good  fit  to  published  values  (Figure  VI- 5  ).  Winter 
and  early  spring  values  follow  the  observed  values  well,  as 
does  the  spring- summer  bloom.  Summer  values  fall  a  little 
lower  than  the  data  might  indicate  bit  are  still  above  the 
minimum  observed  .alues.  In  the  early  fall  the  model  shows  a 
large  phytoplankton  bloom  that  is  not  indicated  in  the  observed 
data.  This  bloom  is  inherent  in  the  simulation.  Higher  numbers 
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of  copepod  grazers  shortens  the  duration  of  the  bloom,  and  de¬ 
lays  its  onset,  but  the  intensity  of  the  phytoplankton  growth 
is  not  depressed  (see  below) .  This  aspect  of  the  simulation  is 
currently  under  study  and  may  be  due  to  inaccurate  nutrient  data, 
missing  drains  or  other  factors.  After  this  bloom,  however,  phyto¬ 
plankton  levels  agree  very  well  with  observed  values. 

CBEM  simulations  of  copepod  biomass  compare  favorably  with  the 
rather  limited  data  available  (Figure  VI- 6  ) .  Acartia  tonsa 
values  are  low  in  the  winter  and  decrease  steadily  through  May, 
after  which  they  increase  rapidly.  The  summer  peak  of  A.  tonsa 
agrees  well  with  the  observed  data.  After  abundant  copepods 
in  July  and  August  there  is  a  decrease  beginning  in  September. 

The  model  simulation  shows  an  elevated  abundance  of  copepods 
midway  through  the  fall  decline,  while  the  observed  data  shows 
the  copepods  to  continue  declining  to  low  abundance.  These  higher 
levels  of  copepods  in  the  fall  reflect  the  interactions  occurring 
between  copepods  and  the  fall  bloom  of  phytoplankton  previously 
discussed.  A  large  amount  of  energy  is  flowing  into  the  cope¬ 
pods  at  this  time  from  the  phytoplankton. 

Data  reported  from  a  less  saline  reach  of  the  Patuxent  estuary 
(Heinle  1966)  indicates  a  somewhat  different  pattern  than  that 
discussed  above  (Figure  VI-  7  ) .  Here  A.  tonsa  (adults)  show  a 
peak  abundance  in  March  and  then  a  slow,  somewhat  irregular 
decline  to  the  end  of  May.  There  seems  to  be  no  winter  decline, 
although  overall  abundance  levels  are  lower.  This  is  not  un¬ 
expected,  however,  since  iu  the  less  saline  areas  of  the  river 
A.  tonsa  is  found  with,  and  presumably  competes  with,  Eurytemora 
af finis  (Heinle  and  Flemcr  1975)  . 

Acartia  claus i  values  from  the  CBEM  simulations  (Figure  VI-7  ; 
were  slightly  higher  than  the  United  observed  data  (Pans 
unpublished  data  -  nor.  shown)  .  A.  claus  i  is  essentially  gone 
from  the  Patuxent  estuary  by  June,  and  this  is  reflected  in 
the  simulation  run.  Herman  et  al.  (1968)  found  that  A.  claus l 
became  increasingly  dominant  in  comparison  to  A  tonsa  ;  r<  •: 

March  to  mid-May  (Figure  VI- 8  )  and  this  is  also  1 1 
the  model  simulation. 
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Figure  21- 6  SIMULATED  AND  OBSERVED  COPEPOD  ABUNDANCE 
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As  previously  stated,  information  on  the  abundance  of  Mnemiopsis 
leidyi  in  the  Patuxent  estuary  is  sparse.  Information  is 
available  on  the  combined  abundance  of  Mnemiopsis  and  Chrysaora 
quinquecirrha  (Mihursky,  McErlean  and  Herman  1967) .  Ziegenfuss 
and  Cronin  (1958,  reported  in  Mihursky  and  Boynton  1978)  reported 
an  increase  in  ctenophore  abundance  beginning  in  early  July  end 
peaking  in  late  summer.  Peak  densities  for  Solomons  and  Broomes 
Island  were  9  and  32  individuals  /m3.  This  would  correspond  to  vol 
umes  (#)  of  ctenophores  of  261  and  928/m3  respectively  (Bishop 
1967)  .  These  values  are  high  compared  to  the  peak  of  33  /m3 
found  in  the  York  River  (Burrell  and  Van  Engel  1976)  and  a  peak 
of  about  50  /m3  reported  from  Narragansett  Bay  (Kremer  1976) . 

The  maximum  values  found  in  Broomes  Island  samples  by  Mihursky 
McErlean  and  Herman  (1967)  for  ctenophores  and  Chrysaora  com¬ 
bined  was  43  /m^.  The  CBEM  simulation  peaked  at  100  cal/m^ (21/m3) 
(Figure  VI-  9 ) .  The  model  simulation  show  a  rapid  increase 

of  this  species.  This  explosive  population  growth  is  character- 

* 

istic  of  Mnemiopsis  (Kremer  1976) . 

The  CBEM  simulation  of  the  lower  Patuxent  estuary  is  in  general 
agreement  with  the  observed  data.  The  early  fall  simulations 
show  a  burst  of  energy  going  through  the  system  which  elevates 
the  phytoplankton  and  copepod  levels  beyond  expected  levels. 
However,  levels  quickly  fall  to  the  range  of  observed  data. 


D.  ALTERNATIVE  SALINITY  SCENARIOS 

Several  very  preliminary  simulations  were  done  under  a  generalized 
scenario  of  "increased  salanity"  corresponding  to  changes  of  2-5%.  These 

included  the  addition  of  B.  ovata,  a  ctenophore  predator  of  Mnemiopsis. 
Although  Beroe  has  been  reported  from  the  Patuxent  estuary  (Hermal  §t  al . 
1968)  ,  it  does  not  occur  there  in  appreciable  numbers  except 
during  low  flow  periods.  It  is  normally  restricted  to  areas 
of  higher  salinity  (Burrell  1972).  Three  levels  of  Beroe 
predation  were  simulated;  10$,  20%  amd  50%.  As  a  preliminary 
investigation  into  the  effects  of  physiological  changes 
associated  with  modified  salinity  regimes  we  altered  the  stand¬ 
ard  physiological  response  of  A.  clausi  (Figure  VI-10) . 
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Figure  31-10  SIMULATED  RESPONSE  OF  Mnemiopsis  fidyj  WHEN  THE 
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(Jeffries  1962) .  At  higher  salinities  A.  clausi  can  apparently 
compete  with  A.  tonsa  more  successfully  as  the  temperature 
increases.  Since  competition  between  developing  stages  of  the 
two  species  might  be  most  intensive  (Conover  1956)  adults  and 
juveniles  were  considered  in  these  preliminary  simulations. 

Adding  Beroe  as  a  predator  of  Mnemiopsis  had  several  interesting 
results  (Figures VI- 11-13) .  Beroe  was  programmed  to  enter  the 
Patuxent  estuary  in  August  (day  246)  and  remain  until  the  end 
of  the  year  (day  365) .  A  daily  predation  rate  of  10%  of  Mnem¬ 
iopsis  had  a  pronounced  effect  reducing  the  peak  biomass  by 
approximately  80%  (Figure  VI-11) .  A  daily  predation  rate  of  20% 
by  Beroe  does  not  reduce  the  peak  standing  crop  any  further 
(because  this  now  occurs  before  Beroe  becomes  a  factor) ,  but 
the  decline  in  Mnemiopsis  abundance  is  very  rapid  (Figure  VI -.12)  . 
At  a  predation  rate  of  50%  the  decline  of  Mnemiopsis  is  even 
more  rapid  (Figure  VI-13) .  Such  drastic  reductions  in  Mnemiopsis 
abundance  are  known  to  occur  in  areas  of  overlap  of  these  two 
species  (Burrell  and  Van  Engel  1976) . 

The  effects  of  the  reduced  standing  crop  of  Mnemiopsis  on  the 
copepod  Acartia  tonsa  are  immediately  evident  (Figure  VI-14, 
VI-15) .  A.  tonsa  abundance  is  increased  and  the  late  summer 
copepod  decline  is  delayed  until  November,  a  difference  of 
about  three  months.  There  is  little  difference  in  A.  tonsa 
biomass  at  the  three  levels  of  Beroe  predation  (the  response  to 
50%  predation  is  not  shown; .  The  other  copepod  in  the  model, 
Acartia  clausi  showed  no  change  during  this  simulation  because 
Beroe  entered  the  run  after  A.  clausi  had  declined. 

The  reduction  of  Mnemiopsis  by  Beroe  had  two  effects  on  the 
phytoplankton.  The  fall  bloom  of  phytoplankton  was  delayed 
several  weeks,  and  the  peak  abundance  of  the  bloom  was  reduced 
(Figure  VI-  16,  vi-  17)  .  An  interesting  result  is  that  the  fall 
phytoplankton  bloom  is  slightly  larger  in  the  run  simulating 
a  20%  predation  rate  than  in  the  10%  predation  run.  This  appar¬ 
ently  resulted  when  the  increased  grazing  released  more  nitrogen. 


STANDARD  SIMULATION 


SIMULATED  ABUNDANCE  OF  Mmmioptit  frity 
PREDATED  AT  10%  OAILY. 


STANDARD  SIMULATION 


SIMULATED  ABUNDANCE  OF  Mnemiopsis  Hidyi 
PREDATED  AT  20%  DAILY. 


p 


Figure  2143  SIMULATED  ABUNDANCE.  OF  Mntm/op$i»  fid/I 
PREDATED  AT  SO%  DAILY. 


N  -2»6- 


A.  fl- 


A.  fonsa 
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Figure  31-17  SIMULATED  ABUNDANCE  OF  PHYTOPLANKTON 

WHEN  Mn«miopsi$  IS  PREDATED  AT  A  DAILY  RATE  OF  20% 


The  phytoplankton  response  to  the  50%  Mneroiopsis  loss  (not 
shown)  was  very  similar  to  the  20%  simulation. 

The  second  preliminary  series  of  simulations  modified  the 
physiological  response  of  Acartia  clausi  to  increasing  temper¬ 
atures.  A.  clausi  naupleii  were  programmed  to  survive  until 
the  temperature  exceeded  24°C  (the  preceeding  runs  terminated 
naupleii  survival  at  20°C  (Jeffries  1962) .  The  assimilation 
to  respiration  ratio  of  adults  was  also  altered.  In  previous 
runs  the  A/R  ratio  was  altered  to  stress  A.  clausi  in  an 
exponential  fashion  beginning  at  approximately  13.5°C.  In 
these  salinity  modification  runs  the  stress  did  not  begin  until 
20°C. 

In  this  simulation,  A.  clausi  did  not  begin  to  decline  in  late 
April  as  previously  noted,  but  increased  its  biomass  several 
orders  of  magnitude  to  a  peak  biomass  in  June,  after  which  the 
temperature- induced  stress  caused  a  decline  in  abundance 
(Figure  VI-18) .  Although  there  appears  to  be  only  a  limited 
effect  on  A.  tonsa  (Figure  VI-18) ,  it  was  enough  to  cause  an 
increase  of  about  25%  in  the  spring  phytoplankton  bloom  (Figure 
VI-19) .  A.  tonsa  was  apparently  kept  low  enough  by  competition 
with  A.  clausi  during  the  growth  stage  of  the  bloom  to  allow 
it  to  grow  very  large.  The  small  decrease  in  the  fall  phyto¬ 
plankton  bloom  is  probably  due  to  a  small  number  of  A.  clausi 
now  present  because  of  the  modified  physiological  response. 

The  final  effect  of  this  salinity  modification  run  was  to 
increase  the  abundance  of  Mnemiopsis  earlier  in  the  year, 
probably  as  a  result  of  feeding  on  A.  clausi. 

E.  MODEL  USES  AND  LIMITATIONS 

The  preliminary  altered-salinity  simulations  discussed  above 
are  examples  of  the  way  in  which  CBEM  can  be  used  as  a  tool  to 
aid  in  understanding  ecosystem  processes.  Specific  responses 
to  salinity  are  being  programmed  into  CBEM  whenever  the  data 
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are  available.  These  responses  typically  include  addition  or 
deletion  of  predators  or  new  food  sources,  change  in  respiration, 
or  change  in  susceptibility  to  other  stresses  (i.e.  temperature). 
These  changes  affect  many  of  the  study  species  depending  on 
specific  conditions.  When  the  data  are  not  available,  the  response 
must  be  estimated  using  the  best  available  information.  Either 
way,  the  possibilities  for  combinations  of  factors  are  large  and 
a  computer  is  needed  to  sort  out  the  various  results.  These 
results  are  not  exact  predictions  of  behavior  of  the  real  eco¬ 
system.  They  are  the  end  product  of  feeding  a  large  quantity 
of  information  into  a  computer  and  asking  for  results  within  a 
highly  structured  framework.  These  results  can  indicate  factors 
and  patterns  which  might  be  important  in  the  real  system.  They 
can  also  indicate  areas  where  the  available  data  do  not  seem  to 
explain  a  simulated  occurrence,  such  as  the  fall  phytoplankton 
bloom  in  the  CBEM  simulation  (see  above) . 

CBEM  can,  within  data  limitations,  be  utilized  to  study  salinity- 
related  low  flow  effects  through  analysis  of  alternative  salinity 
scenarios.  Given  a  change  of  2-3%*  salinity,  for  instance, 
the  first  task  of  the  CBEM  user  would  be  to  identify  major  biolog¬ 
ical  effects  (i.e.  salinity  limitation  of  predators  on  food 
supply,  respiration  changes,  etc.).  These  effects  would  then  be 
entered  into  the  CBEM  input  matrix.  Obviously,  if  the  literature 
defines  no  known  changes,  the  scenarios  will  be  identical;  simul¬ 
ation  can  predict  no  more  than  the  state-of-the-art  in  input 
data.  If,  however,  differences  do  exist,  CBEM  will  give  the 
interactive  results  of  a  trophic  web  of  organisms  acting  under 
the  altered  salinity  conditions  for  one  or  more  years.  This 
can  then  be  checked  and  calibrated  against  known  biological  be¬ 
havior  insofar  as  data  exists.  The  power  of  the  simulation  is 
to  predict  wholeistic  interactions  of  the  individual  components. 

In  a  less  direct  sense,  CBEM  can  be  coupled  to  other  changes 
(i.e.  nutrients)  known  to  change  with  flow,  but  only  insofar  as 
base  data  or  mathematical  calculations  allow. 

It  is  clear  that  without  low  flow  biological  data  to  check  against, 
CBEM  results  will  only  be  verifiable  insofar  as  they  duplicate 


"reasonable"  biological  occurrences.  They  can,  however,  be  used 
as  a  tool  to  analyze  relationships  and  predict  biological  res¬ 
ponse  scenarios . 

Results  of  the  CBEM  can  be  used  then,  to  predict  potential 
differences  in  the  ecosystem  given  two  or  more  sets  of  input 
variables.  Finally,  although  there  are  severe  dif f iciencies  in 
the  information  available  concerning  the  Bay,  there  is  still  a 
great  deal  of  data  available  from  the  research  that  has  been 
done  (e.g.  Chesapeake  Science  vols.  1-18).  The  goal  of  the  CBEM 
is  to  utilize  the  available  data  in  such  a  way  as  to  provide 
insight  into  the  effects  of  reduced  freshwater  inflow  into  the 
Chesapeake  Bay. 
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VII.  DATA  GAPS  AND  IMPACT  ASSESSMENT  METHODOLOGY 


This  chapter  is  intended  to  first  summarize  (Section  VII-A) 
the  major  areas  of  missing  information  or  knowledge  in  the  bio¬ 
physical  data  base  on  the  Chesapeake  Bay.  In  Section  VII-B  the 
summary  will  be  used  to  develop  a  rationale  for  impact  assess¬ 
ment  methodology  which  will  be  used  in  Phase  II  of  the  Biota 
Assessment.  Since  the  information  summarized  in  Section  A 
has  been  previously  referenced  in  Chapters  III  -  VI,  few 
reference  sources  are  re-cited  here.  The  reader  is  referred 
to  the  appropriate  sections  of  the  earlier  chapters. 

A.  DATA  GAPS 
1.  Physical/Chemical 

While  the  underlying  functional  principles  of  Bay  circulation 
have  been  described,  little  is  known  about  certain  current 
mechanisms.  This  lack  is  particularly  severe  in  terms  of  quan¬ 
titative  descriptions  of  currents  at  depth  and  their  precise 
modes  of  interactions  with  surface  currents.  The  velocity  of 
the  longitudinal  Bay  current  which  transports  many  organisms 
up  and  down -bay  can,  at  present,  only  be  inferred  from  the  bio¬ 
logical  data.  Its  interactions  with  currents  from  river  mouth??, 
eddy  currents  or  other  area  specific  phenomena  are  not  at  all 
well-known. 

In  the  tributaries,  the  effects  of  vertical  fronts  have  just 
begun  to  be  investigated,  although  these  may  strongly  affect 
biota  distributions  in  local  situations.  Eddy  patterns 
caused  by  such  fronts,  or  by  bottom  features  or  current  inter¬ 
action  have  only  been  identified  on  a  few  of  the  better-studied 
western  shore  tributaries. 

In  particular,  the  relationship  between  low  or  high  flows  and 
stratification  or  eddy  phenomena  has  not  yet  been  well  docu¬ 
mented  in  a  way  that  permits  generalization  to  other  than  speci¬ 
fic  cases.  Also,  the  effects  of  freshets  on  the  physical  charac 
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teristics  and  their  relationship  to  low  flow  has  not  been  syn- 
theized,  and  although  such  information  may  exist  in  long-term 
monitoring  records,  its  compilation  would  be  a  major  task  in 
itself. 

Correlations  between  flows  and  nutrients  have  been  carried  out 
on  a  few  river  systems;  however,  the  data  base  does  not  yet 
seem  strong  enough  to  formulate  general  rules  which  hold  with 
any  accuracy  in  most  Bay  situations.  The  relationship  between 
other  apsects  of  water  quality  and  flow  is  generally  similar 
to  that  of  nutrients.  It  seems,  however,  that  some  rough-cut 
correlations  between  certain  nutrient  components  and  flows  can 
be  made,  at  least  for  the  major  rivers.  It  may  then  be  possible 
to  extrapolate  these  relations  to  other,  smaller  river  systems. 

The  relationship  between  the  "turbidity  maximum"  and  flow  is 
not  yet  totally  clear,  although  the  location  of  this  zone  has 
been  charted  under  differing  flow  conditions.  One  of  the  com¬ 
plications  here  is  the  lack  of  proportional  reductions  in  flows 
of  different  river  systems  under  historical  drought  conditions. 
Some  river  reduce  flows  more  than  others,  and  similar  patterns 
are  rarely  repeated  twice  in  the  historical  data. 


2.  Biological 

Many  of  the  biological  data  gaps  involve  the  structure  and  func¬ 
tion  of  the  biotic  community  itself, while  others  involve  the 
effects  of  the  physical  or  chemical  system  on  the  biota.  It 
should  be  noted  that  the  level  of  detail  of  information  of  all 
types  is  strongly  species  dependent,  often  favoring  species  of 
commercial  or  recreational  importance.  While  this  is  a  logical 
result.,  it  does  pose  problems  when  one  tries  to  assemble  a  bal¬ 
anced,  detailed  picture  of  the  ecological  system. 

Distribution  of  organisms  tends  to  be  well-known  in  localized 
areas,  usually  only  for  certain  species  or  species  groups. 
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Studies  tend  to  cluster  around  institutions,  areas  of  water 
pollution  problems,  potential  power  plant  sites  and  other  cen¬ 
ters  of  activity  or  potential  environmental  problems.  For  these 
reasons,  some  of  the  smaller  rivers,  especially  those  with 
rural  watersheds,  have  not  been  well-studied.  The  Rappahannock 
and  most  of  the  eastern  shore  rivers  are  examples  of  this. 

In  some  cases,  even  physical  data  such  as  salinity  are  only 
available  for  these  tributaries  from  scattered  time  periods. 

Dietary  information  for  many  of  the  predators  are  not  well- 
known.  Although  stomach  content  analysis  has  been  performed 
for  certain  species  of  waterfowl  or  fish  in  certain  locations, 
variability  of  diet  with  location,  season,  food  availability 
or  other  factors  has  seldom  been  studied ,  probably  due  to  the 
complex  logistics  that  such  studies  would  entail.  In  develop¬ 
ing  the  CBEM  (Chapter  VI)  for  example,  copepod  feeding  rates 
were  found  to  be  fairly  well-defined,  while  feeding  rates  of 
ctenophores  had  to  be  developed  mostly  from  information  out¬ 
side  of  Chesapeake  Bay.  Many  predators  are  apparently  able  to 
switch  food  sources  fairly  easily,  while  others  cannot  suc¬ 
cessfully  do  so.  The  question  of  the  abilities  of  predators 
to  switch  could  potentially  be  one  of  the  more  important  aspects 
of  the  actual  impact  of  low  flows,  or  any  other  stress,  on  the 
ecosystem  as  a  whole.  Predators  with  a  wide  range  of  potential 
food  sources  would  be  much  more  resistant  to  adverse  effects 
than  those  incapable  of  switching.  There  is  evidence  that 
some  species  of  waterfowl  were  able  to  change  diets  during 
the  drastic  changes  in  SAV  populations  that  have  occurred  in 
the  last  two  decades,  although  the  importance  of  this  "switching" 
ability  has  not  yet  been  elucidated. 

Biomass  data  is  often  unknown  or  largely  unpublished  for  such 
major  groups  as  fish.  Sampling  techniques  tend  to  differ  and 
there  are  often  no  obvious  ways  of  standardizing  information  ob¬ 
tained  from  different  catch  techniques.  Often  sufficient  repli¬ 
cation  of  samples  is  not  conducted  to  assure  statistical  relia- 
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bility  of  any  standardization  attempt.  In  other  cases,  even  the 
methods  of  measuring  biomass  (e.g.  with  or  without  shell  in 
shellfish)  is  subject  to  differing  methods  of  reporting.  The 
seasonality  of  biomass  variations  within  differing  areas  of  the 
Bay  and  the  importance  of  seasonal  advances  or  declines  to  migra¬ 
tion  (import  -  export)  data  is  typically  not  defined. 

Often,  entire  organism  groups  have  either  not  been  well- 
studied,  or  have  only  begun  to  be  studied.  Such  groups  as 
nanno-plankton  and  microzooplankton  have  only  recently  begun 
to  receive  the  attention  that  they  seem  to  deserve  based  on 
their  known  roles  in  the  ecological  system.  It  is  hoped  that 
ongoing  studies  in  these  areas  will  help  clarify  information 
on  distribution  and  importance  of  these  organisms. 

Standard  physiological  data  such  as  respiration  rates,  ingestion 
andexcretion  rates,  mortality  and  natality  are  often  known  only 
in  a  very  abstract  sense,  even  for  major  species.  The  functional 
dependence  of  these  rates  on  temperature,  season,  location,  pre¬ 
sence  of  predators,  or  influence  of  salinity  can  of ter  only  be 
defined  in  a  rough  sense  or  by  applying  data  from  areas  other 
than  the  Chesapeake  Bay. 

In  attempting  any  synthesis  of  the  data  on  Chesapeake  Bay,  one 
is  faced  with  a  patchiness  of  information,  geographically,  tem¬ 
porally,  and  by  species.  The  data  were  usually  taken  by  a  variety 
of  methods  and  at  differing  levels  of  detail.  Much  basic  infor¬ 
mation  has  been  generated  by  degree-granting  institutions  since 
these  bodies  are  somewhat  free  of  the  constraint  of  having  to 
focus  on  potential  environmental  problem  areas.  Other  major 
contributions  have  come  from  largescale,  government  studies 
with  uniform  methodologies  (i.e.  SAV  aerial  photographic  studies). 

The  synthesis  of  this  diverse  material  and  its  application  to 
impact  assessment  is  of  prime  importance  to  federal  agencies 
which  have  baywide  responsibilities.  Practical  problems  such 
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as  the  setting  of  tributary  low  flows  must  be  met.  The  available  dati 
must  therefore  be  synthesized  into  a  workable  methodology  for  low 

flow  assessment.  The  conceptual  methodologies  developed  during 
Phase  I  for  such  impact  assessment  are  sketched  in  the  follow¬ 
ing  section.  These  concepts  have  been  shown  during  this  phase 
of  the  Biota  Assessment  to  provide  the  greatest  adaptability 
to  the  discrepancies  and  gaps  in  the  data  base,  while  allowing 
for  the  incorporation  of  new  information  as  Phase  II  progresses. 


B.  IMPACT  ASSESSMENT  METHODOLOGY 

The  assessment  of  biological  impacts  during  various  low  flow 
scenarios  will  be  accomplished  during  Phase  II  of  the  Chesa¬ 
peake  Bay  Low  Flow  Study.  Impact  assessment  will  be  based  on 
changes  in  habitat  available  to  study  species  under  differing 
salinity  scenarios.  The  salinity  data  for  these  scenarios  will 
be  provided  by  the  Corps  of  Engineers  based  on  tests  on  the 
Chesapeake  Bay  Hydraulic  Model  . 

During  Phase  I,  WESTECH  is  quantifying: 

•  known  habitat 

•  potential  habitat 

for  57  major  Chesapeake  Bay  species.  Known  habitat  is  defined 
in  two  ways  depending  on  the  data  from  which  they  are  derived: 

1.  If  data  were  obtained  from  a  mapped  data  set 
showing  actual  areas ,  these  areas  form  known 
habitat  (i.e.  SAV  maps,  see  Map  Atlas). 

2.  If  the  data  is  derived  from  studies  based  on 
sampling  points ,  we  define  known  habitat  as 
the  smallest  area  meeting  basic  depth  salin¬ 
ity  and  substrate  requirements  which  contain 
those  points  (i.e.  benthic  maps  -  see  Map 
Atlas) . 

Potential  habitats  are  identified  by  the  intersection  of  data 
sets  which  define  multidimensional  habitats  for  each  organism. 
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These  data  sets  include  salinity,  seasonal  occurrence  within 
a  geographic  segment,  depth,  substrate  and  requirements  for 
other  organisms  as  outlined  in  Chapters  IV  and  V.  Total 
potential  habitat  area  includes  both  known  habitat  and  other 
potential  habitat  areas. 

Many  maps  in  the  Map  Atlas  show  known  and  potential  habitat 
combined  due  to  geographical  data  gaps  in  distribution  infor¬ 
mation.  However,  in  all  cases  possible, known  habitat  area  is 
mapped  by  organism.  In  some  cases  this  is  broken  down  by 
organism  lifestage  or  seasonality. 

During  Phase  II,  potential  habitat  areas  of  study  species  will  be 
mapped  under  salinity  scenarios  for  an  average  Water  year,  for  the 
1960  drought,  and  consumptive  water  use  (year  2020) .  The  basis  for 
this  mapping  will  be  the  same  habitat  criteria  developed  in  Phase 
I;  depth,  salinity,  seasonality,  substrate  and  presence  of  other 
required  organisms.  The  degree  of  impact  will  be  quantified  by 
a  comparison  between  habitat  under  average  inflow  conditions  (modal 
hydrograph)  and  habitat  available  to  the  organism  under  each  scenario. 

There  are  several  possible  cases  which  will  determine  the 
methods  of  comparison  to  be  used: 

1.  If  known  habitat  is  poorly  defined  under  base 
conditions ,  or 

2.  If  known  habitat  is  essentially  identical 
with  potential  habitat  under  base  conditions, 
impact  will  be  a  ratio  defined  as  the  quotient 
of  potential  habitats  under  future  scenario  and 
under  base  conditions. 

Impacts  =  potential  habitat  (salinity  scenario  -  X) 
potential  habitat  (base  conditions) 

3.  If  known  habitat  has  been  well  defined  on  a  con¬ 
sistent  basis,  it  is  reasonable  to  assume  that  the 
species  distribution  may  not  completely  fill  the 
potential  habitat. 

(The  reasons  for  this  discrepancy  include  the  fact 
that  we  have  not  completely  defined  all  habitat 


variables, as  well  as  the  fact  that  dispersal  mech¬ 
anisms  and  colonization  rates  affect  the  geographi¬ 
cal  distribution  of  a  species.)  In  this  case,  we 
will  assume  that  the  percentage  of  potential  habi¬ 
tat  occupied  will  remain  the  same  under  base  and 
and  future  conditions.  The  amount  of  the  species 
existing  under  a  particular  salinity  scenario  will 
be  defined  by: 

potential  habitat  (base  conditions) _  =  potential  habitat  (scenario-X) _ 

known  (realized)  habitat  (base  conditions)  projected  realized  habitat  (scenario-X) 

where  in  addition  to  projecting  potential  habitat 
through  use  of  the  ratio  defined  for  points  1  and  2 
above,  a  projection  of  the  currently  realized  habitat 
can  be  made  for  appropriate  species.  This  may  prove 
particularly  useful  for  species  such  as  emergent 
vegetation  where  the  coverage  of  known  habitat  is 
well  known  from  aerial  surveys. 


With  any  of  these  projections,  it  is  important  to  realize  that 
there  is  a  considerable  margin  for  error  in  such  prediction, 
this  margin  or  error  is  due  to  the  fact  that  the  impacts  will 
be  dynamic  and  will  evolve  with  time.  Colonization  or  dispersal 
effects  are  not  taken  into  consideration.  Also,  other  habitat 
variables  have  not  been  taken  into  account  by  the  methodology 
due  to  complexity  of  including  other  variables  and  the  lack 
of  good  data. 

In  Phase  II,  known  and  ptoential  habitat  will  be  measured  on  an  in*4 
dividual  species  basis  for  each  low  flow  scenario  for  each  of  the  13 
geographical  subdivisions  so  that  each  geographical  area  can  be 
handled  as  a  unit  (see  Figure  VI-1).  In  the  analysis  of  salinity 
scenarios,,  species  may,  of  course,  transfer  their  range  from  one 
compartment  ot  another,  or  ranges  may  expand  or  contract.  Where 
possible,  impact  ratios  will  be  calculated  for  each  compartment 
as  well  as  being  aggregated  for  the  whole  Bay . 


Also,  factors  which  may  currently  preclude  a  species  movement 
to  another  area  -  even  one  which  it  may  have  occupied  during 
the  past  -  will  be  addressed,  where  known.  For  example,  the 
dramatic  change  in  SAV  distributions  in  recent  years  may  be 
partially  due  to  increased  runoff  of  toxics  (including  herbicides) 
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uu-ut  mill'll,  these  impact  ratios  are  based  on  imcomplete  data  about 
.mu i sum .  Colonization,  species  dispersal  rates  and  other  factors 
i  ibute  to  differential  use  of  potential  habitat  areas  within  a 
.•ii'ci'ios.  One  cannot,  at  this  time  define  all  habitat  parameters 
or  make  completely  accurate  correlations  even  with  identifiable 
parameters  such  as  substrate.  Thus  the  impact  ratios  may  be  subject 
to  considerable  error  margins . 

As  an  example  of  the  impact  assessment  procedure,  consider  species 
"A"  which  is  a  marine-type  species  and  species  "B"  which  prefers 
polyhaline  salinities,  but  can  withstand  marine  salinities.  Consider 
alos  that  A  preys  upon  B.  The  impact  assessment  procedure  will  then 
consist  of  several  steps.  First  the  habitats  of  both  A  and  B  will 
be  mapped  under  "average  inflow"  conditions  from  hydraulic  model 
data.  Secondly,  the  species  distributions  will  be  mapped  under  one 
of  several  salinity  scenarios  (for  our  purposes  here,  we  will  use 
the  1960’s  drought).  Thirdly,  the  difference  in  habitat  areas  will 
be  measured.  This  constitutes  the  impact  ratios  discussed  in  this 
section.  Lastly,  the  conceptual  or  mathematical  (CBEM)  models  will 
be  applied  to  determine  the  effects  of  species  interaction  as  time 
evolves.  This  information  will  be  used  to  modify  the  interpretation 
of  the  impact  ratios,  which  in  themselves,  represent  only  instan¬ 
taneous,  direct  effects  and  take  no  account  of  trophic  interactions. 

To  reduce  the  possibility  of  readers  using  these  ratios  as  highly 
accurate  impact  predictors,  they  will  be  used  only  to  define 
broad  categories  of  effects.  We  have  not  at  this  time  final¬ 
ized  a  category  scheme,  nor  is  it  appropriate  to  do  so  until 
the  types  and  magnitudes  or  errors  inherent  in  the  ratios  can 
be  estimated  based  on  comparisons  with  historical  data  or  other 
means  (i.e.  1960’s  drought).  Attempts  at  calibration  will  be 
made  using  those  organisms  for  which  drought  data  exists  (only 
available  for  a  few  organisms)  and  through  use  of  the  conceptual 
and  CBEM  models.  This  is  expected  to  generate  an  impact  classi¬ 
fication  similar  to  the  example  below  for  the  Impact  Ratio  (IR) 
for  potential  habitats. 


TABLE  VII- 1.  Example  of  Possible  Primary  Impact  Classification 
Scheme. 


IR  =  (salinity  scenario  -  X) 
(base  conditions) 


potential  habitat 


0.7  -  0.9 


0.9  -  1.1 


1.1  -  1.4 


1.4  + 


Classification 
Severe  nagative  impact 
Moderate  negative  impact 
Low  impact 

Moderate  increase  or  enhancement 
Severe  increase  or  enhancement 
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It  should  be  noted  here  that  if  a  species  shows  an  immediate 
positive  effect,  a  net  gain  for  the  ecosystem  is  not  necessarily 
indicated.  This  effect  may  be  short-lived  in  time  or  canceled 
by  net  drops  in  population  of  other  organisms  due  to  shifts  in 
the  Bay's  trophic  systems.  We  will  attempt,  to  the  extent 
possible,  to  analyze  the  changes  in  direct  impact  of  each  species 
on  overall  ecosystem  function  based  on  known  interactions  from 
the  base  period.  Such  analysis  can  only  be  based  on  conceptual 
or  mathematical  models  of  ecosystem  function. 

The  use  of  some  trophic  model  is  necessary  in  order  to  meet 
the  objectives  of  both  Corps  and  Legislative  decision-making 
which  may  result  from  the  low  flow  study  and  other  portions  of 
the  Chesapeake  Bay  study.  The  use  of  impact  ratios  outlined 
above  gives  a  response  considering  each  organisms  as  an  iso¬ 
lated  entity.  In  fact,  organisms  are  linked  to  each  other  through 
feeding,  competition,  predation  and  other  biological  relation¬ 
ships. 

The  eventual  equilibrium  state  of  the  ecosystem  following  a 
drought  or  after  an  extended  period  of  lower  flows  will  depend 
strongly  on  these  interactions.  For  instance,  if  impact  ratios 
of  oyster  predators  show  increase  of  biomass  (or  habitat)  in  one 
area  of  the  Bay,  oyster  biomass  can  be  expected  to  eventually 
decrease  in  these  same  areas  due  to  the  extended  predator  pene¬ 
tration.  Such  effects  can  only  be  ascertained  by  examining 
potential  habitat  overlaps  between  species  which  affect  each 
other,  either  through  making  judganents  based  on  the  conceptual 
relationships,  or  through  modeling  the  effects  of  the  increase 
in  both  species  as  a  predator-prey  relationship  over  time  with 
a  mathematical  model. 

Due  to  lack  of  data,  some  such  instances  must  be  addressed  at 
a  conceptual  model  level  only.  To  accomplish  this,  the  concep¬ 
tual  model  developed  in  Phase  I  and  supporting  compartment 
models  will  be  used  (see  Chapter  VI).  In  cases  where  the  data 
are  adequate  to  define  basic  species  parameters  and  physical 


driving  functions,  it  will  be  possible  to  apply  the  CBEM 
mathematical  model.  Such  modeling  can  be  used  to  identify 
not  only  gross  probable  outcomes  (impacts)  but  may  in  some 
cases  indicate  the  time  evolution  of  the  system  following  the 
impositions  of  drought  conditions.  Even  in  cases  where  pre¬ 
cise  quantitative  analysis  of  stress  effects  on  the  ecosystem 
are  not  possbile,  the  mathematical  model  can  be  used  to  illus¬ 
trate  the  distribution  of  that  stress  on  ecosystem  components 
by  means  of  sensitivity  analysis  (which  entails  measurement 
only  of  percentage  changes  compared  to  the  change  in  a  par¬ 
ticular  stressed  parameter;  see  Patten  1971) . 


VIII.  CONCLUSIONS 


The  object  of  Phase  I  of  the  Biota  Assessment  has  been  to  syn¬ 
thesize  and  standardize  existing  data  and  studies  on  Chesapeake 
Bay  and  develop  methodologies  to  serve  as  a  basis  for  low  flow  impact 
Assessment.  The  existing  literature ,  published  and  unpublished  ,  is 
present  in  a  diversity  of  forms.  Large  amounts  of  information 
on  physical,  chemical  and  biological  aspects  of  the  estuarine 
system  have  been  collected  over  the  past  2  to  3  decades.  How¬ 
ever,  each  study  is  designed  to  meet  its  own  particular  objec¬ 
tive,  resulting  in  differing  study  methods  and  reporting  for¬ 
mats.  Although  several  person-years  have  been  spent  in  Phase 
I  of  the  Biota  Assessment,  we  have  merely  scratched  the  surface 
of  synthesizing  and  standardizing  some  of  this  information, 
limiting  the  study  to  major  information  concerned  with  certain 
"study  species". 

Knowledge  of  the  physical  dynamics  of  the  Bay  is  based  on  many 
detailed  monitoring  and  theoretical  studies  by  researchers 
from  government  and  academic  institutions.  There  are  still 
major  unanswered  questions  about  overall  Bay  circulation,  and 
large  inconsistancies  in  the  understanding  of  tributaries  and 
local  effects  of  altered  flow  regimes.  Although  the  Environ¬ 
mental  Protection  Agency  and  others  have  been  devoting  consider¬ 
able  effort  to  understanding  this  system,  physical  dynamics  are 
complex  to  monitor  or  model  and  the  present  state-of-the-art 
does  not  permit  accurate  quantification,  particularly  of  many 
subsurface  phenomena  involving  circulation. 

The  chemical  and  nutrient  cycles  are  in  a  somewhat  similar 
situation.  Hare,  though,  some  progress  has  been  made  with  the 
major  nutrients,  and  budget  analysis  has  been  carried  out  on 
some  of  the  major  tributaries.  Some  efforts  have  even  been  made 
to  relate  nutrient  levels  to  flow,  although  such  studies  are 
still  in  early  stages. 
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The  biological  systems  of  Chesapeake  Bay  are  of  similar  complex¬ 
ity.  Here,  however,  we  have  been  able  to  reproduce  major  aspects 
of  structure  and  function  of  the  estuarine  community  through 
analysis  of  the  literature,  habitat  classification  and  limitation 
to  major  species,  followed  by  conceptual  and  mathematical  model¬ 
ing.  It  has  been  found  that  major  trophic  responses  can  be  re¬ 
produced  for  purposes  of  investigation  of  salinity  stress  on 
Bay  organisms. 

The  Bay  literature  is  voluminous.  Many  Bay  organisms  have  been 
heavily  studied;  however,  others  are  known  only  from  taxonomic 
collections,  not  ecological  or  distributional  research.  A 
series  of  trophically  oriented  ecological  studies  of  various 
Bay  areas  with  consistant  methodologies  would  be  the  only  rapid 
method  for  closing  many  of  the  data  gaps  that  exist.  For  cer¬ 
tain  organisms  it  is  possible  to  define  habitat;  however,  this 
is  limited  to  a  few  major  variables  (i.e.  depth,  salinity,  sub¬ 
strate,  etc.).  Even  for  these  variables  there  are  some  inconsis¬ 
tencies  in  measurement  techniques.  More  subtle  effects  and 
variables  such  as  dissolved  oxygen,  nutrients,  etc.  have  been 
studied  only  in  isolated  cases.  From  those  organisms  which  have 
a  data  base  encompassing  distribution,  seasonality  and  ecological 
function,  we  have  mapped  distributions  of  57  "study  species" 
during  a  "base  year"  which  represents  average  flow  and  (presum¬ 
ably)  salinity  data. 

It  has  been  possible,  based  on  assimilation  and  synthesis  of 
major  Bay-related  literature,  to  develop  detailed  conceptual 
models  of  functional  ecological  units  (compartments)  and  syn¬ 
thesize  these  into  a  broad  conceptual  model  which  represents 
major  pathways  of  Bay  energy  flow.  From  this  conceptual  model, 
it  has  been  possible  to  select  and  test  a  mathematical  sub¬ 
model  on  one  of  thirteen  geographical  Bay  units  (segments) . 

This  submodel  is  in  good  agreement  with  published  data  for  this 
area.  The  computerized  mathematical  model  (CBEM)  is  stable  but 
dynamic.  It  is  sensitive  to  salinity  variation  and  should  pro- 
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vide  a  powerful  tool  for  future  work  on  ecosystem  stress. 

Many  data  gaps  exist  with  respect  to  organism  distribution,  stress 
tolerance  ,  ecological  functions,  etc.  However,  enough  of  this 
information  exists  to  construct  a  good  first-cut  at  impact  ana¬ 
lysis  of  flow  differences  using  tools  developed  in  this  phase 
of  the  Biota  Assessment. 

The  first  step  in  such  impact  assessment  will  be  to  map  habi¬ 
tats  for  study  species  under  "normal"  and  various  low  flow  con¬ 
ditions  to  identify  species  reactions.  These  reactions  are 
only  immediate,  primary  reactions  to  the  stresses  imposed  by  low 
flow  conditions.  With  the  passing  of  time,  trophic  effects, 
dispersal,  colonization,  competition  and  seasonal  factors  all 
come  into  play.  While  conceptual  models  can  be  used  to  analyze 
and  predict  these  effects  to  some  extent,  mathematical- computer 
models  such  as  CBEM  add  power  and  extend  the  capabilities  of 
such  analysis.  CBEM  can  be  used  to  analyze  either  short-term, 
immediate  reaction  to  stress  (sensitivity  analysis)  or  track 
the  system  over  time.  In  order  to  accomplish  this,  adjustments 
for  the  physical  and  chemical  characteristics  of  each  tributary 
or  Bay  segment  must  be  made.  The  model  has  been  shown  to  serve 
as  an  effective  tool  to  analyze  effects  of  trophic  interactions 
and  is  in  agreement  with  existing  dynamics  on  the  estuary  tested 
(Patuxent) . 

The  purpose  of  Phase  I  of  the  Biota  Assessment  has  been  to  devel¬ 
op  tools  and  methodologies  for  analysis  of  hydraulic  model  data 
in  Phase  II.  Both  the  tools  and  the  impacts  assessment  methods 
have  been  developed  and  (at  least  theoretically)  shown  to  be 
effective  and  workable.  During  Phase  II  of  the  Biota  Assessment 
the  tools  will  be  applied  to  scenarios  involving  differing  flow  and 
salinity  regimes  as  a  part  of  the  Corps  of  Engineers  overall 
Low  Flow  Project  for  the  Chesapeake  Bay. 


END 

/ 

* 

FILMED 


